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Abstract 

Terahertz  (THz)  time  domain  spectroscopy  (TDS)  was  assessed  as  a 
nondestructive  evaluation  technique  for  aircraft  composites.  Material  properties  of  glass 
fiber  composite  were  measured  using  both  transmission  and  reflection  configuration.  The 
interaction  of  THz  with  a  glass  fiber  composite  was  then  analyzed,  including  the  effects 
of  scattering,  absorption,  and  the  index  of  refraction,  as  well  as  effective  medium 
approximations.  THz  TDS,  in  both  transmission  and  reflection  configuration,  was  used 
to  study  composite  damage,  including  voids,  delaminations,  mechanical  damage,  and  heat 
damage.  Measurement  of  the  material  properties  on  samples  with  localized  heat  damage 
showed  that  burning  did  not  change  the  refractive  index  or  absorption  coefficient 
noticeably;  however,  material  blistering  was  detected.  Voids  were  located  by  THz  TDS 
transmission  and  reflection  imaging  using  amplitude  and  phase  techniques.  The  depth  of 
delaminations  was  measured  via  the  timing  of  Fabry-Perot  reflections  after  the  main 
pulse.  Evidence  of  bending  stress  damage  and  simulated  hidden  cracks  was  also  detected 
with  terahertz  imaging. 
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Nondestructive  Evaluation  of  Aircraft  Composites  Using 
Terahertz  Time  Domain  Spectroscopy 

1.  Introduction 

Terahertz  (THz)  radiation  can  be  defined  as  the  region  of  the  electromagnetic 
spectrum  with  frequencies  between  0.1  THz  and  10  THz  or  with  wavelengths  between  3 
mm  and  30  pm.  The  THz  “gap”,  as  shown  in  Figure  1,  falls  between  the  microwave 
region  and  the  infrared/optical  region.  Microwave  electronic  devices  are  typically 
comparable  in  size  to  the  wavelength  of  the  radiation  and  support  single-mode  or  low- 
order  mode  guided  waves.  In  contrast,  optical  and  infrared  devices  have  dimensions 
much  larger  than  a  wavelength  and  often  support  beams  containing  many  modes.  The 
merging  of  laser  and  optics  technology  with  micro  chip  electronics  technology  has  been 
instrumental  in  the  development  of  THz  devices  and  the  advancement  of  THz  research. 
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Figure  1.  Electromagnetic  Spectrum  [14]. 
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1.1  Background 


With  the  rapid  advances  in  the  development  of  laser  technology,  nonlinear  optics 
and  optoelectronic  techniques  for  the  production  and  measurement  of  terahertz  radiation 
have  proliferated.  One  of  the  earliest  and  most  promising  of  these  techniques  is 
terahertz  time-domain  spectroscopy  (THz  TDS),  a  method  which  relies  on  the  generation 
of  broadband  electromagnetic  transients  using  ultrafast  laser  pulses.  Several  novel 
measurement  capabilities  have  become  practical,  including  time-resolved  techniques  like 
the  optical  pump  -  THz  probe  measurement,  THz  emission  spectroscopy,  and  THz 
correlation  spectroscopy.  The  first  reported  use  of  THz  TDS  for  imaging  was  reported  in 
1995  and  this  work  has  sparked  a  flood  of  interest  and  suggested  a  range  of  new 
applications  for  THz  “T-ray”  imaging  [48],  It  has  even  received  media  attention  recently 
in  the  novel  Games  of  State,  by  Tom  Clancy  [15]. 

Composite  materials  such  as  fiberglass,  Kevlar,  and  carbon  fiber  have  high 
strength  to  weight  ratios  and  are  increasingly  being  used  as  structural  components  in  high 
performance  military  aircraft  [89],  The  navy’s  F/A-18A/B  and  the  F/A-18E/F  structural 
weights  are  composed  of  10%  and  19%  composites,  and  approximately  24%  of  the  F- 
22’s  structural  weight  is  composite  in  origin.  The  F-35  Joint  Strike  Fighter  will  also  have 
a  high  percentage  of  its  weight  made  up  of  composites. 

Fiber  composites  not  only  impart  increased  strength  to  aircraft  structures,  but  also 
improve  aerodynamic  performance,  increase  safety,  and  reduce  corrosion.  However, 
composites  can  be  weakened  by  various  defects  that  are  created  during  either  fabrication 
or  during  the  lifecycle  of  the  aircraft.  In  the  layer  build-up  stage  of  the  manufacturing 
process,  foreign  debris  can  become  embedded  within  the  composite  layers,  and  voids  can 


2 


be  formed  by  air  bubbles  or  large  gaps  between  the  fiber  weave.  Defects  that  are  caused 
by  foreign  debris  or  voids  decrease  the  layer  adhesion  and  the  overall  strength  of  a 
composite  component.  Current  quality  control  methods  are  incapable  of  finding  these 
defects  until  curing  has  been  completed,  requiring  that  the  component  be  reworked  or 
thrown  out  and  increasing  the  cost  of  part  production.  Additionally,  routine  maintenance 
of  in-service  composite  materials  requires  the  application  of  relatively  complex 
inspection  and  repair  techniques,  which  often  have  to  be  performed  in  the  field. 

New  techniques  are  being  explored  to  improve  the  efficiency  of  composite  part 
inspection  during  construction  and  repair.  THz  radiation  offers  potential  efficiency 
improvements,  because  it  has  the  unique  ability  to  penetrate  composites  and  identify 
foreign  debris  and  other  defects  such  as  delaminations,  air  bubbles,  and  heat  damage. 

THz  would  offer  a  non-invasive,  non-contact,  non-ionizing  way  of  assessing  composite 
part  condition  and  could  help  overcome  the  shortcomings  of  other  nondestructive 
evaluation  (NDE)  techniques  such  as  x-rays,  ultrasound,  video  inspection,  eddy  currents, 
and  thermographic  techniques. 

1.2  Problem  Statement 

Significant  progress  has  been  made  recently  in  the  development  of  the  sources 
and  technology  for  THz  imaging,  sensing,  and  spectroscopic  and  communications 
applications,  but  the  technology  still  remains  limited  compared  to  that  used  for  the 
microwave  and  optical  regions  of  the  spectrum.  THz  technology  is  being  used  in  an 
increasingly  wide  variety  of  applications  including:  information  and  communications 
technology;  biology  and  medical  sciences;  homeland  security;  quality  control  of  food  and 
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agricultural  products;  global  environment  monitoring;  ultrafast  computing;  and  NDE 
among  others.  Effective  NDE  can  be  performed  with  a  THz  imaging  system;  however,  a 
THz  imaging  system  demands  a  high  repetition  scan  rate  and  a  high  average  power  for 
favorable  signal-to-noise  ratios.  To  fully  take  advantage  of  this  emerging  technology  will 
require  more  compact,  efficient,  and  higher  power  sources.  This  research  effort  attempts 
to  take  advantage  of  commercially  available  THz  sources  in  an  effort  to  more  effectively 
characterize  and  image  composite  materials. 

Innovative  solutions  for  quality  control  of  composite  materials  during  the 
manufacturing  process  and  subsequent  routine  maintenance  are  needed  by  the  Air  Force 
to  save  time  and  lower  cost.  As  a  NDE  technique,  THz  imaging  has  the  potential  to  be  an 
effective  method  of  quality  control  and  maintenance  inspection  of  composites.  Although 
techniques  have  been  developed  to  measure  material  properties  at  THz  frequencies,  work 
still  needs  to  be  done  to  understand  and  characterize  the  electromagnetic  interaction  of 
THz  with  various  materials.  The  optical  and  electrical  properties  of  composites  need  to 
be  fundamentally  understood  to  improve  imaging  techniques  for  NDE.  Additionally,  the 
capability  of  THz  TDS  in  both  transmission  and  reflection  configuration  should  be  tested 
on  representative  damage  samples  for  effectiveness. 

1.3  Research  Objectives 

The  objective  of  this  dissertation  is  to  evaluate  THz  TDS  for  nondestructive 
testing  and  evaluation  of  aircraft  composites  by  using  theory,  modeling,  and 
experimentation  to  analyze  the  interaction  of  THz  radiation  with  composite  materials. 

An  additional  objective  is  to  show  that  THz  imaging  shows  promise  for  NDE,  especially 
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when  imaged  with  a  high  average  power,  high  repetition  rate  source  and  a  fast  scanning 
and  collection  technique  for  inspecting  the  sample. 

1.4  Experimental  Approach 

The  first  step  of  this  research  was  to  set  up  a  THz  TDS  system  in  transmission 
configuration  using  a  photoconductive  switch  for  the  transmitter  and  an  electro-optic 
crystal  for  the  receiver.  Using  this  setup  initially,  data  was  collected  to  make  estimates  of 
optical  and  electrical  parameters  for  semiconductors  and  various  composite  materials, 
such  as  space  shuttle  foam  insulation  and  aircraft  glass  fiber  composites.  The  composite 
material  parameters  could  then  be  used,  along  with  effective  medium  theory  and 
scattering  theory,  to  describe  the  interaction  of  the  THz  radiation  with  the  sample.  Then 
a  series  of  THz  images,  formed  with  a  raster  scan  pattern,  and  depth  measurements  were 
collected  on  samples  prepared  with  various  forms  of  representative  damage. 

The  next  step  in  the  experiment  was  to  change  the  THz  TDS  system  to  a 
reflection  configuration.  With  this  setup,  the  material  parameter  measurements  were 
measured  on  the  same  samples  used  in  the  transmission  experiments.  Finally,  reflection 
mode  images  and  depth  measurements  were  taken  on  the  same  damage  samples  for  a 
comparison  of  the  two  scan  methods. 

1.5  Assumptions/  Limitations 

The  variety  of  aircraft  composite  samples  available  for  this  research  was 
somewhat  limited,  but  encompassed  a  representative  set  of  damages  that  aircraft 
composites  could  encounter.  Additionally,  the  ability  to  obtain  the  limited  set  of 
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representative  aircraft  damage  samples  allowed  this  research  to  be  unique  and  interesting. 
The  research  was  restricted  to  glass  fiber  composites,  however,  because  of  the  inability  of 
THz  radiation  to  penetrate  other  composites  like  carbon  fiber.  It  was  also  difficult  and 
expensive  to  obtain  semiconductor  samples  with  low  resistivities,  due  to  the 
manufacturers’  reluctance  to  sell  wafers  in  small  quantities.  The  assumption  was  made 
that  a  high  resistivity  wafer  was  sufficient  to  calibrate  a  THz  TDS  system.  Another  of 
our  limitations  was  that  we  had  limited  knowledge  about  the  properties,  structure,  and 
volume  concentration  of  the  composite  constituents.  This  resulted  in  making 
assumptions  for  the  effective  medium  and  scattering  approximations. 

The  lab  setup  that  was  used  was  limited  in  its  ability  to  collect  THz  TDS  data  at  a 
fast  enough  rate  to  be  considered  operationally  useful.  The  assumption  was  made  that  the 
data  we  collected  to  make  material  parameter  measurements  and  to  form  images  was  for 
a  proof  of  concept  demonstration  and  that  other  organizations  are  developing  the  means 
to  rapidly  scan  large  samples  in  short  amounts  of  time.  Commercially  available  THz, 
optical,  and  mechanical  components  were  sufficient  to  conduct  this  research,  but 
customized  components  and  techniques  could  allow  the  THz  beam  to  be  more  tightly 
focused  onto  the  sample  and  scanned  at  a  faster  rate  to  maximize  defect  detection. 

1.6  Overview 

The  following  chapters  of  this  dissertation  outline  the  theory,  literature  review, 
experimental  setups,  data  collection  and  analysis,  and  imaging  of  composite  materials 
with  THz  TDS  in  various  configurations.  Chapter  2  begins  by  describing  the  theory  for 
the  emitter  and  detector  technology  used  in  this  research,  followed  by  a  description  of  the 
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components  of  the  THz  TDS  system.  The  next  section  describes  data  analysis,  including 
material  parameter  estimation,  scattering,  and  imaging.  The  chapter  closes  with  a 
summary  of  some  of  the  recent  findings  in  the  literature  regarding  the  use  of  THz  for 
NDE. 

Chapter  3  describes  the  data  collection,  analysis,  and  results  of  the  material 
parameter  estimation  of  the  composite  materials  studied  in  this  dissertation.  The  first 
section  explains  the  setup  and  characterization  of  the  THz  TDS  data  collection  system. 
This  is  followed  by  the  analysis  of  material  parameters  collected  using  the  transmission 
mode,  including  a  comparison  with  the  Drude  Model,  Fabry-Perot  cancellation 
techniques,  and  effective  medium  and  scattering  approximations.  The  next  section  is 
dedicated  to  making  the  same  material  parameter  calculations  using  THz  TDS  in 
reflective  mode.  The  last  section  concludes  with  an  analysis  of  THz  TDS  material 
parameter  measurements  on  burn  damage  as  a  potential  NDE  technique. 

In  chapter  4,  the  imaging  and  analysis  of  composite  damage  is  the  focus,  using 
both  transmission  and  reflection  configurations.  The  first  section  explains  the 
optimization  of  the  THz  optical  system  design  for  a  minimized  spot  size  and  the  image 
processing  technique  used  during  post  processing.  The  next  section  shows  the  results  of 
transmission  mode  imaging  of  the  damage  samples  and  a  depth  of  delamination  analysis. 
The  last  section  of  chapter  4  shows  the  results  of  2-D  imaging  and  depth  analysis  on  the 
same  samples  using  the  reflection  setup.  The  last  chapter  shows  a  summary  of  the  results 
and  recommendations  for  future  THz  research. 
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2.  Theory  and  Background 


In  order  to  estimate  the  optical  and  electrical  parameters  of  a  material  using  THz 
TDS  or  image  that  material  using  THz  TDS,  it  is  necessary  to  establish  the  theoretical 
background  on  the  generation  and  detection  of  THz  radiation  as  well  as  measurement 
techniques.  This  chapter  lays  the  foundation  required  to  explore  the  properties  of 
composite  materials  at  THz  frequencies  and  the  insight  to  investigate  THz  TDS  as  a  NDE 
technique. 

2.1  THz  Pulse  Generation 

In  recent  years  scientists  have  developed  numerous  techniques  for  generating  THz 
radiation  in  an  effort  to  explore  this  spectral  region.  This  research  dealt  solely  with 
pulsed  THz  sources  instead  of  continuous  wave  sources.  A  traditional,  transient  current 
technique  using  a  photoconductive  switch  was  used  exclusively  as  the  THz  emitter. 
Optical  rectification,  a  traditional  nonlinear  technique,  was  also  demonstrated  in  the  lab 
as  a  means  to  generate  THz,  but  was  not  used  in  data  collections  for  this  research. 
Currently,  the  average  power  generated  from  a  typical  photoconductive  switch  is  less  1 
pW,  while  the  peak  power  from  optical  rectification  may  reach  up  to  1  pW  [81].  Larger 
area  photoconductive  switching  devices  have  yielded  electric  fields  as  high  150  kV/cm 
and  a  new  technique,  using  a  four- wave  mixing  process,  showed  a  record  high  field 
strength  of  400  kV/cm  [7,  88],  A  recent  publication  demonstrated  a  high  average  power 
of  0.25  mW  from  a  single  cycle  THz  pulse  using  an  optical  rectification  technique  with  a 
lithium  niobate  crystal  [31]. 
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2.1.1  Photoconductive  Antennas 


Photoconductive  (PC)  switches  are  composed  of  metallic  striplines  deposited  on  a 
semiconductor  material.  PC  switches,  also  called  Auston  switches,  are  based  on  the 
pioneered  work  by  Auston  in  the  mid-1970s,  using  a  mode-locked  Nd:glass  laser  to  pump 
high-resistivity  silicon  [48],  When  an  optical  pulse,  such  as  a  femtosecond  laser,  strikes 
the  semiconductor  at  a  wavelength  above  the  bandgap,  electron  (holes)  are  generated  in 
the  conduction  (valence)  band.  The  carriers  are  then  accelerated  in  phase  by  the  bias 
field  toward  the  anode,  resulting  in  a  pulsed  photocurrent  in  the  PC  antenna.  Time- 
varying  current  occurs  in  the  subpicosecond  regime  and  thus  emits  a  subpicosecond 
electromagnetic  transient  or  THz  pulse  [48],  A  diagram  of  a  photoconductive  switch  is 
shown  in  Figure  2. 

C  op  lanar  Hertzian  dip  ole 


RadiatedField:E(r,t)  #  photoconductive 

Side  View  Top  View  substrate 


Figure  2.  Photoconductive  switch  [48]. 

For  an  elementary  dipole  antenna  in  free  space,  the  radiated  electric  field  E(r,t ) 
at  a  distance  r  and  the  time  t  are  described  by 


E(r,t) 


h  m) 

\nsnc2r  dt 


sin# 


(1) 
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where  J(t )  is  the  current  in  the  dipole,  le  is  the  effective  length  of  the  dipole,  s0  is  the 

dielectric  constant  of  free  space,  c  is  the  speed  of  light  in  a  vacuum,  and  6  is  the  angle 
from  the  direction  of  the  dipole.  The  THz  radiation  amplitude  is  proportional  to  the 
length  of  the  dipole  and  the  time  derivative  of  the  transient  photocurrent.  The 
photocurrent  density  is  described  by 

j(t)ccl(t)®[n(t)qv(t)]  (2) 

where  ®  denotes  the  convolution  product,  1(f)  is  the  optical  intensity  profile,  and  q , 
n(t) ,  and  v(t)  are  the  charge,  density,  and  velocity  of  photocarriers.  The  dynamics  of 
photogenerated  free  carriers  in  a  semiconductor  is  well  described  by  the  classical  Drude 
model.  According  to  the  Drude  model,  the  average  velocity  of  free  carriers  follows  the 
differential  equation 

^l=-vM+±m  (3) 

dt  r  m 

where  z  is  the  momentum  relaxation  time  and  m  is  the  effective  mass  of  the  carrier. 

The  current  density  is  represented  by  n(t)qv(t )  which  is  the  impulse  response  of  the  PC 
antenna.  Most  of  the  THz  pulse  is  emitted  from  the  substrate  side.  This  is  based  on  the 
antenna  theory  that  shows  that  a  dipole  antenna  on  the  surface  of  a  dielectric  material 
emits  roughly  s  12  times  more  power  to  the  dielectric  material  than  to  the  air,  where  s 
is  the  relative  dielectric  constant  of  the  substrate  [67]. 

In  a  PC  antenna,  properties  such  as  a  high  electron  mobility  and  a  high  breakdown 
voltage  should  be  maximized.  Using  a  semiconductor  source  with  a  short  carrier  lifetime 
is  not  a  requirement  to  generate  THz.  Materials  are  desired  which  have  a  high  resistivity 
because  a  high  voltage  is  applied  across  the  striplines  and  a  low  current  is  important  to 
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minimize  the  heat.  Faster  electron  mobilities  lead  to  faster  accelerations  and  higher  THz 
fields.  Two  of  the  most  often  used  substances  are  low  temperature  grown  gallium 
arsenide  (LT-GaAs)  and  radiation-damaged  silicon-on-sapphire  (RD-SOS).  The  PC 
switch  design  used  for  this  research  is  LT-GaAs,  which  has  a  high  electron  mobility,  high 
resistivity,  and  high  breakdown  voltage.  The  amount  of  excess  arsenic  in  the  GaAs  can 
be  controlled  during  epitaxial  growth  to  regulate  the  carrier  lifetimes  (0.3  -100  ps),  and 
longer  carrier  lifetimes  result  from  using  a  higher  temperature  during  growth  [67]. 

2.1.2  Optical  Rectification  in  Nonlinear  Media 

Optical  rectification  was  one  of  the  first  techniques  used  to  generate  THz  and  has 
the  advantage  over  photoconductive  switches  in  that  its  output  saturates  at  much  higher 
pump  powers.  Optical  rectification  is  a  nonlinear  optical  effect  caused  by  the  second 
order  susceptibility  % 12 1  that  produces  a  polarization  Pit)  in  the  material: 

Pit)  =  S0 X(l)m  +  So X(1)mE*  it)  +  s0 %(3)E3  it)  + .  .  (4) 

where  e0  is  the  dielectric  constant  in  free  space  and  E(t)  is  the  electric  field  of  the 

incident  light  [86],  The  electric  field  can  be  represented  by  an  amplitude  and  phase 
component  and  its  complex  conjugate  where  co  is  the  frequency  of  the  light: 

Eit)  =  E0ekot  +c.c.  (5) 

This  second  order  effect  causes  the  creation  of  a  second  harmonic  component  and 
a  DC  component.  If  the  nonlinear  material  is  not  phase-matched  for  the  second  harmonic 
then  the  second  harmonic  output  intensity  is  low.  Since  the  pump  input  is  pulsed,  the  DC 
component  will  have  a  broad  spectral  output  which  is  proportional  to  the  inverse  of  its 
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pulse  width.  For  a  femtosecond  laser  this  falls  in  the  THz  range  [67].  An  optical 
rectification  diagram  is  shown  in  Figure  3. 


Figure  3.  Diagram  of  optical  rectification. 

The  nonlinear  material  used  for  optical  rectification  was  zinc  telluride  (ZnTe). 
ZnTe  is  phase  matched  at  the  frequency  of  the  Ti:Sapphire  laser  and  at  THz  frequencies; 
however,  it  experiences  strong  absorption  from  transverse  optical  phonons  at  around  5.3 
THz  and  1 1  THz.  Optical  rectification  offers  a  potentially  larger  emitted  radiation 
bandwidth  than  the  photoconductive  switch.  With  the  optical  rectification  technique,  the 
generated  THz  radiation  is  proportional  to  the  second-order  time  derivative  of  the 
nonlinear  polarization  P(t),  and  therefore  is  also  proportional  to  the  second-order  time 
derivative  of  the  pump  intensity  I(t).  A  Fourier  transform  gives  the  relationship  between 
the  bandwidth  of  the  THz  radiation  or  and  frequency  distribution  of  the  laser  pulse  I( co) 

ENL  (<u)oc  <u2G(oi)P(oi)  oc  col(co)  (6) 

where  G(co)  is  the  phase  matching  condition.  In  the  third  term  of  the  equation,  G(a>)  is 
approximated  to  be  proportional  to  1/co  in  the  high  frequency  region. 

On  the  other  hand,  THz  generated  from  a  PC  switch  is  proportional  to  the  first 
order  time  derivative  of  the  transient  current  J(t)  and  photocarrier  density  N(t) 
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After  taking  a  Fourier  transform  of  the  pump  intensity  term  I(t)  in  equation  7,  a 
comparison  with  equation  6  shows  that  the  radiation  bandwidth  of  optical  rectification  is 
proportionally  larger  than  a  PC  switch  by  the  laser  pump  frequency  co  [67]. 

2.2  THz  Pulse  Detection 

The  two  traditional  methods  of  THz  pulse  detection  that  were  considered  for  this 
research  were  the  photoconductive  switch  and  the  electro-optic  (EO)  sampling  technique. 
For  this  research,  we  initially  tried  using  LT-GaAs  for  the  THz  detector,  but  then  used  a 
ZnTe  (110)  crystal  as  an  EO  detector  when  the  LT-GaAs  detector  malfunctioned. 

2.2.1  Photoconductive  Antennas 

The  detection  principle  of  a  PC  antenna  is  the  inverse  process  of  the  THz  pulse 
generation  process.  In  the  detection  of  THz  radiation  with  PC  antennas,  a  current  meter 
is  connected  instead  of  a  bias  voltage.  Photocarriers  created  by  a  probe  pulse  are 
accelerated  by  the  electric  field  of  the  incident  THz  radiation  and  are  detected  as  current. 
This  temporal  change  of  the  input  pulse  can  be  traced  by  changing  the  arrival  time  of  the 
optical  gate  pulse. 

The  photocurrent  J(t)  of  the  incident  radiation  at  a  time  delay  t  is  described  by  the 
following  equation: 

/•oo 

J{t)  =  e/j\  E(t')N(t'-t)dt'  (8) 

J—  oo 
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where  Eft  ’),  Nft  ’),  e,  and  ji  are  the  incident  electric  field  of  the  THz  radiation,  number  of 
photocarriers  created,  elementary  electric  charge,  and  electron  mobility,  respectively 
[67].  When  the  PC  antenna  is  considered  as  a  sampling  detector,  the  temporal  increase 
and  decrease  of  Nft)  should  be  as  short  as  possible.  Jft)  would  directly  follow  Eft),  as  the 
delay  time  t  of  the  optical  gate  pulse  were  being  scanned,  if  Nft  ’-t)  were  a  delta  function. 
Nft)  can  be  restricted  by  numerous  factors,  such  as  the  gating  pulse  width  and  the  carrier 
lifetime  and  momentum  relaxation  of  photocarriers.  This  is  one  of  the  reasons  for 
selecting  semiconductor  materials  with  short  carrier  lifetimes.  The  equation  for  the 
photocurrent  can  be  transformed  into  the  frequency  domain  according  to  the  convolution 
theorem  of  the  Fourier  Transform: 

J(o>)*N(co)EIm(<o)  (9) 

The  terms  Jfco),  N(co),  and  EThz(co)  are  Fourier  Transforms  of  each  of  the  components 
J(t),  N(t),  and  Eft)  respectively. 

2.2.2  Electro-Optic  Sampling  in  Crystals 

The  coherent  detection  of  a  THz  pulse  with  EO  crystals  is  based  on  the  linear  EO 
effect,  also  known  as  the  Pockel’s  Effect.  The  incident  THz  pulse  modifies  the  refractive 
index  ellipsoid  of  the  EO  crystal  creating  a  phase  retardation  of  the  linearly  polarized 
optical  probe  beam.  The  field  strength  of  the  THz  pulse  can  be  detected  by  monitoring 
the  phase  retardation  of  the  probe  beam.  A  typical  setup  for  EO  detection  of  THz  pulses 
is  shown  in  Figure  4. 
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A  pellicle  beam  splitter  combines  the  THz  beam  and  the  probe  beam  so  that  the 
polarizations  of  the  THz  and  probe  beam  are  aligned  parallel  to  the  (110)  direction  of  the 
ZnTe  sensor  crystal.  The  phase  retardation  induced  in  the  EO  crystal  is  given  by: 

0  jr 

Ar=Tdn°r,r"Em  (10) 


Quarter 


Figure  4.  Setup  for  EO  detection  of  THz  pulses. 

where  d  is  the  EO  crystal  thickness,  nop t  is  the  group  refractive  index  of  the  EO  crystal  at 
the  wavelength  of  the  probe  beam,  and  r4\  is  the  EO  coefficient  [67].  After  the  sensor 
crystal,  a  quarter-wave  plate  is  used  to  apply  a  jt/4  optical  bias  to  the  probe  beam.  A 
Wollaston  polarizer  is  then  used  to  convert  the  THz-radiation-field-induced  phase 
retardation  of  the  probe  beam  into  an  intensity  modulation  between  the  two  mutually 
orthogonal,  linearly  polarized  beams.  A  pair  of  balanced  silicon  p-i-n  photodiodes  is  fed 
to  a  lock-in  amplifier  referenced  to  the  chopping  frequency,  which  has  been  generated  by 
turning  the  emitter  on  and  off  at  a  high  frequency  (—50  kHz).  This  technique  is  called  EO 
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sampling  because  of  the  sampling  behavior  exhibited  by  the  instantaneous  response  of  the 
Pockel’s  Effect  [67]. 

2.3  Terahertz  Time  Domain  Spectroscopy  System 

A  THz  TDS  system  is  a  coherent  emission  and  detection  system  that  emits  single¬ 
cycle  THz  pulses  and  detects  them  at  repetition  rates  on  the  order  of  100  MHz.  The 
signal  is  detected  in  the  form  of  an  electric  field  and  the  Fourier  transform  of  the  pulse 
signal  results  in  both  amplitude  and  phase  spectra  over  a  wide  spectral  range.  It  has  been 
recognized  that  THz  TDS  is  advantageous  over  Fourier  Transform  Spectroscopy  (FTS), 
since  it  gives  both  phase  and  amplitude  information  and  avoids  the  uncertainties  caused 
by  the  Kramers-Kronig  analysis.  Additionally,  THz  TDS  can  operate  with  a  higher 
signal-to-noise  ratio  than  FTS,  with  better  time  resolution,  and  with  a  lower  minimum 
detectable  power  [67], 

The  first  use  of  THz  TDS  for  imaging  was  reported  by  Hu  and  Nuss  in  1995  [33], 
THz  TDS,  which  can  be  used  to  probe  samples  for  NDE,  can  be  scanned  in  two 
dimensions,  in  either  reflective  or  transmissive  mode.  If  one  adds  time  information  with 
full  field  measurements  of  both  phase  and  amplitude,  then  samples  can  be  scanned  in 
three  dimensions,  or  what  is  called  THz  tomography.  In  the  next  sections,  the 
components  of  a  typical  THz  TDS  imaging  system  will  be  described. 

A  diagram  of  a  typical  setup  for  THz  TDS  in  transmission  mode  using 
photoconductive  switches  is  shown  in  Figure  5.  A  THz  TDS  system  consists  of  a 
femtosecond  laser,  an  optical  delay  line,  an  optically  gated  THz  transmitter,  optics  for 
collimating  and  focusing  the  THz  beam,  the  sample  to  be  imaged,  an  optically  gated  THz 
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receiver,  and  a  lock-in  amplifier.  Various  laser  optics,  including  thin  lenses  and  highly 
reflective  metallic-coated  mirrors,  were  chosen  to  minimize  dispersion  effects. 

The  diagram  in  Figure  5  shows  that  PC  switches  can  be  used  for  both  the  THz 
transmitter  and  receiver.  One  could  also  use  a  ZnTe  crystal  for  optical  rectification  for 
the  THz  transmitter  and  the  EO  sampling  technique  from  Figure  4  for  the  receiver.  The 
technique  of  chopping  the  pump  beam  increases  the  signal  to  noise  ratio  of  the  output, 


but  also  is  a  source  of  additional  loss  to  the  pump  beam.  For  this  experiment  we  have 
chosen  to  use  a  switching  power  supply  (-50  kHz)  for  the  THz  emitter,  rather  than 
chopping  the  pump  beam,  so  that  we  can  operate  at  a  higher  frequency  and  reduce  the 
noise  that  is  experienced  at  lower  frequencies. 
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2.3.1  Femtosecond  Laser 


A  Coherent  Mira  Model  900-F,  Ti:  Sapphire,  mode-locked  laser  that  operates 
around  800  nm  wavelengths  was  used  for  this  research.  It  is  optically  pumped  with  a 
diode-pumped,  frequency-doubled  Nd:YV04  laser.  The  Ti:Sapphire  laser  operates  with  a 
pulse  repetition  rate  of  76  MHz  and  has  an  output  power  of  over  500  mW.  The  pulse 
width  of  the  Ti:  Sapphire  laser  was  measured  with  a  Pulse  Check  autocorrelator  and 
autocorrelation  widths  ranged  between  150-  170  femtoseconds.  Since  Mira  pulses  are 
best  described  by  a  sech2  function,  a  factor  of  0.648  should  be  applied  to  convert 
observed  autocorrelation  widths  to  actual  pulse  widths.  This  converts  the  observed  pulse 
width  to  between  97-110  femtoseconds. 

Dispersion  of  the  Ti:Sapphire  laser  beam  was  also  a  setup  consideration, 
especially  when  using  femtosecond  pulses.  Lenses  had  to  be  chosen  that  were  both  thin 
and  had  material  that  provided  the  minimum  amount  of  dispersion.  Mirrors  had  to  have  a 
reflective  coating  which  had  the  maximum  amount  of  reflectivity  and  the  minimum 
amount  of  dispersion. 

2.3.2  THz  Beam  Optics 

A  set  of  THz  beam  optics  was  used  to  collimate  the  output  from  the  emitter,  focus 
the  THz  on  to  the  sample,  and  then  focus  the  THz  on  to  the  detector.  A  key  ingredient  to 
a  high  performance  THz  TDS  system  is  an  optical  system  that  allows  one  to  focus  the 
THz  waves  to  a  diffraction-limited  focal  point  at  the  object  while  allowing  the  highest 
possible  transmission.  At  THz  frequencies,  the  wavelength  is  not  negligible  compared 
with  the  size  of  the  optical  elements,  and  diffraction  effects  can  dominate  ray 
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propagation.  Additionally,  because  of  the  large  spectral  bandwidth  of  THz  emitters,  the 
optical  system  needs  to  be  achromatic  and  exhibit  a  flat  phase  response  over  the 
frequency  range  of  the  pulse. 

2.3.2. 1  Substrate  Lens 

A  hyperhemispherical  substrate  lens  design  made  of  silicon  was  used  to  direct  the 
THz  radiation  to  or  from  the  substrate.  Silicon  has  practically  no  dispersion  over  the 
entire  THz  frequency  range  and  therefore  is  achromatic.  There  are  two  types  of 
hyperhemispherical  designs  that  are  typically  used  for  the  substrate  lenses:  collimating 
and  aplanatic.  In  a  collimating  lens,  the  rays  emitted  near  the  optic  axis  emerge  as  a 
collimated  beam,  while  the  rays  emitted  at  larger  angles  emerge  at  substantial  angles  or 
are  internally  reflected  and  lost.  The  aplanatic  hyperhemispherical  lens  design,  which 
ensures  that  the  critical  angle  for  emerging  rays  is  at  90  degrees,  causes  rays  to  diverge  as 
a  Gaussian  beam  with  a  half-angle  of  sin'1  (1/n)  (~  15°  for  silicon).  It  has  three 
advantages  over  the  collimated  lens.  It  exhibits  no  astigmatism,  it  suffers  from  no  losses 
due  to  internal  reflection,  and  it  allows  for  a  larger  effective  aperture  for  the  lens  and  no 
additional  diffraction  of  the  emerging  beam  [48], 

For  this  research  project,  only  the  collimating  lens  was  available,  since  it  was 
manufactured  as  one  assembly  with  the  Ekspla  photoconductive  antennas,  although  the 
output  beam  wasn’t  really  collimated  in  that  it  initially  converged  to  a  minimum  spot  and 
then  diverged.  Either  one  of  the  two  substrate  lens  designs  can  be  employed  effectively, 
and  can  be  useful  depending  on  the  type  of  application.  The  collimating  lens  can  be  more 
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useful  in  imaging  applications,  since  the  focal  spot  of  a  THz  TDS  system  can  be 
frequency-independent  when  re-imaged  through  a  second  lens  or  parabolic  mirror. 


23.2.2  Off-Axis  Parabolic  Mirrors 


A  series  of  off-axis  parabolic  mirrors  can  be  used  to  collimate  and/or  focus  the 
THz  beam  once  it  departs  the  silicon  lens.  For  this  research,  gold  coated  were  chosen 
initially  to  maximize  reflectivity,  but  aluminum  also  worked  effectively  later  in  the 
research.  THz  beams  can  be  considered  Gaussian  in  nature  and  can  therefore  be 
numerically  modeled  as  Gaussian  beams.  The  effective  diameter  of  a  focused  Gaussian 
beam  can  be  defined  in  terms  of  the  diameter  which  contains  86%  of  the  focused  energy 
and  at  the  edges  of  which  the  focused  intensity  is  already  down  to  1/e2  of  its  peak  value. 
Based  on  this  criteria,  the  effective  diameter  of  the  focused  Gaussian  spot,  d,  is  given  by 


d 


(11) 


where  X  is  the  wavelength,  /  is  the  focal  length  of  the  focusing  element,  and  d0  is  the 

diameter  of  the  collimated  beam  [72],  This  diameter  is  an  approximation  of  the 
diffraction  limit  of  the  THz  imaging  system.  Corrections  must  be  applied  if  the  THz 
beam  incident  on  the  focusing  mirror  is  not  collimated. 


2.3.23  Gaussian  Beams 

Terahertz  radiation  can  be  modeled  as  a  Gaussian  beam,  with  a  representation  of 
that  beam  shown  in  Figure  6.  The  equation  for  the  electric  field  E(a>,r,z )  of  a  Gaussian 
beam  in  the  TEMoo  mode  can  be  expressed  as 
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Figure  6.  Diagram  showing  the  components  of  a  Gaussian  beam  [14]. 

In  this  equation,  w(a>,  z)  is  the  spot  size  of  the  beam  (e1  point  of  the  electric  field),  w 


is  the  beam  waist,  R  (a>,  z)  is  the  curvature  of  the  phase  front,  and  z0  is  the  Rayleigh 
range.  The  spot  size,  radius  of  curvature,  and  Rayleigh  range  are  defined  as 
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A  THz  beam  is  quasi-optical  and  is  strongly  impacted  by  diffraction,  since  the  Rayleigh 
range  at  0.3  THz  is  only  3  cm  and  at  1.0  THz  it  is  10  cm.  Characterizing  THz  radiation 
as  a  Gaussian  beam  is  important  because  it  allows  one  to  numerically  propagate  the  beam 
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through  a  system  of  lenses  given  an  initial  spot  size  and  radius  of  curvature  at  the  exit  of 
the  emitter’s  hemispherical  lens. 


2.3.2.4  Gaussian  Beam  Expansion 


A  Gaussian  beam  can  be  represented  in  terms  of  the  curvature  of  its  phase  front  R 
and  the  spot  size  of  the  beam  w  by  the  complex  beam  parameter  q.  The  equation  for  this 
relationship  is 

1 

q{co,z) 

The  complex  beam  parameter  can  then  be  used  to  propagate  the  Gaussian  beam  through 
an  optical  system  by  the  use  of  the  ABCD  matrix  formalism.  The  output  complex  beam 
parameter  qout  can  represented  in  terms  of  the  input  qin : 
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The  ABCD  matrix  has  the  following  format  for  propagation: 
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The  ABCD  matrix  formulism  allows  one  to  represent  an  entire  optical  system  by 
multiplying  a  series  of  propagation  matrices  together. 
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2.3.3  Scanning  Optical  Delay  Line 


The  THz  TDS  system  requires  a  means  of  varying  the  delay  of  one  optical  beam 
relative  to  a  second  beam  in  order  to  move  the  sampling  gate  across  the  waveform  to  be 
sampled.  This  is  accomplished  by  varying  the  optical  path  length  traversed  by  the  two 
beams  by  using  a  scanning  optical  delay  line  on  one  of  them.  The  optical  delay  is 
achieved  by  mounting  a  retro-reflective  mirror  on  a  mechanical  scanner.  For  this 
application,  a  Newmark  NLS4  mechanical  stage,  remotely  controlled  by  a  computer,  was 
used  to  scan  an  8  cm  diameter  retro-reflective  mirror.  The  delay  line  has  a  maximum 
speed  of  2  inches/second,  a  resolution  of  0.125  pm,  and  a  repeatability  of  5  pm. 

In  many  imaging  applications  it  is  desirable  to  scan  the  optical  delay  as  rapidly  as 
possible  in  order  to  increase  the  waveform  acquisition  rate.  However,  there  are  not  many 
mechanical  devices  that  can  scan  faster  than  a  few  tens  of  Hertz.  Typical  point-by-point 
images  formed  by  imaging  systems  have  a  data  acquisition  rate  of  20  pixels  per  second, 
which  translates  to  a  few  minutes  per  image,  if  one  uses  a  raster  scanner  to  move  the 
sample  through  the  THz  beam  [48], 

2.3.4  Raster  Scanner 

A  Newmark  NLS4-X-25  horizontal  and  an  NLS4-X-16  vertical  stage  were 
connected  to  form  an  X-Y  stage  used  to  perform  a  raster  scan  on  various  samples  for 
imaging.  The  X-Y  stage  was  computer-controlled  and  could  be  synchronized  with  the 
scanning  optical  delay  line  through  the  Galil  motion  control  software.  The  vertical 
component  had  a  maximum  speed  of  0.5  inches/second,  a  resolution  of  0.06  pm,  and  a 
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repeatability  of  5  pm,  while  the  horizontal  component  had  the  same  specifications  as  the 
scanning  optical  delay  line. 

2.3.5  Lock-in  Amplifier 

A  Stanford  Research  Systems  Model  SR850  DSP  lock-in  amplifier  was  used  to 
collect  electronic  data  off  of  the  THz  detector.  The  time  domain  data  of  the  THz  pulse 
could  be  measured  with  either  current  or  voltage.  The  lock-in  amplifier  time  constant 
and  optical  delay  line  speed  had  to  be  carefully  chosen  so  that  the  TDS  system  could 
accurately  collect  the  entire  pulse,  maximize  the  signal-to-noise  ratio,  and  minimize  the 
total  scan  time.  The  sensitivity  and  dynamic  reserve  had  to  be  chosen  to  maximize  the 
dynamic  range  of  the  lock-in  amplifier.  Furthermore,  the  data  collection  rate  and  the 
optical  delay  line  had  to  be  synchronized  to  control  the  bandwidth  and  frequency 
resolution  of  the  THz  TDS  system.  The  bandwidth,  B,  is  determined  by 

B  =  (2A  t)"‘  (20) 

in  which  At  is  the  time  step  between  data  points  collected  in  the  THz  pulse.  The 
frequency  resolution  A/  is  given  by 

Af  =  (2AtN)-l=(2Tyl  (21) 

where  N  is  the  total  number  of  points  and  T  is  the  total  time  of  the  scanned  pulse. 

Finally,  “auto-phase”  was  chosen  when  collecting  the  pulse  data  to  adjust  the  reference 
phase  shift  to  eliminate  any  residual  phase  error  and  eliminate  the  imaginary  component. 
This  maximized  the  amplitude  of  the  peak  of  the  THz  pulse.  Data  collection  could  be 
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accomplished  by  either  manual  entry  of  the  data  through  the  interface  on  the  front  panel 
of  the  lock-in  amplifier  or  through  a  computer  interface  using  LabView. 


2.4  THz  TDS  Data  Analysis  on  Composites 


After  the  data  has  been  collected  from  the  lock-in  amplifier,  various  mathematical 
software  tools  to  analyze  that  data.  Figure  7a  shows  the  time  domain  response  of  a 
typical  THz  TDS  system  transmitting  through  ambient  air  as  the  optical  delay  line  is 
scanned.  The  time  domain  data  can  then  be  transformed  into  the  frequency  domain  via  a 
numerical  fast  Fourier  Transform  (FFT).  A  FFT  uses  an  algorithm  to  solve  a  discrete 
Fourier  Transform  (DFT)  more  quickly.  A  DFT  is  given  by 


N-l  2iri 

E((o)  =  YjE{t)e  N  k  =  0,l,...,N-l 

t= 0 


(22) 


Frequency,  THz 


Figure  7.  (a)  THz  pulse  time  waveform,  (b)  THz  frequency  spectrum  [26]. 


where  E(a>)  is  the  frequency  response  of  the  received  THz  signal,  E(  t)  is  the  time 


domain  response  of  the  received  THz  signal,  and  N  is  the  total  number  of  samples.  A 
discrete  FFT  is  equivalent  to  a  plane  wave  expansion  that  returns  the  complex  amplitude 
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of  the  discrete  frequency  components  con  of  the  THz  pulse.  This  representation  of  the 


THz  electric  field  E  ( a>n )  can  written  as 

E(a>n)  =  E0  [con)ei(k{o,")x+03j)  (23) 

where  E0  is  the  amplitude  and  k  is  the  wave  number  Figure  7b  represents  the  magnitude 

of  the  Fourier  transform  of  the  THz  pulse  signal  in  7a. 

In  THz  TDS,  one  measures  a  sample  Esamp(t j  and  a  reference  pulse  Ere/t)  in  the 
time  domain  and  then  converts  them  both  into  the  frequency  domain,  Esamp  (oj)  and 

Eref  (oj  )  ,  with  an  FFT.  The  pulse  that  is  measured  at  the  receiver  is  a  convolution  of  the 

transfer  functions  of  the  emitter  R(t)  and  the  receiver  S(t).  The  sample  also  imparts  a 
transfer  function  H(t)  to  the  measured  THz  pulse  when  it  is  included  in  the  THz  path.  A 
diagram  showing  the  effect  of  the  sample  on  the  transfer  functions  in  the  time  and 
frequency  domains  is  shown  in  Figure  8.  The  convolution  of  the  sample  transfer  function 
in  the  time  domain  equates  to  a  multiplication  in  the  frequency  domain.  An  example  of 
the  resulting  waveforms  in  the  time  and  frequency  domain  is  illustrated  in  Figure  9. 

It)  0  M  J)  (I  M 

R(a>)S(a))  =  R(t)*S(t)  R(m)H(a>)S(a>)  =  R(t)*H((o)*S(t) 

Figure  8.  Components  of  the  THz  waveform  at  the  detector  in  the  time  and 
frequency  domain  [14]. 
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Figure  9.  Waveforms  resulting  from  the  collection  of  data  at  the  THz  detector  [14]. 

2.4.1  Material  Parameter  Estimation 

Material  parameters,  like  the  index  of  refraction  and  the  absorption  coefficient, 
can  be  measured  using  a  THz  waveform.  The  THz  system  response  of  a  sample  can  be 
determined  by  dividing  the  frequency  response  of  system  with  a  sample  E  (oj)  by  the 

system  response  of  the  reference  Eref  (oj)  .  The  system  response  can  be  shown  with  the 
following  equation 


E„mp  (a>)  =  =  e-, (»(.)-*„), 

Eref(a)  R0  (aj)S(oj)e{oJt-k(,x) 


(24) 


in  which  R0  is  the  amplitude  of  the  THz  emitter  transfer  function,  ko  is  the  wave  number 

in  free  space,  k(a>)  is  the  wave  number  for  the  sample  material,  and  x  is  the  length 
propagated.  The  value  for  k(co)  can  be  written  in  terms  of  its  complex  parts 


(25) 


and  then  added  to  the  version  of  the  sample  THz  system  response  given  in  terms  of 
amplitude  A  and  phase  tj) 
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E sample  j*0)  _  -i[a/c[n\(o)-\)\x  _ 

ErefH 


(26) 


2.4.2  Drude  Model 

By  measuring  the  THz  system  response  in  either  transmission  or  reflection  mode, 
one  can  obtain  the  complex  index  of  refraction  of  the  material.  The  complex  index  of 
refraction  n(co)  can  be  written  as 


n(co)  =  n(oj)  -  i/c(oj) 


(27) 


where  n(a>)  represents  the  real  refractive  index  and  rc(a>)  is  the  extinction  coefficient 


and  is  proportional  to  the  absorption  coefficient  a{co) ,  which  in  cgs  units  is 


k(co)  = 


a(co)c 
2  co 


(28) 


One  of  the  ways  of  predicting  the  optical  properties  of  a  THz  TDS  system  is  to  use  the 
Drude  Model.  The  Drude  Model  offers  one  method  of  calibrating  a  THz  TDS  system  so 
that  the  accurate  performance  of  the  system  can  be  verified. 

Usually  the  optical  constants  of  a  semiconductor  in  the  far-infrared  region  are 
analyzed  with  the  Drude  model.  According  to  the  Drude  model,  the  complex 
conductivity  is  expressed  as 


&(a>)  =  e0  ■ 


co„ 


1 

uo  +  — 

T 


(29) 


and  in  this  expression  r  is  the  momentum  relaxation  time  of  free  carriers  and  cop  is  the 
plasma  angular  frequency  defined  by  cop  =  ^nce2  / £(]m'  .  In  this  plasma  frequency 
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expression  nc  is  the  carrier  density,  e  is  the  charge  of  an  electron,  and  m  is  the  effective 
mass  [67],  The  complex  electrical  conductivity  can  be  written  in  terms  of  its  real  and 
imaginary  parts:  a(co)  =  cr,  (co)  -  icr2  (co)  and 

crj  (co)  =  s0cos2 (co) ,  (30) 

cr2(co)  =  (co)-sJ  (31) 

with  £f  being  the  relative  dielectric  constant  of  the  undoped  semiconductor  and  £0  being 

the  permittivity  of  free  space  [67],  In  a  nonmagnetic  medium,  the  complex  dielectric 
constant  £(co)  is  given  as 

£(oo)  =  n2(co)  =  £-i^^-  (32) 

CO£0 

and  cr(a>)  is  the  complex  conductivity. 

2.4.3  Scattering 

Scattering  is  the  process  by  which  electromagnetic  radiation  is  forced  to  deviate 
from  a  straight  trajectory  by  one  or  more  localized  non-uniformities  in  the  medium 
through  which  it  passes.  There  are  two  types  of  elastic  scattering:  Rayleigh  scattering 
and  Mie  scattering.  Rayleigh  scattering  is  the  process  by  which  radiation  is  scattered  by  a 
small  element,  treated  as  a  small  sphere,  in  which  the  diameter  d  of  the  sphere  is  much 
smaller  than  the  wavelength  A  of  the  radiation.  The  upper  limit  for  scattering  of  this  type 
is  d/X  ~  1/10.  For  elements  with  diameters  comparable  to  the  wavelength,  Mie  scattering 
describes  the  scattering  process,  and  the  shape  of  the  scattering  elements  becomes  more 
relevant.  The  study  of  the  propagation  and  scattering  of  THz  pulses  through  media  has 
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been  studied  recently  [13,  27,  55,  56,  57,  58,  60,  71].  For  nondestructive  evaluation  of 
composites  with  THz  radiation,  one  would  be  interested  in  isolating  air  bubbles,  voids, 
and  delaminations  in  the  material.  One  could  look  at  the  scattering  properties  of  thin 
samples  and  study  the  multiple-scattering,  diffusive,  phenomena  that  can  be  measured 
with  coherent  pulses.  Additionally,  the  scattering  effects  of  the  glass  fibers  and  bundles 
within  the  composite  could  be  studied. 

The  scattering  of  coherent  THz  radiation  Ecoh  (  oj  )  can  be  described  using  the 


scattering  mean  free  path  ls  (co)  as  a  function  of  frequency  is  given  by  the  equation: 


Ecoh  H 

-z 

-a[(o)z 

r  (  \  =GXP 

Ein  (®) 

_24  M_ 

exp 

2 

4  n 


(»  + 1) 


(33) 


where  a  is  the  absorption  coefficient  and  z  is  the  thickness  of  the  sample.  The 
scattering  mean  free  path  is  the  inverse  of  the  scattering  coefficient  ps.  In  order  to  make 
predictions  about  scattering  of  THz  radiation  from  particles  one  would  attempt  to 
determine  the  scattering  cross  section  <js  (  oj  j  of  the  scattering  centers: 


4  H  =  [w  («)] 


(34) 


in  which  n0  is  the  number  density  of  the  scatterers  in  the  material. 

If  the  scattering  centers  are  much  smaller  than  the  THz  wavelength  then  the 
Rayleigh  approximation  can  be  used.  The  scattering  cross  section  for  Rayleigh  scattering 
is  independent  of  the  particle  shape  and  is  given  by 


o’.  =  ■ 


2  x5d6 


3  r 


n2- 1 
y  n  +2  j 


(35) 
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where  n  is  the  index  of  refraction  of  the  scattering  centers  and  d  is  the  diameter  of  the 
spheres.  The  intensity  scattered  per  unit  volume  for  an  unpolarized  electromagnetic  wave 
is  given  by  the  following: 


1=1, 


71 V 


0  8 r2A4 


n_- 1 
y  n  +  2  j 


(l  +  cos2  (#)) 


(36) 


where  6  is  the  angle  measured  from  the  direction  of  propagation  and  r  is  the  distance  to 
the  particle.  The  scattering  efficiency  Qeff  is  commonly  used  in  place  of  the  cross  section 

and  is  related  by 


Qeff  = 


(37) 


in  which  A  is  cross  sectional  area  of  the  scattering  center.  The  equation  for  spheres  is 
more  complicated  when  the  size  of  the  scatterers  are  closer  in  size  to  the  wavelength  of 
the  radiation  [42], 

In  composites,  there  are  strands  of  fiber  woven  together  and  then  held  together 
using  a  resin.  The  fibers  most  closely  resemble  an  infinitely  long  circular  cylinder. 

When  the  diameter  of  the  cylinder  is  sufficiently  small  compared  to  the  wavelength, 
certain  approximations  can  be  used  for  the  scattering  coefficients.  The  scattering  cross 
section  per  unit  length  for  the  TE  and  TM  modes  are 

=^r-(™2-1)2  (38) 


4;rV 


Gte 
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m 


y  111  +  1  j 


(39) 


where  a  is  the  radius  of  the  cylinder  and  m  is  the  ratio  of  the  index  of  refraction  of  the 
scattering  center  to  the  index  of  refraction  of  the  containing  medium.  The  samples  tested 
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in  this  research  had  an  overlapping  weave  pattern  of  glass  strands  with  all  of  the  strands 


oriented  in  one  of  two  directions  at  90°  angles  to  each  other.  The  assumption  was  made 
that  the  overlapping  pattern  would  result  in  an  average  of  the  scattering  cross  section  of 

the  TE  and  TM  components,  a  ~  +  °TE  The  intensity  of  the  TE  and  TM  modes 

at  normal  incidence  are  given  by 

Itm=I (40) 
A  r  v  ' 


J  _  T  4a4 a4  f  m2  - 1  )  n^2(a\ 

1 TE~h  33  2,i  C0S  (^) 

A  r  m  +1 


when  the  imaginary  portion  of  the  refractive  index  is  small  compared  to  the  real  part. 

This  illustrates  that  the  intensity  of  the  scattered  radiation  from  infinitely  long  circular 
cylinders  varies  as  the  fourth  power  of  the  cylindrical  radius  and  inversely  as  the  cube  of 
the  wavelength.  This  limiting  case  for  infinite  cylinders  is  compared  to  the  intensity  of 
Rayleigh  scattering  from  small  spheres,  which  varies  as  the  sixth  power  of  the  radius  and 
inversely  as  the  fourth  power  of  the  wavelength  [42], 

The  radiation  from  a  photoconductive  switch  is  linearly  polarized  according  to  the 
orientation  of  the  emitter  dipole.  For  this  research,  the  emitter  was  oriented  for  vertical 
polarization.  If  the  orientation  of  the  glass  fibers  is  known,  then  one  could  predict  the 
intensity  of  THz  radiation  scattering  from  the  glass  fiber  composite.  The  TM  mode  has 
its  electric  field  parallel  to  the  axis  of  the  cylinder  while  the  TE  mode  has  its  axis 
perpendicular  to  the  axis  of  the  cylinder.  Using  this  orientation  where  the  axis  of  the 
cylinder  is  vertical,  one  can  determine  the  intensity  of  THz  radiation  that  is  scattered 
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given  an  incident  linear  polarization  angle.  The  intensity  for  the  vertically  and 
horizontally  polarized  components  can  be  defined  as 

lv  =  ~1)2  s™2  (x)  (42) 

A  r 

.44/  2  1 V 

cl  m  1  9  /  „\  9  /  \  N 

4=/ 0— ^ - 2 — 7  cos  l6*)008  U)  (43) 

/t  r  ym  +1 J 

in  which  %  is  defined  as  the  angle  from  the  positive  horizontal  axis  and  r  is  the  radial 
distance  from  the  x-axis.  A  diagram  of  the  configuration  is  shown  in  Figure  10  with  the 
scattered  radiation  propagating  as  a  cylindrical  wave  [42],  The  glass  weave  strands  could 
be  oriented  at  any  angle  to  the  THz  polarization  direction,  at  normal  incidence. 


x 


Figure  10.  Diagram  showing  incident  and  scattered  ray  upon  an  infinitely  long  cylinder. 

An  alternate  approach  has  been  presented  which  models  scattering  effects  within 
a  densely  distributed  array  of  infinite  cylinders.  Its  model  utilizes  an  effective  field  and  a 
quasi-crystalline  approximation  for  the  multiple  scattering  interactions  within  the  dense 
medium  [45]. 
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2.4.4  Effective  Medium  Approximations 

The  propagation  of  electromagnetic  waves  in  heterogeneous  media  can  be 
analyzed  using  various  effective  medium  approaches  [1,  13,  27,  29,  49,  74],  Among 
these  approaches,  there  are  variations  of  three  that  have  been  used  most  frequently, 
including  the  Maxwell-Gamett  theory,  the  Bruggeman  effective  medium  model,  and  the 
simple  effective  medium  approximation.  All  of  these  theories  assume  that  spherical 
inclusions  have  been  dispersed  into  a  host  matrix,  the  inclusions  are  much  smaller  than 
the  wavelength  in  the  host  material  (<  0.17,  -  0.27.),  and  that  the  material  can  therefore  be 
treated  like  an  effective  medium  as  shown  in  Figure  11  [6].  The  model  can  also  be 
extended  to  cover  cylindrical  inclusions  by  using  dielectric  tensors. 


Inhomogeneous  Material  Effective  Medium 


Figure  11.  An  inhomogeneous  system  is  replaced  by  a  (virtual)  effective  medium. 

The  simplest  effective  medium  approach  is  with  linear  interpolation,  using  the 

filling  factor  f  which  defines  the  volume  fraction  of  the  particles: 

£eff  (®)  =  A  O)  +  (1  -  f)£h  O)  (44) 

in  which  st  is  the  dielectric  constant  of  the  inclusions  and  sh  is  the  dielectric  constant  of 
the  host  material.  Another  model  that  is  used  is  called  the  Maxwell-Gamett  effective 
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medium  theory: 


Seff{0))  =  £h(C0) 


1  +  2/[(g,-g/i)/(g/+2g/,)]] 

l-f[(e,-eh)/(et+2sh)]  J 


(45) 


This  theory  assumes  that  the  total  polarization  should  be  the  same  as  the  sum  of 
polarization  contributions  from  the  embedded  particles.  The  Maxwell-Gamett 
approximation  holds  up  well  at  low  filling  factors  / until  the  filling  factor  reaches  a  level 
of  approximately  0. 1 5.  At  this  point,  the  theory  developed  by  Bruggeman  must  be  used 
to  achieve  an  accurate  determination  of  the  effective  dielectric  constant.  The  formula  is 
given  by 
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£i  Seff 

\£i+K£eff  J 


+  (1  -/) 


g/i  ge# 

\£k  +  ^£eff  J 


=  0 


(46) 


where  K  is  a  geometric  factor.  The  value  for  K  is  1  if  the  inclusions  are  an  array  of 
cylinders  with  their  axes  collinear  with  the  incident  radiation  and  the  value  for  K  is  2  if 
they  are  spherical  nanoparticles  [8], 


2.4.5  Imaging 

The  THz  wave  distribution  function  and  the  target’s  refractive  index  can  be 
described  by  Maxwell’s  Equation: 

V2E(co,r,t)  +  ^fi  j^E(co,r,t)  =  0  (47) 

where  r  denotes  a  point  in  space,  c  is  the  speed  of  light,  and  n  is  the  complex  refractive 
index  of  the  medium.  In  the  case  where  the  polarization  effects  of  the  wave  can  be 
neglected  in  the  temporal  frequency  domain,  this  equation  can  be  written  as  the  scalar 
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Helmholtz  equation 


W2  E(co,r)  +  klh(co,r)2  E(co,r)  =  0  (48) 

in  which  E(co,r)  is  the  complex  amplitude  function  of  the  electromagnetic  field  and  ko 
represents  the  wave  number.  This  Helmholtz  equation  is  valid  in  imaging  and 
tomography  applications.  The  purpose  of  THz  imaging  is  to  determine  the  refractive 
index  function  n(co,r )  using  the  measured  THz  wave  amplitude  and  phase  distribution 
function  E(co,r)  [75], 

THz  imaging  can  be  performed  in  either  transmission  or  reflection  mode.  A  THz 
TDS  imaging  system  is  constructed  with  a  focus  at  the  desired  location  of  the  object  to  be 
imaged.  The  waveform  of  the  THz  pulse  that  goes  through  the  object  at  that  focus  can 
then  be  measured.  By  translating  the  object  and  measuring  the  transmitted  waveform  for 
each  position  of  the  object,  an  image  can  be  built  pixel  by  pixel.  The  two  dimensional 
image  that  is  formed  may  be  represented  by  such  aspects  of  the  waveform  as  amplitude, 
phase,  or  a  combination  of  the  two.  Phase,  or  time  delay,  imaging  consists  of  the 
extraction  of  the  arrival  time  of  the  waveform  at  each  pixel.  For  amplitude  imaging,  one 
could  use  the  peak-to-peak  amplitude  of  the  time-domain  waveform  for  a  measure  of 
bandwidth-averaged  transmission.  Alternatively,  one  could  perform  an  FFT  of  each 
pixel’s  waveform  and  then  integrate  over  any  desired  frequency  band  to  display  some 
part  of  the  spectral  information  [12,  48], 

Three-dimensional,  THz  tomography  refers  to  the  cross-sectional  imaging  of  an 
object  by  measuring  either  the  transmitted  or  reflected  illumination.  One  way  of 
accomplishing  this  is  to  measure  the  time-domain  THz  pulse  transmitted  through  or 
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reflected  by  the  material.  The  resulting  waveform  will  consist  of  a  series  of  reflections 
from  the  sample  of  the  various  dielectric  interfaces  within  that  sample.  The  polarity  and 
magnitude  of  each  reflection  are  given  by  the  reflection  coefficient  at  each  interface  and 
are  related  to  the  size  and  sign  of  the  corresponding  index  step  [48].  This  form  of  THz 
tomography  is  restricted,  though,  in  its  applicability.  This  type  of  tomography  is  based 
on  three  assumptions:  (1)  the  samples  have  no  dispersion,  absorption,  or  diffraction 
affects,  (2)  the  reflections  are  so  weak  that  multi-reflections  can  be  neglected,  and  (3)  the 
refractive  index  is  uniform  throughout  each  layer  [75].  Any  one  of  the  affects  listed  can 
significantly  alter  the  shape  of  the  reflected  waveform.  A  numerical  Fourier 
deconvolution  in  the  time-domain  could  be  applied  to  the  waveform  to  more  accurately 
measure  the  position  of  each  of  the  reflections.  Other  methods  of  tomography  are 
diffraction  tomography,  computed  tomography,  binary  lens  tomography,  optical 
coherence  tomography,  and  digital  holography  [75]. 

2.4.6  Image  Processing 

The  image  of  a  composite  sample  can  be  formed  by  arranging  the  pixels  that  were 
collected  during  the  THz  TDS  scan.  Since  THz  images  are  formed  from  a  limited 
number  of  samples,  their  images  appear  “pixelated.”  Interpolation  techniques  can  be 
employed  which  enable  one  to  remove  the  sharp  edges  of  the  pixels  to  make  the  image 
easier  to  view  visually.  Interpolation  is  the  estimate  of  a  signal  at  locations  other  than 
where  it  was  sampled.  If  a  two  dimensional  array’s  pixels  are  expanded  equivalently  in 
both  dimensions,  the  interpolation  is  called  zoom  or  electronic  zoom.  Various 
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interpolation  techniques  in  the  time  domain  include  the  ideal  sine,  zero  order,  triangular, 
or  cubic  spline  interpolator  [32], 

A  frequency  domain  interpolation  technique  was  applied  for  the  THz  images  in 
this  research.  The  sampling  rate  was  increased  by  zero-padding  in  the  spatial  frequency 
domain.  First,  the  two-dimensional  image  can  be  transformed  to  the  frequency  domain 
with  a  discrete  cosine  transform  (DCT).  A  DCT  is  similar  to  the  discrete  Fourier 
Transform  (DFT)  except  that  it  uses  only  the  real  numbers,  and  is  frequently  used  in  data 
compression  algorithms  because  it  tends  to  concentrate  information.  Then  the  resulting 
matrix  can  be  zero-padded  in  the  spatial  frequency  domain  which  increases  the  sampling 
rate  in  the  spatial  domain.  Finally,  an  inverse  DCT  can  be  performed  on  the  data  which 
results  in  the  display  of  much  smaller  pixels  and  an  apparent  smoothing  of  the  original 
pixel  edges.  By  using  the  interpolation  technique  outlined  above,  zero-padding  out  to  ten 
times  the  length  and  width  in  the  frequency  domain,  the  new  image  will  appear  to  look 
the  same  as  if  the  original  image  was  shrunk  ten  times. 

MATLAB  also  has  a  function  called  ‘interp’  which  is  designed  to  smooth 
boundaries  between  pixels.  However,  when  this  function  is  used  on  its  own,  it  tends  to 
create  a  blurry  image  if  the  pixels  are  too  large.  Occasionally,  the  DCT  can  leave  a 
periodic  pattern  in  the  image  from  sharp  boundaries  which  can  be  further  smoothed  out 
with  a  filter.  A  filter  that  has  been  used  successfully  for  THz  imaging  by  Zhang  et  al.  is 
the  Wiener  filter  [75], 
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2.5  THz  TDS  Literature  Review 


This  section  summarizes  some  of  the  recent  work  that  has  been  published  on  the 
NDE  of  materials  using  THz  TDS  and  on  the  progress  in  modeling  the  THz  electrical  and 
optical  properties  of  those  materials.  This  is  followed  by  a  summary  of  some  of  the  work 
done  on  imaging  and  on  characterizing  the  THz  properties  of  composites,  ceramics,  and 
semiconductors,  as  well  as  paint  and  foam  adhering  to  metallic  surfaces. 

2.5.1  Composites 

A  composite  is  defined  in  the  American  Heritage  Dictionary  as:  “A  complex 
material,  such  as  wood  or  fiberglass,  in  which  two  or  more  distinct,  complementary 
substances,  especially  metals,  ceramics,  glasses,  and  polymers,  combine  to  produce 
structural  or  functional  properties  not  present  in  any  individual  component”  [3]. 
Composite  designers  can  choose  from  a  variety  of  fiber  reinforcements  and  resins  to 
develop  the  laminate  structure.  The  reinforcement  material  provides  mechanical 
properties  such  as  stiffness,  tension  and  impact  strength,  while  the  resin  system  (matrix) 
provides  physical  properties  including  resistance  to  fire,  weather,  ultraviolet  light  and 
corrosive  chemicals.  Three  factors  must  be  considered  when  choosing  reinforcements. 
One  must  first  choose  the  type  of  fiber:  usually  fiberglass,  aramid  (in  Kevlar),  or  carbon 
fiber  (also  called  graphite).  Then  one  must  choose  the  form:  either  roving  strands,  mat, 
or  fabrics.  Finally,  an  orientation  of  fiber  direction  is  chosen,  including  parallel, 
circumferential,  or  helical  along  the  length  of  the  part,  and/or  with  random  continuous 
strands.  The  three  types  of  resins  that  are  available  are  polyester,  vinyl  ester,  and  epoxy. 
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A  form,  either  wet  lay-up  or  prepreg  (a  reinforcement  saturated  with  resin),  must  be 
carefully  chosen  to  ensure  a  successful  design  [78]. 

Several  articles  have  been  written  describing  the  response  of  THz  radiation  when 
transmitted  through  various  composites.  One  of  the  types  of  composites  that  has  been 
studied  involves  polymers,  which  can  be  mixed  with  other  additives  or  fdlers  to  yield 
compounds  with  modified  physical  properties.  In  the  manufacture  of  these  composites,  it 
can  be  difficult  to  verify  the  degree  of  homogeneity  of  the  resulting  compound,  especially 
when  using  nanoscale  fillers.  One  research  group  studied  samples  of  low-density 
polyethylene  (LDPE)  containing  titanium  dioxide  nanospheres  (diameter  ~  270  nm) 
coated  with  silver  nanoparticles  (~  20  nm)  as  well  as  glass-fiber  reinforced  with  epoxy 
[64,  66].  For  the  first  case,  they  were  able  to  see  inhomogeneous  areas  in  the  LDPE 
polymer,  but  were  not  able  to  unambiguously  determine  whether  the  inhomogeneities 
were  a  result  of  the  internal  content  or  the  surface  reflections  from  the  rough  surface.  In 
the  case  of  the  glass  fiber  matrix,  THz  TDS  was  able  to  locate  materials  and  defects 
hidden  inside  the  sample.  A  separate  article  characterizes  the  refractive  index  and 
absorption  coefficient  of  polymers,  including  high-density  polyethylene,  Picarin,  and 
cyclic  olefin  polymer  [69].  All  of  these  polymers  can  be  used  to  fabricate  THz 
communication  systems.  The  authors  showed  that  up  to  90%  of  the  incident  THz 
radiation  passes  through  2  mm  of  each  of  the  materials  with  very  little  attenuation, 
making  them  good  candidates. 

Another  composite  that  has  been  studied  is  the  polymer-graphite  composite  [70], 
These  types  of  composites  are  excellent  absorbers  of  THz  radiation  with  good  shielding 
efficiency  and  small  reflection.  THz  TDS  was  used  to  measure  the  complex  refractive 
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index  of  this  composite,  showing  that  with  a  35.7%  graphite  concentration,  reflection  was 
less  than  10%,  while  the  absorption  was  over  80%.  Finally,  a  separate  study  has  been 
conducted  on  aerospace  carbon- fiber  composites  and  space  shuttle  insulating  tiles, 
exploring  the  effectiveness  of  using  THz  radiation  for  NDE  of  heat  damage  and  voids 
[62],  The  article  was  focused  primarily  on  comparing  continuous  wave  (CW)  THz 
techniques  to  pulsed  techniques,  showing  that  CW  THz  systems  have  sufficient 
resolution  to  detect  voids  and  heat  damage. 

2.5.2  Ceramics 

Although  ceramics  can  be  considered  a  subset  of  composites,  a  ceramic  is  more 
specifically  defined  in  Wikipedia  as  “an  inorganic  non-metallic  material  whose  formation 
is  due  to  the  act  of  heat”  [78],  Ceramics  are  used  in  a  wide  range  of  applications 
including  jet  engines  and  mirrors.  Silicon  carbide  (SiC)  is  a  ceramic  material  which  is 
strong  and  light-weight  and  is  an  excellent  choice  to  be  used  for  space-based  mirrors.  It 
is  a  difficult  material  to  mount,  and  the  unrelieved  stress  that  results  from  mounting  can 
cause  failure  during  launch.  The  photoelastic  effect,  measured  in  the  THz  spectrum,  can 
be  used  to  measure  the  stress  induced  in  the  SiC.  Recent  literature  has  shown  that  the 
photoelastic  effect  has  been  used  effectively  to  characterize  defects  in  silicon  carbide 
wafers  in  the  optical  regime  [43,  46], 

Studies  have  recently  been  conducted  using  THz  TDS  on  aluminum  oxide 
bearings  and  glass  ceramics  to  measure  their  optical  constants  and  to  probe  them  for 
existing  fractures  [62,  63,  68],  Measurements  have  also  been  made  on  the  dielectric 
constants  and  loss  tangents  of  low-loss,  high  dielectric  constant  ceramics  at  THz 
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frequencies  [11].  The  materials  tested  included  alumina,  steatite,  titania  loaded 
polystyrene,  and  zirconium-tin-titanate.  The  real  part  of  the  dielectric  constants  of  these 
materials  were  6.5,  9.3,  16,  and  37-93  respectively  for  each  of  them  and  could  be  used  for 
cost  effective  photonic  crystal  components.  Finally,  material  parameters  were  measured 
on  various  building  materials,  like  plaster,  glass,  and  wood  [59]. 

2.5.3  Paint  Thickness  over  Metallic  Surfaces 

The  external  coating  of  nearly  all  military  aircraft  is  stripped  down  to  bare  metal 
during  programmed  depot  maintenance  cycles.  This  paint  stripping  process  has  become 
cost  prohibitive  in  recent  years  and  continues  to  be  a  major  problem  for  the  sustainment 
of  aging  aircraft  fleets.  One  of  the  main  contributing  factors  is  the  requirement  to 
visually  inspect  aircraft  for  corrosion  and  fatigue  damage.  Several  nondestructive 
techniques  have  been  used  to  avoid  visual  inspection,  including  infrared,  eddy  current, 
ultrasonic,  radiographic,  and  near- field  microwave.  THz  offers  the  advantage  of  being 
non-contact,  penetrating  nonconductive  materials,  and  reflecting  from  both  metallic 
materials  and  from  impedance  discontinuities  like  voids  and  inclusions  [16]. 

One  paper  presented  an  approach  to  use  THz  tomography  to  measure  the 
thickness  and  degree  of  drying  in  the  painting  process.  The  thickness  resolution  dmm  is 
the  limiting  factor  when  using  THz  radiation  and  is  given  by 

dmia=(cAT)/(2ng )  (49) 

where  c  is  the  speed  of  light  in  a  vacuum,  AT  is  the  temporal  width  of  the  THz  pulse, 
and  ng  is  the  group  index  of  refraction  for  the  medium  [87],  This  paper  reported  that  the 
resolution  was  limited  to  ~  10  pm.  Another  paper  reported  that  surface  roughness 
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differences  below  ~  10  pm  were  indistinguishable  from  the  metallic  background  when 
searching  for  corrosion  under  paint  [5]. 


2.5.4  Semiconductors 


THz  TDS  can  also  be  used  to  measure  the  optical  and  electrical  properties  of  a 
semiconductor  by  the  technique  outlined  in  section  2.4. 1 .  Such  techniques  have  been 
applied  to  various  semiconductor  materials  and  their  nanostructures,  including  quantum 
dots  [54],  photonic  crystals,  and  carbon  nanotubes,  which  can  have  both  metallic  and 
semiconductor  properties.  Photonic  crystals,  made  from  semiconductor  materials,  have 
been  studied  recently  using  THz  TDS  to  determine  both  defects  and  group  velocity  and 
group  velocity  dispersion  characteristics  [37,  38,  39,  53],  The  dielectric  constant  of 
carbon  nanotubes  has  been  measured  with  a  combination  of  the  Drude  and  Lorentz 
harmonic  oscillator  models  [41,  47].  In  the  dielectric  function  s(a>)  shown  in  the 
following  equation,  the  second  term  in  the  sum  represents  the  Drude  model  and  the  third 
term  represents  the  Lorentz  model: 


e{a>)  = 


CD„ 


CD 


CD  +  iT  CD 


-  +  - 


p,L 


(< orL  -  or  )  -  iVLoj 


(50) 


The  term  srfj  is  the  dielectric  constant  at  infinity,  cop  is  the  plasma  frequency,  T  is  the 


damping  rate  in  the  Drude  model,  col  is  the  phonon  frequency,  T L  is  the  spectral  width, 

and  cOp,L  is  the  strength  of  the  Lorentz  oscillator.  Further  characterization  using  the  Drude 
model  could  also  allow  one  to  determine  the  number  of  photocarriers  which  have  been 
generated. 
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Carrier  dynamics  can  also  be  determined  by  using  a  technique  called  Optical- 
pump  THz-probe  spectroscopy  to  measure  the  conductivity  of  photocarriers  as  a  function 
of  time  after  photoexcitation  [2,  10].  In  this  photo  luminescence  setup,  an  additional 
probe  beam  from  the  laser  source  is  used  to  optically  illuminate  the  sample  and  generate 
carriers  so  that  it  can  then  be  probed  with  a  THz  beam.  The  optical-pump  THz-probe 
technique  has  been  used  extensively  to  study  semiconductor  carrier  dynamics.  It  has  also 
been  used  to  measure  the  optical  and  electronic  characteristics  of  the  material  by  using 
the  free  electron  Drude  model  or  the  Lorenz  model  to  estimate  the  dielectric  constant. 

2.5.5  Measuring  Material  Properties  with  THz  Radiation 

Significant  efforts  have  been  expended  on  the  development  of  THz  imaging 
systems;  however,  somewhat  less  attention  has  been  paid  to  the  issue  of  the  propagation 
through,  or  scattering  from,  realistic  target  materials  which  might  be  investigated  with 
such  systems.  Understanding  the  interaction  of  THz  radiation  with  a  wide  range  of  both 
homogeneous  and  inhomogeneous  materials  is  an  essential  part  of  the  design  of  systems 
for  THz  sensing  and  imaging.  This  section  attempts  to  provide  a  summary  of  some  of  the 
recent  work  done  on  modeling  the  scattering  of  THz  radiation  and  on  the  estimation  of 
material  parameters  using  THz  TDS. 

2.5.5.1  Transmission  Mode  Measurement  Technique 

One  method  of  determining  the  material  parameters  of  a  homogeneous,  planar 
material  in  using  THz  TDS  in  transmission  mode  is  outlined  by  Domey  et  al  [23].  A 
schematic  of  the  interaction  of  a  THz  pulse  with  a  sample  is  shown  in  Figure  12.  An 
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incident  pulse  contacts  the  sample  where  part  of  the  wave  is  reflected  and  part  is 
transmitted.  As  the  THz  wave  propagates  through  the  material  it  interacts  with  that 
material.  Multiple  reflections  can  also  occur  at  the  internal  interfaces  of  the  material, 
creating  multiple  output  pulses. 


Figure  12.  Transmission  and  reflection  of  a  THz  wave  through  a  planar,  homogeneous  material. 
The  Fresnel  equations  describe  the  transmission  and  reflection  of  the  THz  wave  at 

each  interface 


lJ0J)  = 


2  na  ( co )  cos  6 


ha  ( co )  cos  /?  +  hb  ( co )  cos  6 


(51) 


raM)  = 


nb  (< co )  cos  6  -  na  (co)  cos  (3 
na  (co)  cos  [3  +  hb  (co)  cos  6 


(52) 


where  tab(co)  is  the  transmission  coefficient  of  a  wave  at  an  incident  angle  0  from  region 
a  to  region  b  and  rab(co)  is  the  reflection  in  region  a  at  the  a-b  interface.  The  angle  6  is 
estimated  and  the  angle  /?  is  an  approximation  of  Snell’s  law  for  nearly  transparent 


materials: 
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I 3  =  arcsin 


n  sin  6 


(53) 


As  the  wave  moves  through  the  material  along  the  path  designated  d,  it  accumulates  a 
phase  change  through  given  by  pd(a>,d )  according  to  a  transfer  function  or  THz  system 
response 

,  ( -in,(oj)ojd )  .... 

pd(co,d)  =  ex p  - - -  (54) 

\  c  J 

where  c  is  the  speed  of  light  and  the  ray  length  d  is  determined  by 

d  =  ~ — •  (55) 

cos  p 

Scattering  is  neglected  in  this  model  including  such  effects  as  the  interface  roughness 
[23], 

The  THz  path  both  with  and  without  the  sample  must  be  considered.  For  the  free- 
air  path,  or  reference  path,  we  have 

Eref(a>)  ^ initial  (56) 

with  the  distance  x  being  the  distance  between  the  transmitter  and  receiver  and  na;r  being 
the  index  of  refraction  of  air  (1.00027).  For  the  propagation  through  a  sample,  the 
distance  m  through  the  sample  must  be  subtracted  from  the  distance  traveled  through  air. 
This  distance  varies  with  the  angle  of  incidence  and  the  material  parameters.  This 
distance  m  is  defined  as: 

m  =  d  cos(#  -  /?) .  (57) 

In  Figure  12  there  are  only  two  additional  paths  shown  representing  the  Fabry-Perot 
effect,  but  several  additional  paths  could  also  be  shown.  The  equation  for  the  primary 
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received  signal  through  the  sample  is  given  by 

Eprimaryi®)  =  E initial  i™)  P , ■*[«,(*  "  ™)]  X 'oi# !ample(CO,d)tl0  .  (58) 

In  this  expression,  t0l  is  the  transmission  coefficient  at  the  air-sample  interface  and  tl0  is 

the  transmission  coefficient  at  the  sample-air  interface.  Now  that  we  have  both  the 
reference  wave  and  the  sample  wave  defined,  we  can  define  the  ratio  of  the  two 
E  (w)  4n  .  (co)n  ,  (a>)  co  Sco  /§ 

primary  x  /  air  ^  ^  sample  v  7  7 

£„:/ M  [njoj) co  p  +  nsample(w) co  <S]2 

This  provides  the  transfer  function  for  the  model,  where  only  the  complex  function  for 
nsamPie(C0)  >  U  and  9  are  free  variables  [23], 

This  procedure,  outlined  by  Domey  et  al.,  involves  a  three  step  iterative  process 
of  determining  the  thickness  and  the  complex  index  of  refraction.  The  first  step  is  to 
make  an  initial  guess  at  the  thickness  of  the  sample.  Second,  the  algorithm  calculates  the 
beginning  functions  for  the  complex  index  of  refraction,  assuming  a  non-dispersive 
medium.  Finally,  a  gradient  descent  algorithm  iterates  the  complex  index  of  refraction  in 
frequency  until  the  total  error  no  longer  decreases  monotonically.  The  complex  index  of 
refraction  is  then  computed  in  this  manner  for  a  range  of  thicknesses.  A  metric  called 
total  variation  is  then  used  to  select  which  pair  of  refractive  index  and  thickness  should 
be  selected  [23]. 

2.5.5.2  Literature  Summary  of  THz  Material  Parameter  Measurements 

The  technique  for  modeling  the  material  parameters  of  a  homogeneous,  planar, 
optically  thin,  low-absorption  material  is  shown  in  section  2.5.5. 1.  This  technique  was 
designed  by  Duvillaret  et  al.  and  refined  by  Domey  and  Mittleman  [22  -  25].  Although 
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Duvillaret  et  al.  claim  that  an  accurate  initial  guess  at  the  sample  thickness  is  very 
important  to  determine  the  material  parameters,  Domey  et  al.  claim  that  they  are  able  to 
accurately  determine  sample  thickness  and  refractive  index  by  initially  bounding  the 
upper  and  lower  limits  of  the  thickness  / 


in  which  At  is  the  time  delay  between  the  reference  pulse  and  the  first  pulse  from  the 
sample.  The  parameters  «,  and  n2  are  the  limits  of  the  range  of  refractive  index,  of 


which  the  range  1.2-8  covered  most  practical  cases  [23],  Withamayachumnankul  et  al. 
were  able  to  use  a  fixed-point  iterative  method  to  determine  the  complex  refractive  index 
of  the  sample  within  a  few  iterations,  given  the  sample  thickness.  Using  the  same 
measurement  techniques  as  Domey  and  Duvillaret,  their  iterative  process  gave  the  correct 
parameters  for  a  silicon  wafer  sample  when  the  material  thickness  exceeded  200  pm  [79, 
80], 

Xu  et  al.  prepared  a  simulation  with  a  Lorentz  oscillator  model  to  evaluate  the 
accuracy  of  the  extracted  optical  constants  in  THz  TDS  [85].  They  analyzed  the 
influence  of  the  angle  deflection,  the  delay-line  timing  error,  and  background  and  1/f 
noise  on  the  THz  TDS  measurement  precision.  They  found  that  the  deflection  of  the 
incident  angle  can  influence  both  the  real  and  imaginary  parts  of  the  index  of  refraction 
and  the  spectral  features  of  the  errors  show  a  parabolic  shape.  They  also  found  that  the 
delay-line  timing  uncertainty  introduced  an  exponential-dependent  deviation  curve  of  the 
imaginary  part  in  the  frequency  domain  and  a  dramatic  feature  in  the  resonant  region  for 
the  real  part. 
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A  separate  study  reported  on  the  propagation  of  a  single-cycle  THz  pulse  in  a 
random  media.  Pearce  et  al.  demonstrated  a  unique  method  for  characterizing  the 
diffusive  radiation  emerging  from  a  random  medium,  involving  the  measurement  of  both 
its  amplitude  and  phase  [56].  They  were  able  to  compute  a  time- windowed  temporal 
correlation  between  fields  measured  at  different  spatial  locations  which  enabled  them  to 
identify  individual  scattering  events  for  imaging  applications.  Pupeza  et  al.  have 
developed  a  refined  approach  to  measuring  optical  material  parameters  [61]  while  Kuzel 
et.  al.  has  extracted  both  the  complex  refractive  index  and  the  wave  impedance  using  THz 
TDS  [44], 

2.5.6  THz  TDS  for  NDE 

This  subsection  outlines  some  of  the  recent  literature  reported  specifically  on  the 
use  of  THz  for  NDE  and  quality  control.  In  composite  research,  fiber  glass  samples  from 
the  automotive  industry  have  been  analyzed  for  their  material  property  characteristics  and 
have  been  imaged  to  show  damage  from  stress  [64,  66],  Some  of  the  polymers  that  are 
typically  used  as  the  matrix  material  in  composites  have  also  been  researched  [76,  77], 
The  first  article  looks  at  the  use  of  THz  imaging  to  inspect  plastic  weld  joints.  They  were 
able  to  use  THz  TDS  to  reveal  contaminates  like  metal  or  sand  within  a  weld  joint  of  two 
polyethylene  sheets  and  to  identify  small  air  gaps  where  the  welding  process  had  failed. 

In  the  second  article,  the  authors  were  able  to  determine  the  additive  content  of  polymers 
used  in  composites,  including  magnesium  hydroxide  in  polyethylene,  glass  fibers  in 
polyamide,  and  calcium  carbonate,  silicon  or  silicon  dioxide  in  polypropylene. 
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THz  NDE  is  also  being  proposed  as  a  method  of  inspecting  the  adhesive  bond¬ 
line  of  space  shuttle  tiles  for  defects  [4,  84,  92],  Space  shuttle  tiles  are  manufactured 
from  silica  into  a  very  porous  material.  These  tiles  are  then  bonded  to  the  bottom  of  the 
space  shuttle  to  provide  insulation  for  re-entry.  A  triangular  shaped  defect  was  identified 
using  data  from  a  THz  NDE  scan  and  post-processing  with  a  Hilbert-Huang  Transform 
[4],  This  process  successfully  highlighted  the  boundary  between  the  triangular  shaped 
defect  and  the  background. 

There  are  other  articles  which  have  been  written  to  describe  the  use  of  THz  TDS 
for  NDE  of  materials  other  than  composites.  One  of  the  authors  presents  work  using  THz 
TDS  in  reflection  mode  to  inspect  mural  paintings  for  damage.  They  are  able  to  image  a 
sub-surface  painting  of  a  butterfly  hidden  beneath  various  colors  of  paint  [34],  A 
separate  application  uses  THz  TDS  to  detect  foreign  bodies  in  chocolate  bars,  including 
metallic  objects,  stone,  glass,  and  plastics.  It  can  also  distinguish  these  objects  from  nuts 
and  other  food  products  within  the  bars  [40]. 
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3.  Material  Parameter  Estimation  with  Terahertz  TDS 

In  the  nondestructive  testing  of  aircraft  composites,  one  is  ultimately  interested  in 
identifying  damage,  including  damage  from  heat  and  mechanical  stress,  and  the  detection 
of  voids  and  delaminations.  Measuring  the  material  properties  of  a  composite  provided 
insight  into  the  fundamental  interaction  of  THz  radiation  with  the  material  and  could  be 
tested  as  a  technique  to  see  if  the  material  properties  of  a  damaged  composite  were 
changed.  A  major  focus  of  this  chapter  is  to  present  the  results  obtained  when  THz  TDS, 
in  transmission  and  reflection  configuration,  was  used  to  measure  the  index  of  refraction 
and  absorption  coefficient  of  various  materials.  In  addition,  scattering  and  effective 
medium  theory  were  studied  for  their  applicability  to  composites.  Finally,  this  material 
parameter  measurement  technique  was  applied  to  composite  samples  with  bum  damage. 

3.1  THz  TDS  System  Characterization 

There  were  several  challenges  associated  with  setting  up  a  THz  TDS  system, 
including  the  accurate  alignment  and  measurement  of  laser  pump  and  probe  beams,  THz 
beam  alignment,  and  the  maximization  of  THz  output  signal  quality.  A  diagram  of  our 
first  setup  for  material  parameter  measurements  is  shown  in  Figure  13,  using  only  two 
parabolic  mirrors  since  the  beam  did  not  need  to  be  focused  onto  the  samples.  To 
determine  the  point  at  which  the  TDS  system  measurements  were  of  sufficient  quality, 
we  started  with  an  initial  set  of  criteria  for  comparison.  The  next  sections  will  detail 
these  criteria  and  the  methodology  for  analyzing  the  THz  output,  including  the 
measurement  of  water  absorption  spectra.  This  will  be  followed  by  a  subsection  on  the 
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modeling  of  the  THz  pulse  and  associated  spectrum.  This  will  be  followed  by  an 
explanation  of  some  of  the  difficulties  that  we  had  initially  with  the  PCA  detector  and  of 
the  techniques  that  we  used  to  help  align  the  THz  TDS  system  more  quickly. 


Figure  13.  THz  TDS  setup  for  material  parameter  measurements. 


3.1.1  THz  TDS  Output  Specifications 

The  output  specifications  that  we  used  to  help  evaluate  our  THz  TDS  system  were 
detailed  by  Cheville  et.  al.  in  their  THz  tutorial  [14].  These  criteria  are  compared  in 
Table  1  along  with  the  measured  results  from  our  THz  TDS  setup.  Although  one  could 
measure  the  signal-to-noise  ratio  from  a  logarithmic  plot  of  the  THz  amplitude  spectrum, 
Cheville  et.  al.  recommend  measuring  the  signal-to-noise  ratio  and  the  dynamic  range  in 
the  time  domain.  The  signal-to-noise  level  was  measured  by  fixing  the  optical  delay  line 
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at  the  position  of  peak  pulse  amplitude  and  then  collecting  data  for  30  minutes  at  that 
position.  One  could  then  measure  the  mean  of  the  peak  amplitude  and  observe  how  much 
the  signal  drifted.  The  signal-to-noise  ratio  is  given  by 


Signal 

Noise 


v 

pk-mean 

NsP 

v  2  J 


(61) 


when  the  THz  signal  is  measured  in  volts.  The  dynamic  range  is  then  measured  by 
collecting  a  scan  of  the  pulse,  making  sure  to  take  sufficient  time  domain  data  in  front  of 
the  pulse  (t  <  0).  The  data  is  then  scanned  to  find  the  smallest  voltage  step  in  this  period 
prior  to  the  pulse  and  is  then  compared  to  the  peak  of  the  pulse 

V 

dynamic  range  = - — — .  (62) 

^  ^min(/<0) 

The  measurement  accuracy  specification  was  made  by  comparing  10  separate  THz  TDS 
scans  and  comparing  the  absolute  variations  in  the  amplitude  spectra  between  them. 


Table  1.  THz  Time  Domain  Spectrometer  specifications  and  measurements. 


Specification 

Typical  Value 

Conditions 

AFIT  Typical 
Measured  Value 

Amplitude 
Dynamic  Range 

20,000:1 

Measured  in  time  domain 

14,000:1 

Power  Dynamic 
Range 

400,000,000:1 

Amplitude  value  squared 

196,000,000:1 

Median 

Frequency 

0.654  THz 

Determined  by  spectral 
amplitude 

0.704  THz 

Bandwidth 

0.06  to  2.51  THz 
0.12  to  1.64  THz 
0.26  to  0.92  THz 

2%  points  of  peak  amplitude 
10%  points  of  peak  amplitude 
50%  points  of  peak  amplitude 

0.01  to  2.5  THz 
0.05  to  1.9  THz 
0.15  to  1.35  THz 

Signal-to-Noise 
(Signal  Drift) 

±  0.5%  from  mean 

Measured  sitting  at  peak  of 
pulse  over  30  minute  period 

±  2.75%  from  mean 
standard  deviation 
=  0.9% 

Measurement 

Accuracy 

±  0.5%  @  0.2  THz 
±  0.5%  @  0.5  THz 
±  0.25%  @  1  THz 
±  2%  @  2  THz 

No  signal  averaging, 

10  minute  scan  length 

±3.23%  @  0.2  THz 
±0.22%  @  0.5  THz 
±  5.14%  @  1  THz 
±  20%  @  2  THz 
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3.1.2  Water  Vapor  Absorption  Spectra 


Another  method  that  was  used  to  determine  if  our  THz  TDS  setup  was  producing 
accurate  results  was  to  look  at  water  vapor  absorption  spectra.  The  water  vapor  in  the 
atmosphere  leaves  fluctuations  after  the  main  pulse  in  the  time  domain  and  it  leaves 


characteristic  notches  in  the  frequency  spectrum.  A  full  scan  of  the  THz  pulse  through 
air  was  taken  by  moving  the  optical  delay  line  at  a  slow  rate.  The  results  of  the  scan  are 
shown  in  Figure  14.  It  should  be  noted  that  the  smaller  pulse  in  14a  is  a  result  of 
reflections  in  the  ZnTe  crystal  which  cause  the  Fabry-Perot  etalon  effect  in  14c.  The 


phase  plot  in  14b  shows  how  the  phase  unwraps  linearly  with  frequency. 


Freq  (THz) 


(c) 

Figure  14.  THz  TDS  scan  of  ambient  air:  (a)  THz  time  domain  pulse,  (b)  THz 
phase  spectrum,  and  (c)  THz  amplitude  spectrum. 
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The  absorption  lines  in  14b  can  now  be  identified  in  the  THz  frequency  spectrum  and  can 
be  compared  with  those  measured  and  calculated  in  the  literature.  As  shown  in  Table  2, 
our  absorption  lines  line  up  fairly  closely  with  published  values  between  1.0-  2.4  THz. 
Several  articles  have  been  written  measuring  the  absorption  spectrum  of  water  at  various 
relative  humidity  readings  and  studying  the  various  molecular  transitions  and  phenomena 
associated  with  absorption  [83,  90]. 

Table  2.  Absorption  spectrum  rotational  transitions  for  water  vapor  between  0.2  -  2.4  THz  [82]. 


Transition 

Number 

Published 

Frequency 

(THz) 

Measured 

Frequency 

(THz) 

Av  (THz) 
Published  - 
Measured 

1 

0.557 

0.56 

-0.003 

2 

0.753 

0.75 

0.003 

3 

0.989 

0.99 

-0.001 

4 

1.098 

1.10 

-0.002 

5 

1.164 

1.16 

0.004 

6 

1.2264 

1.23 

-0.0036 

7 

1.4116 

1.41 

0.0016 

8 

1.603 

1.60 

0.003 

9 

1.671 

1.67 

0.001 

10 

1.718 

1.71 

0.008 

11 

1.798 

1.80 

-0.002 

12 

1.869 

1.87 

-0.001 

13 

1.921 

1.92 

0.001 

14 

2.076 

2.07 

0.006 

15 

2.264 

2.27 

-0.006 

16 

2.346 

2.35 

-0.004 

17 

2.393 

2.40 

-0.007 

3.1.3  THz  Pulse  and  Spectrum  Modeling 

The  shape  of  the  THz  pulse  was  very  sensitive  to  the  position  and  alignment  of 
the  laser  pump  and  probe  beams  as  well  as  the  THz  optics.  One  could  attempt  to  predict 
the  shape  of  the  THz  beam  in  both  space  and  time  by  approximating  it  as  a  Gaussian 
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pulse.  In  one  article,  the  authors  have  already  attempted  to  model  the  effects  of  EO 
sampling  on  a  Gaussian  pulse  emitted  from  a  GaAs  crystal  and  detected  by  a  ZnTe 
crystal  [9].  Beside  the  effect  of  the  THz  emitter  and  detector,  both  the  alignment  of  the 
lenses  and  mirrors  as  well  as  the  sample  have  an  effect  on  the  eventual  shape  of  the  pulse. 
In  this  exercise,  the  emphasis  was  on  modeling  the  effects  of  a  glass  fiber  composite 
sample  on  the  THz  pulse  and  comparing  modeled  data  to  measured  data. 

When  pulses  reach  the  single  cycle  regime,  like  they  do  with  THz  TDS,  one  must 
use  the  following  equation  to  model  the  electric  field  E(r,t)  of  the  pulse: 


E(r,t)  =  exp 


-w0^irG0r 


mrG02T](t) 


where  w0  is  the  spot  size  of  the  THz  beam,  cox  is  the  central  frequency  of  the  THz 
spectrum,  rG0  is  the  pulse  duration,  r  is  the  distance  from  the  center  of  the  beam,  L  is  the 
distance  the  pulse  has  traveled,  and  v  is  the  speed  of  the  pulse  through  the  medium 
modified  by  the  complex  index  of  refraction  [19].  The  terms  n(t)  and  ?G  are  defined  as: 


L  r2 

rl(t)  =  t  — —  and 

v  2  Lv 


r  2+  ^ 

rco  lc 


In  modeling  the  THz  pulse  shape  in  the  time  domain,  the  spot  size  was  assumed  to 


have  a  width  %  of  0.6  mm  which  became  mostly  irrelevant  when  considering  the  center 


of  the  THz  spot  at  r  =  0.  Both  the  THz  pulsewidth  and  frequency  center  of  mass  had  to 
be  measured  from  actual  data.  Data  was  chosen  carefully  so  that  the  majority  of  the  pulse 
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had  a  positive  polarity  like  a  Gaussian.  The  full  width  at  half  maximum  (FWHM)  value 
of  tgo  was  measured  to  be  0.62  psec  and  the  THz  central  frequency  was  0.55  THz. 

After  choosing  the  THz  signal  characteristics,  the  Gaussian  pulse  through  the  air 
reference  was  modeled  in  the  time  domain  and  frequency  domain.  The  glass  fiber 
complex  index  of  refraction  was  applied  to  the  frequency  domain  data  using  2.05  -iK(co) 

where  K(co)  is  the  imaginary  part  of  the  index  of  refraction.  Then  an  inverse  Fourier 
transform  was  applied  to  return  the  pulse  to  the  time  domain.  The  comparisons  of  the 
modeled  and  measured  pulses  and  frequency  spectra  can  be  seen  in  Figure  15.  Although 
a  Gaussian  distribution  for  the  pulse  was  assumed,  the  FWHM  value  rG  -  t pj\.\ll  for  a 
Gaussian  was  not  small  enough  to  match  the  measured  frequency  response.  A  factor  of 
tg  =  r;;/l  .763  ,  similar  to  a  sech  distribution,  was  more  accurate  in  modeling  the  THz 

pulse  and  spectrum  characteristics.  The  magnitude  of  the  modeled  air  reference  pulse 
was  chosen  to  match  the  peak-to-peak  amplitude  of  the  measured  THz  pulse.  It  is 
apparent  that  factors  other  than  the  sample  have  an  effect  on  the  THz  signal,  such  as  the 
THz  emitter,  detector,  optics,  and  water  vapor  absorption. 


.lodded  Gaussian  Pulse  Through  Air 
Modeled  Gaussian  Pulse  'llvough  Glass  Fiber 
ircd  THz  Pulse  Tlnough  Air 
Measured  THz  Pulse  Through  Glass  Fi 


2  3 

Time  (sec) 


x  10 


Figure  15.  Comparison  of  a  Gaussian  modeled  THz  pulse  through  air  and  glass 
fiber  compared  to  THz  TDS  measurements  for  air  and  glass  fiber  in  the  (a)  time 
domain  and  (b)  frequency  domain. 
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3.1.4  System  Setup  Challenges 


The  most  challenging  part  of  collecting  data  from  a  THz  TDS  system  is  the 
alignment  of  the  optical  components  with  the  near-infrared  laser  and  the  terahertz  beam. 
The  laser  has  to  be  aligned  with  the  photoconductive  switch  antenna  and  with  the  detector 
(either  a  photoconductive  switch  or  an  electro-optic  crystal).  Then  the  THz  beam  had  to 
be  aligned  through  the  optical  system  of  parabolic  mirrors  to  maximize  the  THz  pulse 
received  at  the  detector.  One  piece  of  equipment  that  was  essential  to  the  THz  beam 
alignment  was  a  pyroelectric  detector.  This  pyroelectric  detector  could  be  used  to 
maximize  the  beam  intensity  at  each  of  the  optical  components  before  locking  them  into 
place  and  helped  to  manage  the  many  degrees  of  freedom  that  existed  in  the  alignment 
problem.  A  HeNe  laser  could  be  used  to  initially  help  align  the  near-infrared  Ti:Sapphire 
laser  with  the  emitter  and  detector.  Then  fine-tuning  of  the  position  of  the  laser  spot 
could  be  steered  by  measuring  the  resistance  of  the  emitter  and  detector.  The  laser  power 
sent  to  the  pump  and  probe  beams  also  had  to  be  closely  monitored  so  that  the  maximum 
specifications  were  not  exceeded.  A  list  of  the  specifications  that  were  measured  are 
listed  in  Table  3. 


Table  3.  List  of  Specifications  for  Ekspla  photoconductive  switch  components. 


Measurement 

Measured  Value 

Maximum  Limit/ 
Expected  Value 

Emitter  Dark  Resistance 

Infinite 

Infinite 

Emitter  Illuminated  Resistance 

3.1  UQ 

8MQ 

Detector  Dark  Resistance 

68  UQ 

67  MQ 

Detector  Illuminated  Resistance 

1.3  UQ 

0.8  MQ 

Pump  Laser  Power 

125  mW 

125  mW 

Probe  Laser  Power 

50  mW 

50  mW 

Power  Supply  Voltage 

50  Volts 

50  Volts 

Power  Supply  Modulation 

120  kHz 

50  kHz 
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Significant  effort  was  expended  to  align  the  Ekspla  photoconductive  detector, 
before  it  was  considered  damaged  or  unusable.  The  pulse  widths  that  were  measured 
were  much  wider  in  the  time  domain  and  much  lower  in  amplitude  than  the  advertised 
values  in  the  Ekspla  manual.  Figure  16  shows  the  (a)  THz  pulse  and  its  (b)  amplitude 
spectrum  measured  at  the  photoconductive  switch  receiver.  In  the  time  domain,  the  main 
pulse  was  spread  out  over  6  psec  and  several  oscillations  proceeded  for  a  long  duration 


after  the  main  pulse. 


Figure  16.  THz  pulse  measured  with  photoconductive  (PC)  switch  detector:  (a)  time 
domain  pulse  and  (b)  frequency  domain  spectrum  compared  to  electro-optic 
detector. 


It  was  also  two  orders  of  magnitude  lower  in  amplitude  than  the  pulse  measured  with  the 
ZnTe  crystal.  This  translated  to  a  narrower  bandwidth  (~0.7  THz)  and  a  reduced 
amplitude  in  the  frequency  domain.  The  energy  in  the  pulse  was  determined  by 
estimating  the  noise  floor  for  each  spectrum  and  then  calculating  the  area  under  the  curve 
above  the  noise  floor  for  each  spectrum  separately.  The  total  energy  under  the  PC  switch 
spectrum  was  only  2%  of  that  under  the  EO  detector  spectrum. 

3.2  Transmission  Mode  Measurements 

This  section  outlines  the  results  of  our  transmission  mode  measurements  of 
optical  parameters  using  THz  TDS.  The  first  two  subsections  show  the  calculation 
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techniques  used  to  determine  the  material  parameters  and  the  results  of  those 
calculations.  The  last  subsection  shows  how  to  use  the  Drude  model  to  calibrate  the 
material  parameter  measurements  from  the  THz  TDS  system  using  a  semiconductor. 


3.2.1  Calculation  Technique 


Although  the  theory  for  calculating  materials  parameters  was  discussed  in  section 
2.4.1,  the  practical  technique  used  in  this  research  for  determining  indices  of  refraction 
and  absorption  coefficients  is  outlined  in  this  section.  The  index  of  refraction  «(<u)  was 
calculated  using  the  following  equation 

n  (a>  =  1  +  — - -  (66 

y  ’  cod 

sam 

where  (f>  is  the  phase  of  the  sample  (sam)  or  the  reference  (ref),  c  is  the  speed  of  light, 
and  dsam  is  the  sample  thickness.  One  of  the  challenges  of  using  this  formula  was  to 
accurately  unwrap  the  phase  data  linearly  between  the  sample  and  reference  since 
software  routines  alternate  between  -ji  and  +n.  The  absorption  coefficient  a(oj)  could 
be  calculated  using  the  following  formula 


«(<») 


-In 


■M 


T(co)\Eref{co) 


(67) 


in  which  | Esam  (o)|  is  the  magnitude  of  the  THz  field  collected  through  the  sample,  and 

Eref  (<y)|  is  the  magnitude  of  the  THz  field  collected  through  air.  T  (&>)  is  the  fraction  of 
the  power  transmitted  through  the  air-sample  interface  given  by 
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(68) 


,  .  4 n(a>) 

T(a>)  = - ^-5-. 

(n(®)  +  l) 

The  index  of  refraction  and  the  absorption  coefficient  could  be  determined  with  an 
iterative  process,  solving  both  of  them  simultaneously  at  each  frequency  component. 

Since  the  absorption  coefficient  calculations  were  somewhat  noisy,  a  windowing  filter 
could  be  applied  which  averaged  the  ten  nearest  data  points  at  each  frequency. 

When  making  material  parameter  measurements,  one  must  also  pay  attention  to 
the  dynamic  range  DR  of  the  THz  TDS  setup.  The  largest  absorption  coefficient  that  can 
be  measured  reliably  with  a  given  DR  corresponds  to  the  situation  in  which  the  sample 
signal  is  attenuated  to  a  level  equal  to  the  noise  floor: 

amd  =  2 In  DR-^  .  (69) 

L  ("  +  0  J 

If  the  absorption  coefficient  exceeds  «max  then  the  analysis  will  show  an  apparent 
absorption  coefficient  equal  to  amax .  The  index  of  refraction  is  not  subjected  to  the  same 
limitations  [35], 

Experimental  uncertainty  or  error  is  another  consideration  in  making  calculations 
with  measured  data.  One  must  first  estimate  the  measurement  error  for  all  of  the 
potential  sources,  given  as  Sx,  Sy,  etc.  If  the  source  of  uncertainty  is  from  a  standard 

deviation,  it  can  be  adjusted  by  the  factor  1  /yfn  if  n  is  the  total  number  of  runs.  Once 
this  is  accomplished  then  the  total  error  Sf  of  a  function / can  be  expressed  by 
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There  are  several  potential  sources  of  error  with  a  summary  of  a  few  of  them  listed  in 
Table  4. 


Table  4.  Sources  of  error  for  THz  material  parameter  measurements. 


Source  of 

Measurement  Error 

Error  Type 

Error 

Material  Parameter 
Measurement  Error 

Amplitude  Accuracy  of 
Frequency  Spectrum 

Amplitude 
(std  dev  -  10 
measurements) 

±  3.29%  @  0.2  THz 
±  5.7%  @  0.5  THz 
±  4.46%  @  1  THz 
±  40.7%  @  2  THz 

Absorption 

coefficient 

Phase  Accuracy  of 
Frequency  Spectrum 

Phase 

(std  dev  -  10 
measurements) 

±  1.01%  @0.2  THz 
±  2.5%  @  0.5  THz 
±  0.15%  @  1  THz 
±  14.2%  @  2  THz 

Index  of  refraction 

Sample  thickness 

Length 

±0.01  mm 
(w/calipers),  ±1  mm 
(w/ruler) 

Index  of  refraction 
and 

Absorption 

coefficient 

Sample  angle 

Angle  from 
normal 

±5  degrees 

Absorption 
coefficient 
±0.4%  amplitude 

Optical  delay  line  timing 

Length 

repeatability 

±5  pm 

Index  of  refraction, 
Absorption 
coefficient 

±16.66  femtoseconds 

3.2.2  Measurement  Results  on  Various  Materials 

Using  the  THz  TDS  setup  in  Figure  13,  data  was  collected  from  numerous 
samples,  including  space  shuttle  tile,  silicon  wafers,  aircraft  glass  fiber  composite  strips, 
an  aircraft  panel  made  with  composites,  plexiglass,  and  paper.  The  measurement 
technique  from  the  previous  section  was  then  applied  and  the  results  are  summarized  in 
the  following  paragraphs.  The  first  material  that  was  measured  was  space  shuttle  tile,  a 
sample  that  was  borrowed  from  AFRL/RXLP.  The  sample  was  5  cm  in  thickness  and 
was  composed  of  a  mixture  of  70-80%  air  and  the  remainder  of  silica.  The  sample  was 
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scanned  and  the  resulting  pulses  and  amplitude  spectra  are  shown  in  Figure  17.  One  can 
observe  that  the  amplitude  of  the  pulse  is  reduced  and  the  pulse  broadens  after  traveling 


Figure  17.  THz  TDS  results  through  air  &  space  shuttle  tile:  (a)  THz  pulse  and  (b) 
amplitude  spectra. 


through  the  space  shuttle  tile.  The  bandwidth  of  the  pulse  is  reduced  from  2.5  THz  down 
to  0.7  THz.  The  next  step  was  to  calculate  the  index  of  refraction  and  the  absorption 
coefficient.  The  results  are  shown  in  Figure  18.  The  index  of  refraction  remained 
constant  at  1.16  over  the  range  of  0.1  -  0.7  THz,  while  the  absorption  coefficient  rose  to 
0.9  cm"1  at  0.5  THz  and  peaked  at  about  1.4  cm'1  at  0.7  THz.  However,  a  calculation  of 
amax  in  Figure  18b  shows  that  the  absorption  measurements  are  only  accurate  until  0.72 
THz.  At  this  frequency,  the  dynamic  range  of  the  sample  spectrum  is  insufficient  to 


Figure  18.  Space  shuttle  tile  calculations  for  (a)  index  of  refraction  and  (b) 
absorption  coefficient  with  a  calculation  for  a,,,,,. 
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measure  the  absorption  coefficient  and  in  the  rest  of  the  graph  noise  dominates  the 
sample  spectrum.  Measurements  in  the  THz  frequency  range  on  silica  have  been 
reported  in  the  literature,  ranging  between  1.93  -  1.975  for  refractive  index  and  between 
0.15  -  0.62  cm'1  for  absorption  coefficient  at  0.5  THz  [51]. 

Since  space  shuttle  tile  is  a  mixture  of  silica  and  air,  one  would  have  to  use  an 
effective  medium  approximation  to  estimate  the  index  of  refraction  for  the  space  shuttle 
tile.  Assuming  that  the  true  index  of  refraction  for  air  is  1  and  silica  is  1.96,  an  estimate 
was  made  for  the  filling  factor / for  silica  given  our  calculated  index  of  1 . 16.  The  results 
are  summarized  in  Table  5.  Since  the  Maxwell-Gamett  Model  is  generally  considered 
valid  when  f  <  0.15,  the  Bruggeman  Model  appears  to  be  the  most  accurate  choice  since 
it  yields  a  filling  factor  estimate  within  the  bounds  of  70  -  80%  air  composition.  Silica 
aerogel,  a  similar  thermal  insulating  material  consisting  of  silica  and  air,  was  measured 
with  an  index  of  refraction  of  ~1.03  and  an  absorption  that  ranged  from  0.5  cm'1  at  0.5 
THz  to  1.3  cm'1  at  1.0  THz  [91]. 


Table  5.  Summary  of  results  for  effective  medium  approximations  on  space  shuttle  tile  at  0.5  THz. 


Effective  Medium 
Model 

Mathematical  Formula 

Silica  Filling 
Factor  (/)  Estimate 

Linear  Model 

£eff  (ffl)  =  f  S{  (< CD) )  +  (1  -  f)Eh  (< CD) ) 

12.5% 

Maxwell-Gamett 

Model 

(Valid:  f<0.15) 

,  ,  ,  Jl  +  2 /[(*, -**)/(*, +2**)]1 

«,<»)-».«»)  j  1_/[(4i_,t)/(4i+24t)]j 

21.5% 

Bruggeman  Model 

/  \  /  \ 

£ .  £  rr  £  i  £  rr 

/  *  +0-/)  *  =o 

V£i+K£eff)  {£k+K£effJ 

24% 

In  the  earlier  calculation  of  the  absorption  coefficient  we  assumed  that  all 
reduction  in  amplitude  of  the  signal  was  due  to  absorption.  It  is  also  realistic  to  assess 
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that  scattering  plays  an  important  role  in  a  sample  with  a  mixture  of  two  substances  with 
different  dielectric  constants.  Using  an  average  value  of  0.36  cm'1  from  the  reported 
values  for  the  absorption  coefficient  of  silica  and  the  equation  in  section  2.4.3,  one  could 
estimate  the  scattering  mean  free  path  to  be  approximately  1 .2  cm  at  1  THz. 

A  slice  of  a  glass  fiber  composite  from  was  also  scanned  with  our  THz  TDS 
system  so  that  material  parameters  could  be  estimated.  The  THz  pulse  was  measured  on 
the  glass  fiber  sample  (1.2  mm  thickness)  and  then  an  FFT  was  performed  to  convert  it 
into  the  frequency  domain  as  illustrated  in  Figure  19a  and  b.  The  graphs  showing  the 
calculated  index  of  refraction  and  absorption  coefficient  are  shown  in  Figure  19c  and  d. 
The  measured  value  of  refractive  index  is  constant  across  the  frequency  spectrum  at  2.05, 
while  the  absorption  coefficient  increases  to  7  cm'1  at  0.5  THz  and  25  cm'1  at  1.0  THz.  A 
plot  is  also  shown  in  19d  which  shows  amax  for  the  glass  fiber.  The  highest  frequency  at 
which  there  is  sufficient  dynamic  range  in  the  sample  spectrum  is  at  about  1.25  THz. 

This  calculation  neglects  any  affects  of  scattering.  Other  attempts  to  measure  fiber  glass 
have  yielded  a  refractive  index  of  2.05  -  2. 1  and  an  absorption  coefficient  of  7  cm'1  at  0.5 
THz  and  25  cm'1  at  0.8  THz  [66], 

Another  material  that  was  measured  was  an  entire  aircraft  panel.  The  panel  was 
approximately  2.8  cm  in  thickness  and  two  thin  layers  of  glass  fiber  surrounding  a  thicker 
layer  of  honeycomb  material  in  the  middle.  Only  a  very  limited  amount  of  THz  energy  at 
the  lower  end  of  the  THz  spectrum  penetrated  the  material  and  no  further  analysis  was 
attempted.  Plexiglass  was  also  scanned,  and  the  resulting  index  of  refraction  was 
calculated  as  1.55,  compared  with  a  published  optical  value  of  1.48.  A  paper  notebook 
was  also  measured,  resulting  in  an  index  of  1 .45  compared  with  wood,  which  was  been 
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(C)  (d) 


Figure  19.  THz  pulse  measured  through  glass  liber:  (a)  time  domain,  (b)  frequency 
domain,  (c)  index  of  refraction,  and  (d)  absorption  coefficient  with  amax  included. 


measured  in  the  lower  end  of  the  THz  region  as  1.35.  When  sheets  were  removed  from 
the  notebook,  the  spectral  energy  under  the  curve  changed  in  proportion  to  the  number  of 
sheets  that  were  removed. 

3.2.3  Effective  Medium  Approximation  and  Scattering  Analysis 

Using  a  combination  of  effective  medium  theory  and  scattering  theory  in 
conjunction  with  THz  TDS  measurements,  the  propagation  of  THz  through  a  glass  fiber 
composite  can  be  explained.  First,  the  material  properties  of  weaved  glass  fiber, 
polyimide  resin,  and  the  glass  fiber  composite  were  measured  separately  using  THz  TDS. 
Next,  effective  medium  approximations  for  the  composite  were  used  to  show  the 
correlation  between  the  composite  and  its  components.  Finally,  scattering  theory 
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calculations  were  made  to  show  the  minimal  contribution  of  scattering  to  the  extinction 
coefficient. 

In  our  system,  the  glass  fiber  weave  was  considered  as  the  inclusion  and  the 
polyimide  was  the  host  material.  The  glass  fiber  weave  is  made  up  of  individual 
cylinders  of  glass  that  are  approximately  10  pm  in  diameter  and  much  smaller  than  the 
THz  wavelength.  These  cylinders  are  then  bundled  together  into  groups  that  were  about 
520  pm  in  width  which  are  then  used  to  weave  the  overlapping  pattern  used  in  the 
composite.  The  polyimide  resin  is  added  and  permeates  throughout  the  individual  glass 
fiber  strands.  The  glass  fiber  volume  concentration  is  typically  between  40  -  60% 
depending  on  the  type  of  technique  used  for  combining  the  two  components. 

In  an  effort  to  determine  the  scattering  effects  of  THz  radiation  in  a  glass  fiber 
composite,  samples  of  polyimide  and  glass  fiber  weave  were  prepared.  The  material 
parameters  were  determined  first  using  THz  TDS  and  thickness  measurements.  The 
results  for  polyimide  (0.85  mm  thick)  in  Figure  20  show  a  refractive  index  that  remained 
fairly  constant  around  1.81  ±  0.03  and  absorption  coefficients  of  3  ±  2  cm'1  (0.5  THz),  8 
±  2  cm'1  (0.8  THz),  and  11  ±  2  cm'1  (1.0  THz). 


Fn\|  (THz) 


Figure  20.  THz  TDS  material  parameter  measurements  using  actual  thickness 
measurements  from  polyimide  showing  (a)  index  of  refraction  and  (b)  absorption 
coefficient. 
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It  was  difficult  to  obtain  a  thickness  measurement  on  the  weaved  glass  and 
therefore  the  refractive  index  and  absorption  coefficient  in  the  THz  region  could  not  be 
accurately  measured.  Other  measurements  have  been  made  on  a  variety  of  glasses  in  the 
THz  frequency  range  [51],  but  no  information  is  yet  available  on  the  glass  used  in  this 
research.  One  of  the  most  common  types  of  glass  used  in  aerospace  glass  fiber 
components  is  S-2  glass.  A  literature  source  reported  an  index  of  refraction  of  2.21-2.30 
at  1  MHz  [18].  However,  since  the  glass  fiber  composite  could  be  measured,  effective 
medium  theory  was  used  to  estimate  the  approximate  dielectric  properties  of  the  glass. 
Then  the  glass  and  polyimide  measurements  were  used  to  estimate  the  effects  of 
scattering  in  the  glass  fiber  composite.  Using  the  Bruggeman  model  with  a  glass  fiber 
refractive  index  of  2.03,  a  polyimide  index  of  1.80,  and  assuming  the  composite  is  made 
up  of  60%  glass  by  volume  (common  for  high  performance  composites),  the  refractive 
index  for  the  glass  was  estimated  to  be  2.22.  Another  paper  extracted  the  index  of 
refraction  for  glass  in  a  glass  fiber  composite  in  the  THz  frequency  range  to  be  2.17  [65]. 

It  would  also  be  useful  to  determine  the  amount  that  scattering  contributes  to  the 
extinction  coefficient  versus  absorption.  Since  the  diameters  of  the  individual  strands  of 
glass  (10  pm)  are  much  smaller  than  the  THz  wavelengths  used  in  this  laboratory  setup, 
the  scattering  analysis  of  infinitely  long  circular  cylinders  from  section  2.4.3  is 
comparable  to  Rayleigh  scattering  for  small  spheres.  Using  the  values  of  the  refractive 
index,  the  value  of  m,  which  appears  in  section  2.4.3,  was  calculated  to  be  1.23.  An 
estimate  of  the  scattering  coefficient  for  the  glass  fiber  composite  was  then  determined, 
assuming  that  the  two  perpendicular  components  of  the  scattering  coefficient  contributed 

equally,  G  ~  ^G™  +  Gte  )/2  .  This  assumption  was  made  based  on  the  overlapping  weave 
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pattern  of  the  glass.  The  results  of  the  scattering  approximation  are  shown  in  Figure  21. 
Since  the  estimate  was  low  in  comparison  to  the  extinction  coefficient  presented  in 
Figure  19,  scattering  was  considered  to  be  insignificant  in  the  glass  fiber. 


Figure  21.  Scattering  coefficient  calculation  for  glass  fiber  in  the  THz  frequency  range. 

A  further  attempt  was  made  to  use  the  Bruggeman  effective  medium 
approximation  to  show  the  relationship  between  the  absorption  properties  of  the 
polyimide  and  the  glass  weave  within  the  composite  and  to  explain  the  scattering 
properties  of  the  glass  weave  within  the  composite.  Although  the  glass  weave  (in  air) 
thickness  could  not  be  measured,  two  assumptions  were  made  about  its  composition.  The 
first  was  that  the  glass  had  an  index  of  refraction  of  2.22  and  the  second  was  that  the 
weave  was  60%  glass  and  40%  air.  By  inserting  these  assumptions  into  the  Bruggeman 
model,  the  glass  and  air  mixture  was  calculated  to  have  a  refractive  index  of  1 .73.  One 
could  then  analyze  the  data  and  modify  the  thickness  until  this  index  was  obtained.  The 
extinction  coefficient  could  then  be  determined  at  this  thickness  with  the  results  shown  in 
Figure  22a.  One  can  see  that  there  are  two  valleys  located  at  about  0.45  THz  and  0.9 
THz.  This  corresponds  to  THz  wavelengths  of  approximately  330  pm  and  660  pm  in  air, 
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which  is  the  same  order  of  magnitude  of  the  width  of  fiber  bundles.  This  could  indicate 
that  Mie  scattering  from  the  fiber  bundles  is  present  in  the  glass  weave  and  air  mixture. 

Using  the  absorption  coefficient  data  collected  for  polyimide,  one  could  then  use 
the  Bruggeman  model  to  estimate  the  absorption/extinction  coefficient  of  the  composite. 
A  comparison  of  the  data  from  the  Bruggeman  approximation  and  the  actual  glass  fiber 
composite  data  (without  coating)  is  displayed  in  Figure  22b.  The  two  are  a  close  match 
and  Bruggeman  approximation  shows  the  valleys  from  the  glass  weave/air  mixture  data. 
Then  a  calculation  of  the  scattering  coefficient  was  made  for  the  glass  weave  and  air 
mixture  in  Figure  22c,  assuming  the  glass  fiber  is  much  smaller  than  1  and  ignoring 


Figure  22.  (a)  Extinction  coefficient  estimation  on  glass  weave  assuming  an  air-glass 
index  mixture,  (b)  Comparison  of  the  Bruggeman  EMA  for  glass  fiber  extinction 
using  polyimide  and  glass  weave  data  vs.  actual  glass  fiber  extinction  data,  (c) 
Scattering  coefficient  calculations  for  the  air  and  glass  weave  mixture,  (d) 
Comparison  of  EMA  for  glass  fiber  using  polyimide  and  glass  weave  (scattering 
removed)  vs.  glass  fiber  extinction  coefficient. 
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possible  Mie  scattering  from  the  fiber  bundles.  Finally,  the  scattering  from  the  glass-air 
mixture  was  removed  in  Figure  22d  and  the  Bruggeman  approximation  is  again  compared 
to  the  actual  composite  data.  The  contribution  of  scattering  to  the  extinction  coefficient 
would  be  more  significant  for  the  glass  weave/  air  mixture  than  the  composite  due  to  the 
larger  index  gradient  between  glass  and  air.  However,  since  the  Bruggeman 
approximation  for  the  extinction  coefficient  in  22b  and  22d  produced  such  close  results 
with  the  actual  glass  fiber  extinction  data,  it  could  not  be  determined  whether  the 
scattering  from  the  individual  cylinders  of  glass  fiber  contributes  to  the  overall  extinction 
coefficient. 

3.2.4  Material  Parameter  Calibration 

For  calibration  of  a  material  parameter  measurement  technique,  one  can  use  a 
material  with  absorption  coefficient  and  refractive  index  properties  that  are  well  known. 

A  material  with  a  low  absorption  coefficient  like  teflon  could  be  used  to  make 
measurements  on  the  refractive  index  only.  However,  a  semiconductor  material  like 
silicon  could  be  used  to  calibrate  both  refractive  index  and  absorption  coefficient  through 
the  use  of  the  Drude  Model.  One  can  use  silicon  in  its  high-resistivity  form  to  calibrate 
the  index  of  refraction  and  then  use  a  more  highly  doped  silicon  sample  to  calibrate  the 
absorption  coefficient. 

3.2.4.1  Drude  Model  Predictions 

Silicon  wafers  have  several  specifications  that  are  measured  prior  to  sale, 
including  thickness,  resistivity,  dopant,  orientation,  and  whether  it  was  polished.  We 
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were  looking  for  tight  tolerances  on  resistivity  and  thickness  so  that  we  could  model  the 
optical  parameters  accurately.  The  first  step  was  to  use  the  Drude  Model,  explained  in 
section  2.4.2,  to  predict  the  index  of  refraction  and  absorption  coefficient  for  silicon.  A 
series  of  samples  with  a  wide  range  of  resistivities  were  modeled  and  the  results  are 
shown  in  Table  6.  The  value  used  for  s,r ,  the  relative  dielectric  constant  of  undoped 

silicon,  was  chosen  to  be  1 1.68,  even  though  published  values  ranged  from  1 1.6  -  1 1.9. 
The  carrier  density  of  the  dopant  could  be  estimated  from  graphs  in  the  Sze  text 
Semiconductor  Devices:  Physics  and  Technology  or  from  the  website  for  Solecon 
Laboratories,  while  values  for  mobility  of  carriers  could  be  estimated  from  graphs  in  the 
Sze  text  [73],  The  fraction  of  THz  that  would  be  transmitted  is  given  by 

f -a  I 'T  4  n  1 

Fraction  =  exp  -  - — -  (71) 

^ 2  Ak+i)  J 

which  is  important  when  choosing  a  wafer,  since  a  thicker  wafer  with  higher  absorption 
rates  will  result  in  no  THz  radiation  being  detected.  The  thickness  of  high-resistivity 


Table  6.  Drude  Model  parameters  and  predictions  for  Silicon  at  1  THz. 
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silicon  required  to  remove  the  Fabry-Perot  oscillations  is  determined  by  the  optical  path 
length  through  silicon.  A  thickness  of  at  least  6  mm  is  required  to  remove  the  oscillations 
from  a  133  psec  scan,  assuming  nSi  =  3.4176. 


3.2.4.2  Silicon  Measurements 


The  first  silicon  that  we  measured  was  a  silicon  wafer  with  an  unknown  resistivity 
and  a  measured  thickness  of  approximately  0.37  pm.  The  pulse  was  measured  first  and  is 
shown  in  Figure  23a.  One  can  see  the  first  and  second  Fabry-Perot  reflections  from  the 
sample.  The  Fourier  transform  of  the  pulse  is  shown  in  23b.  The  Fabry-Perot  reflections 
in  23a  led  to  the  oscillations  in  the  refractive  index  and  absorption  coefficient  spectra. 


(c)  (d) 

Figure  23.  Material  parameter  measurements  on  a  silicon  wafer:  (a)  THz  pulse 
showing  FP  reflections,  (b)  THz  spectrum,  (c)  index  of  refraction,  and  (d)  absorption 
coefficient. 
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The  second  piece  of  silicon  that  was  tested  was  a  piece  of  polished  high  resistivity 
silicon  with  a  thickness  of  1.016±0.02  mm  and  a  resistivity  of  30-100  Q-cm.  From  this 
resistivity,  the  Drude  Model  prediction  would  be  3.4176  for  the  refractive  index  and 
0.0008  -  0.0025  cm’1  for  the  absorption  coefficient.  The  pulse  and  spectrum  data  are 
shown  in  Figure  24a  and  24b  with  the  strong  Fabry-Perot  etalon  effect  in  the  frequency 
domain  being  a  result  of  the  reflections  visible  in  the  time  domain.  The  index  of 
refraction  in  24c  oscillates  between  3.375  and  3.475  with  an  average  of  3.425.  The 
absorption  coefficient  24d  oscillates  around  0  cm'1. 


Freq  (THz) 


Frcq  (THz) 


Figure  24.  Material  parameter  measurements  on  high  resistivity  silicon:  (a)  THz 
pulse,  (b)  THz  spectrum,  (c)  index  of  refraction,  and  (d)  absorption  coefficient. 

The  third  piece  of  silicon  was  a  thick  piece  of  unpolished  float  zone  (FZ)  silicon 

with  a  thickness  of  6.0±0. 1  mm  and  an  average  measured  resistivity  of  4200  Q-cm. 

Using  the  Drude  model,  the  predicted  index  of  refraction  would  be  3.4176  and  the 

absorption  coefficient  would  be  0.02  cm’1.  Figure  25a  and  25b  shows  the  pulse  and 
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spectrum  data  for  the  FZ  silicon  collected  with  the  THz  TDS  setup  with  the  first  Fabry- 
Perot  reflection  shown  at  the  far  right  of  the  time  domain  plot.  The  index  of  refraction  in 
25c  was  measured  to  be  3.416  and  the  absorption  coefficient  in  25d  was  approximately  0 
cm'1.  These  calculations  and  measurements  can  be  compared  to  those  made  by  Dai  et.  al. 
where  their  results  were  3.4175  for  the  index  of  refraction  and  0.01  -  0.23  cm'1  (0.5  -  4.5 
THz)  for  the  absorption  coefficient  for  FZ  silicon  (39  mm  thickness)  from  the  same 
manufacturer  [17]. 
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Figure  25.  Material  parameter  measurements  on  6  mm  thick  float  zone  (FZ)  silicon: 
(a)  THz  pulse,  (b)  THz  spectrum,  (c)  index  of  refraction,  and  (d)  absorption  coefficient. 

3.2.5  Fabry-Perot  Resonance  Removal  Techniques 


One  can  attempt  to  remove  the  Fabry-Perot  etalon  effect  illustrated  in  the  previous 
section  with  one  of  three  different  techniques.  One  such  technique  is  to  window  the  pulse 
in  the  time  domain,  removing  all  of  the  etalon  reflections  that  occur  after  the  pulse.  The 
problem  with  this  technique  is  that  the  reduction  in  the  total  collection  time  of  the  pulse 
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causes  the  resolution  of  the  THz  spectrum  in  the  frequency  domain  to  be  significantly 
reduced.  This  caused  an  inaccuracy  in  the  unwrapping  of  the  phase  change  required  to 
calculate  the  index  of  refraction. 

A  second  method  of  removing  the  Fabry-Perot  etalon  effect  is  to  subtract  out  the 
reflections  from  the  time  domain  pulse.  One  can  obtain  a  THz  field  in  the  time  domain 
that  has  removed  the  first  two  reflections,  assuming  that  there  is  little  or  no  absorption. 
In  the  following  equation  the  measured  THz  field  Em  (  t  )  is  multiplied  by  the  reflected 
field  at  the  air-silicon  interfaces  for  the  first  two  reflections  to  obtain  the  resultant  field 
E ,(!)-. 


E,(t)  =  Em(t)-REm(t+d)+R2Em(t  +  2d). 


(72) 


In  this  equation,  R  is  the  reflectance  and  dt  is  the  time  between  reflections,  given  by 

2nd 


R  = 


'n- \y 

v«  +  l  j 


,  dt  =  ■ 


(73) 


where  d  is  the  thickness  of  the  material.  One  can  also  observe  the  data  and  measure  the 
time  difference  between  the  main  pulse  and  the  reflected  pulses.  An  example  of  this  FP 
etalon  removal  technique  is  demonstrated  in  Figure  26.  A  similar  technique  was  used  by 
Naftaly  et.  al.  [52], 

Another  technique,  introduced  by  Duvillaret  et.  al.,  uses  an  iterative  process  to 
remove  the  Fabry-Perot  reflections  in  the  frequency  domain  [79].  A  block  diagram 
describing  the  technique  is  in  Figure  27.  The  first  iteration  starts  with  the  basic  equations 
to  determine  the  index  of  refraction  and  absorption  coefficient  from  section  3.2.1.  The 
Fabry-Perot  FP  term  represents  the  accumulation  of  all  the  etalon  reflections  through  the 
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Figure  26.  Diagrams  showing  the  removal  of  FP  oscillations  in  the  time  domain:  (a) 
index  of  refraction  with  FP  oscillations,  (b)  index  of  refraction  without  first  2  FP 
oscillations,  (c)  absorption  coefficient  with  FP  oscillations,  (d)  absorption  without 
first  2  FP  oscillations. 
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Figure  27.  Iterative  process  for  removal  Fabry-Perot  etalon  effect  in  frequency 
domain  described  by  Duvillaret  et.  al. 

sample  and  is  an  expansion  of  the  transmission  equation  from  section  2.5.5. 1 


E sam 

Eref  (ffl) 


This  Fabry-Perot  term  is  defined  as: 
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(75) 


r~ ,  .  -i2 

n[a>)- 1 

co 

L  xexp 

-i2n(co)-d 

n(a>)- l-l 

c 

where  the  index  of  refraction  is  a  complex  number,  including  the  absorption  coefficient. 

On  the  next  iteration,  one  divides  Esam/Eref  by  FP  and  then  re-calculates  the  optical 
parameters.  The  FP  term  is  then  readjusted  and  the  process  continues  until  the  optical 
parameters  stop  changing.  Results  are  shown  in  Figure  28  for  the  refractive  index  (a)  and 
(b)  and  the  absorption  coefficient  (c)  and  (d)  after  5  iterations  of  the  routine  described 
above.  In  addition  to  Fabry-Perot  removal  techniques,  efforts  have  been  made  to  remove 
noise  peaks  caused  by  water  vapor  absorption  lines  and  to  remove  spurious  effects  caused 
by  the  early  truncation  of  a  THz  signal  [28], 


Figure  28.  Fabry-Perot  oscillation  removal  technique  showing  the  real  part  of  the  index 
of  refraction  (a)  before  &  (b)  removal,  followed  by  the  imaginary  part  of  the  index  of 
refraction  (c)  before  &  (d)  after  the  Fabry-Perot  removal  technique  was  applied. 
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3.3  Reflection  Mode  Material  Parameter  Measurements 


Most  maintenance  techniques  limit  the  inspection  of  materials  from  only  one  side. 
Additionally,  some  materials  are  optically  thick  and  therefore  do  not  transmit  radiation  in 
the  THz  region.  A  different  setup  is  needed  to  characterize  these  types  of  materials.  The 
reflection  mode  THz  TDS  system  used  in  this  research  is  shown  in  Figure  29.  For  THz 
optics,  it  uses  a  25  mm  diameter  (f  =  25  mm)  parabolic  mirror  first,  followed  by  four  50 
mm  diameter  (f  =  100  mm)  parabolic  mirrors.  The  equation  used  to  model  reflection 
spectroscopy  at  normal  incidence  is  given  by 


^ s ample 

yjR(O)) 

exp[-iA<f>(co)]  [1  -  n(®)][  1  +  href  («)] 

Eref(o)) 

R,eM) 

exp[-z  A<f>ref  (of)]  [1  +  n(«j)\ [1  -  href  (of)] 

reference  mirror  must  be  positioned  at  precisely  the  same  location  as  that  of  the  sample. 
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Any  misalignment  in  the  reference  mirror  will  cause  a  severe  error  in  phase  and  should 
be  suppressed  to  less  than  1  pm.  The  misalignment  problem  in  reflection  measurements 
can  be  resolved  by  aligning  the  THz  TDS  system  with  attenuated-total-reflection  and 
ellipsometry  methods.  One  author  provides  dispersion  relations  in  the  analysis  of 
experimental  data  from  reflection  spectroscopy  to  help  remove  the  sample  misplacement 
error  [30].  A  second  paper  presents  the  effects  of  surface  roughness  on  the  reflection 
spectra  obtained  from  THz  TDS  [20].  A  final  paper  shows  a  technique  for  measuring  the 
alcohol  content  of  an  aqueous  solution  in  a  closed  bottle  using  a  reflection  geometry  [36]. 


3.3.1  Calculation  Technique 


If  one  uses  a  metallic  mirror  for  the  reference  scan  then  the  reflection  coefficient 
is  close  to  unity.  In  this  case,  the  equation  in  the  previous  section  reduces  to 


|r|exp(^)  =  -^— J-  (77) 

n  + 1 

in  which  r  is  the  amplitude  of  ratio  of  the  sample  spectrum  to  the  reference  spectrum  and 
the  relative  phase  is  ^  [35],  This  equation  can  then  be  inverted  and  simple  equations  can 
be  obtained  for  the  index  of  refraction  n  and  the  absorption  coefficient  a: 


1- 1 

r(co)2 

1+ 

r(&>)2  -2 

;  r(<u)  cos^(®) 

(78) 


4xo)  2|r(ffl)|sin^(®) 

2 

c  1+  r(a>)  -2  r(m)  cos^(<») 


(79) 
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A  correction  term  is  also  needed  to  account  for  the  angle  of  reflection  for  the  THz  TDS 
apparatus.  The  following  term  should  be  multipled  by  both  Equation  70  and  71  to 
account  for  the  incident  angle,  6 : 


/.Factor  =  cos  < 


sin  (6*) l/ 

sin  1 

/ n  . 

cos  (#) 

(80) 


Both  the  index  of  refraction  and  the  absorption  coefficient  have  terms  with  a 
phase  dependence.  In  addition,  the  relative  phase  would  have  to  unwrap  in  a  manner  in 

which  the  phase  would  stay  relatively  close  to  zero  so  that  cos(^) « 1  and  sin(^)  « (f) . 

Unfortunately,  in  practice  the  relative  phase  unwrapped  in  a  way  that  quickly  became 
greater  than  zero.  It  became  apparent  after  several  attempts  that  it  would  be  impossible, 
with  the  current  lab  equipment,  to  maintain  a  consistent  position  of  the  sample  and 
reference  and  therefore  the  required  phase  accuracy  could  not  be  obtained. 

An  alternate  approach  was  attempted  where  the  relative  phase  was  assumed  to 
change  at  a  much  slower  rate  by  dividing  the  relative  phase  by  105  so  it  remained  close  to 
zero,  but  would  still  unwrap  linearly.  Since  the  relative  amplitude  |r(<y)|  was  more 

consistent  across  the  THz  frequency  spectrum,  the  index  of  refraction  could  be 
approximated  using  equation  78.  The  index  of  refraction  for  FZ  silicon,  determined 
using  reflection  mode,  is  shown  in  Figure  30a  where  the  mean  value  is  3.40±0.04, 
assuming  an  angle  of  20±2°.  Although  the  mean  value  is  similar  to  the  transmission 
measurement,  the  figure  clearly  shows  a  greater  deviation  about  this  mean.  The 
absorption  coefficient  calculated  using  the  small  relative  phase  (unknown  accuracy)  is 
shown  in  30b. 
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Figure  30.  Material  parameter  measurements  using  reflection  configuration  on  6  mm 
thick  float  zone  (FZ)  silicon:  (a)  index  of  refraction  and  (b)  absorption  coefficient. 

3.3.2  Practical  Approximation  of  Material  Parameters 


Since  the  phase  accuracy  requirement  was  so  difficult  to  meet  when  comparing  a 
sample  to  a  reference  piece  of  aluminum,  an  alternate  approach  was  attempted.  Instead, 
the  reflection  off  the  back  surface  of  the  sample  was  compared  with  the  reflection  off  the 
front  side  of  the  sample.  An  assumption  was  made  that  the  index  of  refraction  spectrum 
had  a  flat  response  and  could  be  determined  via  a  time  of  flight  measurement.  Then  the 
absorption  coefficient  could  be  determined  using  this  index  of  refraction  measurement. 
The  easiest  way  to  test  this  method  was  by  using  the  thickest  piece  of  FZ  silicon  so  that 
the  effects  of  the  ZnTe  Fabry-Perot  reflections  would  not  interfere  with  the  back 
reflection.  The  results  of  the  pulse  are  shown  in  Figure  31a.  The  index  of  refraction  n 

was  calculated  using  the  following  formula:  n  =  {Th-Tf}c  2d  where  Tb  is  the  time  of 


arrival  of  the  back  pulse  and  Tf  is  the  time  of  arrival  of  the  front  pulse.  One  can  notice 
that  the  polarity  of  the  back  reflection  is  the  opposite  of  the  front  reflection  due  to  the  n 
phase  shift  encountered  when  reflection  from  the  first  air-to-glass  fiber  interface.  Using 
the  measured  thickness  d  of  6.0  mm,  the  calculated  value  of  n  was  3.418  which  was  very 
close  to  the  values  reported  for  transmission  configuration  and  for  the  Drude  model. 
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Figure  31.  THz  TDS  reflection  setup  for  thick  FZ  silicon:  (a)  pulses  off  front  and 
back  surface  and  (b)  absorption  coefficient  using  Fourier  transform  from  each  pulse. 

Since  the  pulses  were  sufficiently  separated  in  time,  one  could  split  the  200  psec 
time  period  in  half,  and  take  the  Fourier  transform  of  each  half.  The  absorption 
coefficient  is  then  calculated  like  in  transmission  setup,  except  that  the  Fresnel  reflections 
are  different.  A  diagram  showing  the  geometry  of  the  Fresnel  reflections  in  reflection 
setup  is  shown  in  Figure  32.  The  following  equation  shows  how  to  determine  the 
absorption  coefficient  a  from  the  two  reflected  pulses: 


a 


-In 


\Esam  (®)| 


The  term  R(co)  is  defined  based  on  the  geometry  in  Figure  32  as 


R(a>)  = 


^12  1  )r23  (^1  ) 


rn  l) 


where  the  coefficients  are  defined  according  to 


rn  = 


nx  cos  (0J  )-n2  cos  (#2 ) 
cos  (#, )  +  n2  cos  (#2 ) 


,  tl2  1  +  f\2  and 


'23 
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The  absorption  coefficient  could  then  be  calculated  as  shown  in  Figure  32b.  It  remained 
close  to  0  cm'1  and  shows  some  similarities  to  the  transmission  configuration  graph. 


Figure  32.  Diagram  showing  the  geometry  of  the  first  Fabry-Perot  reflection  from 
the  sample  in  reflection  configuration. 

3.3.3  Other  Material  Parameter  Measurements 

In  the  previous  two  sections,  the  thick  piece  of  FZ  silicon,  with  a  significant  time 
separation  before  the  first  FP  reflection,  was  used  to  test  two  methods  of  extracting 
material  measurements  using  a  THz  TDS  reflection  configuration.  In  this  section, 
material  parameter  measurements  were  attempted  on  the  high  resistivity  silicon  and 
finally  on  the  glass  fiber  in  an  effort  to  find  a  useable  reflection  technique. 

The  first  material  analysis  was  performed  on  the  high  resistivity  silicon.  Figure 
33  shows  the  THz  (a)  pulse  return  and  the  (b)  spectrum  with  the  effects  of  the  Fabry- 
Perot  reflections  from  both  the  sample  and  the  ZnTe.  This  made  it  difficult  to  determine 
the  index  of  refraction  from  Equation  78.  To  determine  the  material  parameters,  the  first 
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Figure  33.  THz  TDS  reflection  (a)  pulse  and  (b)  spectrum  for  1.016  mm  high 
resistivity  silicon.  Isolation  of  the  (c)  front  reflected  pulse  and  the  (d)  back  side 
reflected  pulse  in  the  time  domain.  The  (e)  absorption  coefficient  for  high  resistivity 
silicon  measured  using  the  spectrum  of  the  front  and  back  pulse. 

technique  that  was  attempted  was  to  measure  the  index  of  refraction  with  the  time  delay 
of  the  THz  pulse  between  the  front  and  back  reflections.  The  index  was  measured  to  be 
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3.405  from  either  the  peak  of  the  pulse  or  the  zero  crossing.  In  order  to  measure  the 
absorption  coefficient,  a  separate  Fourier  transform  would  have  to  be  taken  for  both  the 
front  and  the  back  reflection.  First,  the  front  pulse  was  multiplied  by  the  factor  from 
Equation  82  (-0.686)  and  then  subtracted  from  the  back  reflection  shown  in  33c.  Then 
the  reference  pulse  was  multiplied  by  the  Fresnel  reflection  coefficient  (0.546)  and 
subtracted  from  the  front  reflected  pulse  for  33d.  The  Fourier  transform  of  each  was 
taken  and  the  result  is  shown  in  33e.  The  result  oscillates  about  zero  from  0.1  -  1.2  THz. 

The  material  measured  next  was  the  glass  fiber  composite  without  an  exterior 
coating.  The  glass  fiber  was  unique  from  the  silicon  samples  in  that  it  had  a  considerably 
higher  absorption.  The  THz  pulse  in  Figure  34a  data  was  taken  with  the  TDS  reflection 
setup.  A  first  attempt  was  made  to  measure  the  index  of  refraction  using  the  timing 
between  the  front  and  back  reflections,  however,  the  index  measurements  were  lower 
than  in  transmission  configuration.  An  alternate  method  was  tried  using  the  relative 
amplitude  of  the  pulses  measured  peak-to-peak.  The  amplitude  of  the  first  glass  fiber 
pulse  (1.5133)  to  the  aluminum  reference  pulse  (4.2124)  was  measured  to  be  0.3592. 
Using  the  equation  for  the  Fresnel  reflection  coefficient  from  the  front  surface,  the  index 
of  refraction  was  calculated  to  be  2.03  ±  0.02.  With  this  index  of  refraction,  the 
absorption  coefficient  was  calculated  using  the  method  outlined  in  section  3.3.2.  The 
graph  in  34b  shows  the  Fabry-Perot  etalon  effect  in  the  glass  fiber,  while  the  graph  in  34c 
is  filtered  using  a  windowing  technique.  The  coefficient  was  measured  to  be  7  ±  2  cm'1 
at  0.5  THz  and  18  ±  2  cm'1  at  1.0  THz  which  were  close  to  the  results  for  transmission. 
The  measurement  appears  to  run  out  of  dynamic  range  at  approximately  1.0  THz,  which 
makes  it  slightly  worse  than  transmission  which  was  at  about  1 .25  THz. 
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(a)  (b) 
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Figure  34.  Plot  of  the  (a)  THz  pulses  reflected  from  the  front  and  back  of  a  glass 
fiber  coupon.  Absorption  coefficient  for  the  glass  fiber  (b)  without  filtering  and  (c) 
with  filtering. 

3.4  Burn  Diagnostics  Using  THz  TDS 

One  of  the  problems  with  designing  aircraft  is  that  the  heat  from  jet  engines  can 
cause  damage  to  the  external  structure  of  the  aircraft.  THz  TDS  could  be  used  to 
measure  the  material  properties  of  the  aircraft  glass  fiber  composites  to  determine  if  the 
composition  of  the  material  has  fundamentally  changed.  Polyimide  is  a  resin  which  has 
the  characteristic  of  being  very  resilient  to  thermal  loads,  maintaining  its  properties  for 
long  periods  of  continuous  use  at  230°C,  and  for  short  excursions  as  high  as  480°C. 


87 


Composite  samples  were  prepared  that  were  heated  with  various  temperatures  close  to  the 
maximum  for  polyimide  for  short  durations. 


3.4.1  Transmission  Mode 


One  of  the  first  steps  that  was  taken  in  the  nondestructive  evaluation  of  damaged 
samples  was  to  consider  whether  the  material  properties  were  altered  in  any  way  by  the 
damage.  Since  the  samples  with  heat  damage  on  them  were  different  than  the  previously 
measured  glass  fiber  samples  in  Figure  19,  we  measured  the  THz  pulse  and  amplitude 
spectrum  for  the  sample  and  reference  shown  in  Figure  35a  and  35b.  These  damage 
samples  had  a  coating  applied  to  the  exterior  of  the  glass  fiber  which  changed  the 
properties  slightly.  The  index  of  refraction  for  the  unaltered  composite  remained 
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Figure  35.  THz  TDS  (a)  pulse  and  (b)  amplitude  spectra  for  air  reference  and 
coated  glass  fiber,  (c)  Index  of  refraction  and  (d)  absorption  coefficient 
measurements  for  coated  glass  fiber  in  the  THz  frequency  range. 
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relatively  constant  in  the  THz  frequency  range  at  2.05  ±  0.01  in  35c  while  the  absorption 
coefficient  was  measured  as  17.5  ±  0.8  cm'1  (0.5  THz),  27.5  ±  0.8  cm'1  (0.8  THz)  and  38 
±  0.8  cm'1  (1  THz)  shown  in  35d.  This  is  comparable  to  another  measurement  made  on 
glass  fiber  samples,  where  the  index  of  refraction  was  2.05-2.1  and  the  absorption 
coefficient  was  12  cm'1  (0.5  THz)  and  26-  32  cm'1  (0.8  THz)  [4], 

Using  the  setup  shown  in  Figure  40,  we  measured  an  undamaged  glass  fiber 
sample  followed  by  a  series  of  damaged  samples  with  heat  damage.  Thickness 
measurements  were  made  on  each  of  the  samples  and  their  burn  areas,  where  the  most 
visually  noticeable  difference  was  in  the  largest  bum  area  (440°C  for  4  minutes).  Using 
the  measured  thickness  of  each  blistered  area,  we  obtained  the  results  in  Figure  36  for  the 
(a)  index  of  refraction  and  (b)  absorption  coefficient.  Figure  36  shows  that  there  is  a 
variation  based  on  the  thickness  of  the  sample  blister. 

For  comparison,  the  burned  material  parameter  measurements  were  calculated 
with  the  undamaged  sample  thickness,  such  that  the  optical  path  length  through  each 
damage  spot  was  constant.  As  is  shown  in  Figure  37,  two  of  the  bum  areas  have 


Figure  36.  THz  TDS  material  parameter  measurements  using  actual  thickness 
measurements  for  burn  damage  areas  showing  (a)  indices  of  refraction  and  (b)  absorption 
coefficients. 
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relatively  the  same  optical  path  length  as  the  undamaged  area.  The  third  damage  area, 
burned  at  430°C  for  6  minutes,  may  have  been  different  due  to  a  small  difference  within 
the  bum  spot.  There  is  an  area  in  the  center  of  the  bum  spot  where  there  appears  to  be 
either  a  removal  of  the  outer  coating  or  a  residue  created  during  the  burning  process. 
Since  it  is  in  the  center  of  the  small  bum  spot  it  was  not  possible  to  ensure  a  THz  TDS 
measurement  could  be  made  outside  of  the  difference  area.  These  measurements  suggest 
that  burning  the  samples  at  these  high  temperatures  does  not  alter  the  material 
parameters,  but  merely  introduces  a  blistering  of  the  glass  fiber  and/or  coating. 


Figure  37.  THz  TDS  material  parameter  measurements,  assuming  the  same 
thickness,  for  burn  damage  areas  showing  (a)  indices  of  refraction  and  (b) 
absorption  coefficients. 

3.4.2  Reflection  Mode 

Reflection  configuration  could  also  potentially  be  used  for  material  parameter 
bum  diagnostics  on  composites.  A  reflection  technique  would  assume  that  the  back  side 
reflection  of  the  THz  would  be  large  enough  to  use  the  technique  described  in  section 
3.3.3.  However,  an  exterior  coating  was  present  on  all  three  of  bum  spots  available  for 
this  research,  and  this  prevented  most  of  the  back  reflection  from  being  detected.  A 
comparison  of  a  reflected  pulse  from  glass  fiber  both  with  and  without  the  coating  is 
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shown  in  Figure  38.  Nearly  double  the  amplitude  is  reflected  from  the  coated  sample  due 
to  the  two  reflections  at  the  coating  and  composite  surfaces.  A  very  small  back  side 
reflection  is  evident  in  38a  for  the  coated  sample.  A  burned  area  of  the  glass  fiber 
without  a  coating  would  be  required  to  test  this  technique. 


Figure  38.  Comparison  of  the  reflected  signal  from  a  coated  and  uncoated  sample, 
showing  the  (a)  time  domain  and  (b)  frequency  spectrum  returns. 
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4.  THz  Imaging  for  Nondestructive  Evaluation 


This  chapter  explains  the  results  and  analysis  of  imaging  using  THz  TDS  for  the 
nondestructive  evaluation  of  composites.  It  is  divided  primarily  into  a  section  on 
transmissive  imaging  and  another  section  on  reflection  imaging.  However,  prior  to 
analyzing  the  imaging  results,  there  is  an  explanation  of  the  samples  and  the  setup, 
including  sections  on  the  numerical  modeling  and  measurements  of  the  THz  beam  at  the 
sample  position.  Using  each  of  the  two  imaging  techniques,  two-dimensional  imaging 
was  performed  on  the  test  samples.  Finally,  the  last  portion  of  each  technique  described 
the  analysis  of  the  depth  of  delaminations  using  THz  TDS. 

4.1  Sample  Preparation 

The  composite  samples  that  were  provided  by  AFRL/RXLP  to  test  representative 
forms  of  exterior  aircraft  damage  are  shown  in  Figure  39.  The  samples  that  were  heated 
with  various  temperatures  for  short  durations  are  labeled  2  and  3  in  Figure  39.  The  first 
sample  was  burned  at  440°C  for  4  minutes,  creating  a  blister  on  the  sample  about  2  cm  X 
1.5  cm.  The  second  sample  was  burned  in  two  places:  430°C  for  6  minutes  and  425°C 
for  20  minutes. 

Another  goal  was  to  determine  if  various  forms  of  damage  could  be  imaged  with  a 
THz  TDS  system  in  either  transmission  or  reflection  setup.  Three  additional  samples 
were  prepared  to  test  the  ability  of  a  THz  TDS  system  to  detect  the  damage.  The  first 
composite  specimen,  labeled  1  in  Figure  39,  was  used  as  a  thickness  standard  coupon  to 
test  the  response  of  THz  TDS  imaging  techniques  to  various  sample  thicknesses. 
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Figure  39.  Photograph  showing  the  5  glass  fiber  samples:  (1)  thickness  calibration 
sample,  (2)  &  (3)  burn  samples,  (4)  mechanical  stress  sample,  and  (5)  hidden  defect 
sample  showing  the  hidden  location  of  two  of  the  eight  defects. 

Another  coupon,  labeled  5,  was  prepared  to  search  for  flaws.  This  coupon  consisted  of 

two  laminated  pieces  and  is  described  in  another  reference  [50].  The  bottom  piece  was 

prepared  by  etching  four  layers  and  then  creating  a  6  mm  slit  and  a  3  mm  flat  bottom  hole 

(—70  pm  depth)  in  each  of  the  four  layers.  A  top  layer  was  created  by  etching  4  layers 

with  the  same  thicknesses  and  then  attaching  to  the  bottom  layer  in  the  opposite  direction 

with  epoxy,  being  careful  not  to  allow  any  epoxy  into  the  flaws.  Finally,  a  third  sample, 

labeled  4  in  Figure  39,  was  prepared  by  bending  a  piece  of  glass  fiber  about  a  fixed  axis  a 

total  of  6240  cycles  to  investigate  mechanical  damage  with  THz  TDS  imaging. 


4.2  Transmission  Mode  Setup 

The  setup  that  was  used  for  the  THz  TDS  nondestructive  evaluation  and  imaging 
is  illustrated  in  Figure  40.  A  photoconductive  switch  was  used  for  the  THz  emitter  and 
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an  electro-optic  crystal  (ZnTe)  was  used  for  the  detector.  This  setup  differed  from  the 
one  used  for  material  parameter  measurement  in  that  it  used  four  off-axis  parabolic 
mirrors  to  steer  the  THz  beam,  focusing  it  tightly  at  the  sample.  The  first  and  last 
parabolic  mirror  had  a  diameter  and  focal  length  of  25  mm  while  the  center  two  parabolic 
mirrors  had  a  focal  length  and  diameter  of  50  mm. 


Figure  40.  THz  TDS  Setup  for  Imaging  Composites  in  Transmission  Mode. 


4.3  THz  Beam  Propagation  and  Modeling 

This  section  shows  the  results  of  modeling  a  THz  beam  spatially  as  a  Gaussian 
through  our  THz  TDS  system  and  the  minimization  of  the  spot  size  at  the  sample.  The 
numerical  results  are  presented  for  both  the  transmission  and  reflection  configurations. 
Finally,  actual  THz  spot  size  measurements  are  compared  to  the  modeling  predictions. 
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4.3.1  Spot  Size  Minimization 


The  first  step  in  calculating  the  propagation  of  a  THz  beam  through  an  optical 
system  was  to  obtain  the  measured  radius  of  curvature  and  spot  size  measurements  at  the 
output  of  the  hemispherical  lens  from  Ekspla,  the  THz  emitter  and  detector  manufacturer. 
These  parameters  are  listed  in  Table  7  below.  The  Ekspla  hemispherical  lens  is 
considered  a  collimating  lens  and  the  beam  actually  narrows  initially  as  it  propagates  out 
of  the  lens  before  it  spreads.  Many  of  the  hemispherical  designs  result  in  the  spreading  of 
the  THz  beam  at  a  given  angle,  which  can  then  be  collimated  by  a  parabolic  mirror  and 
then  focused  by  a  second  parabolic  mirror. 


Table  7.  Ekspla  measurements  at  the  exit  of  the  silicon  hemispherical  lens. 


Frequency 

Radius  of  curvature 

Spot  size 

(THz) 

(mm) 

(mm) 

1.0 

-0.224 

4.1641 

2.0 

-0.225 

4.16395 

3.0 

-0.2256 

4.1639 

The  transmission  mode  setup  for  this  experiment  requires  a  unique  mirror 
focusing  arrangement  to  achieve  the  diffraction-limited  spot  size  at  the  sample.  One  has 
to  ensure  that  the  THz  beam  fills  the  entire  second  mirror,  is  minimized  at  the  sample, 
does  not  expand  beyond  the  diameter  of  any  of  the  mirrors,  and  then  is  focused  at  the 
detector.  A  diagram  showing  the  propagation  of  the  THz  beam  through  the  optical 
system  is  shown  in  Figure  41.  The  final  design  is  able  to  achieve  the  diffraction- limited 
spot  size  within  a  limited  path  length  and  with  available  equipment.  A  summary  of  the 
diffraction- limits  and  the  calculated  values  for  the  optical  design  is  in  Table  8.  The  beam 
diameter  at  the  focus  (*)  at  0.5  THz  was  estimated  by  extrapolating  the  radius  of 
curvature  and  spot  size  from  the  Ekspla  data. 
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Figure  41.  Diagram  showing  the  transmission  optical  system  with  THz  beam 
focusing  on  sample  at  1  THz. 


Table  8.  Comparison  of  diffraction-limited  spot  sizes  with  THz  beam  optical 
design  focus  spot  sizes  for  the  transmission  configuration. _ 


Frequency  (THz) 

Diffraction-Limited 
Beam  Diameter  (pm) 

THz  Beam 
Diameter  at  Focus 
for  Design  (pm) 

0.5 

1200 

1181* 

1.0 

600 

617 

2.0 

300 

312 

3.0 

200 

209 

A  similar  analysis  can  be  performed  to  design  the  optical  system  for  a  THz  TDS 
reflective  configuration.  The  goal  of  that  design  is  to  minimize  the  spot  size  at  the 
sample  with  available  equipment  within  the  shortest  path  length.  Using  the  Ekspla  data 
from  Table  7,  one  can  design  an  optical  system  in  which  five  parabolic  mirrors  are  used 
to  steer  the  THz  beam  as  shown  in  Figure  42.  The  maximum  diameter  of  the  THz  beam 
at  the  third  parabolic  mirror  is  3 1  mm,  limiting  the  focusing  ability  of  the  system.  The 
diffraction-limited  spot  size  and  the  calculated  Gaussian  beam  size  at  the  sample  are 
shown  in  Table  9. 
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Figure  42.  Diagram  showing  the  reflection  optical  system  with  THz  beam  focusing  on 
sample  at  1  THz. 


Table  9.  Comparison  of  diffraction-limited  spot  sizes  with  THz  beam  optical 
design  focus  spot  sizes  for  the  reflection  configuration. _ 


Frequency  (THz) 

Diffraction-Limited 
Beam  Diameter  (mm) 

THz  Beam 
Diameter  at  Focus 
for  Design  (mm) 

0.5 

2.42 

2.48* 

1.0 

1.21 

1.27 

2.0 

0.60 

0.69 

3.0 

0.40 

0.52 

4.3.2  THz  Beam  Measurements 

The  next  task  was  to  build  and  align  the  THz  TDS  setup  similar  to  the  numerical 
calculations,  tweaking  the  parabolic  mirrors  and  hemispherical  mirrors  to  minimize  the 
spot  size.  Although  the  final  setup  did  not  follow  the  model  exactly,  it  provided  a  good 
starting  point.  The  first  setup  that  was  measured  for  spot  size  at  the  sample  was  the 
transmission  configuration.  The  first  technique  that  was  used  was  the  knife  edge  method, 
scanning  a  knife  blade  through  the  THz  beam  horizontally  and  vertically.  The  knife  was 
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scanned  through  the  beam  from  0%  blockage  to  100%  blockage  and  beamwidth  was 
measured  between  88%  of  peak  voltage  and  12  %  of  peak  voltage  of  an  unobstructed 
beam  with  the  electro  optic  detector  [21].  The  smallest  beam  measurements  were  1.65 
mm  in  the  horizontal  direction  and  2.38  mm  in  the  vertical  direction,  making  it  elliptical 
in  nature.  A  circular  hole  of  roughly  1.5  mm  in  diameter  was  then  moved  into  the 
position  of  the  minimum  beam  size.  The  result  was  that  approximately  90%  of  the  peak 
pulse  voltage  was  transmitted  through  the  hole. 

Before  imaging  the  composites  in  reflection  configuration,  the  spot  size  was 
measured  again.  The  same  knife  edge  technique  was  attempted  in  reflection 
configuration,  except  that  the  knife  was  used  to  reflect  the  THz  and  not  block  it.  In  the 
horizontal  direction,  a  spot  width  of  3.0  mm  was  obtained,  while  a  spot  width  of  2.75  mm 
was  obtained  for  the  vertical  direction.  The  positioning  of  the  knife  was  more  difficult  in 
reflection  setup  so  another  method  was  attempted.  This  second  method  included  the  use 
of  the  circular  hole  with  a  piece  of  aluminum  behind  the  hole.  The  geometry  of  this  setup 
is  shown  in  Figure  43  where  0  is  the  angle  of  incidence,  X  is  the  diameter  of  the  incident 
THz  beam,  and  Y  is  the  spot  size  of  the  THz  beam  on  the  aluminum  reflector.  The 
method  of  calculating  the  spot  size  on  the  sample  (Y)  was  to  close  the  pupil  and  the 
measure  the  diameter  at  which  92%  of  the  pulse  peak  was  detected.  The  following 
equation  was  used  to  determine  the  spot  size  (Y): 

Y  =  D-2d  sin(<9)  (85) 

where  d  =  5  mm.  Given  an  incident  angle  of  20°  and  a  measured  pupil  diameter  of  7.2 
mm,  the  spot  size  was  calculated  to  be  3.78  mm.  Since  it  was  at  an  angle,  the  actual  THz 
beam  width  would  be  around  3.55  mm. 
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Figure  43.  Diagram  showing  the  geometry  of  the  incident  and  reflected  beam  off  the 
apparatus  for  measuring  the  THz  spot  size  in  TDS  reflection  configuration. 

4.4  Imaging  of  Composites  with  Transmission  Setup 

This  section  shows  and  explains  the  results  for  the  first  method  of  THz  imaging 
using  the  TDS  transmission  setup  shown  in  Figure  40  to  obtain  the  images.  Although  this 
setup  is  not  likely  for  maintenance  operations,  it  provides  insight  into  the  capabilities  and 
limitations  of  THz  TDS  imaging  in  general.  In  the  first  part  of  this  section,  there  is  a 
description  of  some  of  the  results  using  thickness  measurements  and  a  comparison  of 
some  image  processing  techniques.  The  next  portion  of  the  section  shows  the  results  of 
THz  TDS  scans  on  various  types  of  damage,  including  heat  damage,  bending  damage, 
and  hidden  defects.  Finally,  the  last  section  shows  a  measurement  analysis  of  the  depth 
of  delaminations. 

4.4.1  THz  TDS  Thickness  Measurements 

One  area  of  interest  for  imaging  is  whether  the  characteristics  of  the  pulse  change 
with  different  thicknesses  and  heat  damaged  areas  of  glass  fiber  material.  Figure  44 
shows  how  the  THz  pulse  is  not  only  reduced  in  amplitude  through  thicker  samples,  but 
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is  also  spread  out  in  time  as  higher  frequency  components  are  attenuated  more  rapidly 
than  lower  ones.  The  arrival  time  of  the  pulse  is  also  longer  due  to  a  longer  optical  path 
length  through  thicker  materials.  This  illustrates  how  the  THz  TDS  system  could  be  used 
to  measure  sample  thickness  and  hence  a  simple  phase  technique  could  be  used  to  image 
composites  for  nondestructive  evaluation. 


Figure  44.  Relative  delay  of  THz  pulses  for  various  thicknesses  of  glass  fiber 
sample,  including  those  that  were  etched  (*). 

4.4.2  Image  Processing  Comparison 

The  image  processing  technique  that  was  described  in  section  2.4.6  was  tested  on 
a  piece  of  glass  fiber  with  a  circle  etched  in  it.  The  circle  was  12.5  mm  in  diameter  and 
approximately  70  pm  in  depth.  A  series  of  images  are  shown  in  Figure  45,  showing  the 
progression  of  techniques  from  the  original  1  X  1  mm  pixel  image  in  a.  Figure  45c 
shows  the  image  after  zero-padding  in  the  spatial  frequency  domain  and  Figure  45d 
shows  the  additional  of  the  MATLAB  ‘interp’  function  to  further  smooth  out  the  image  in 
Figure  45c.  Figure  45b  was  added  to  show  the  effect  of  the  MATLAB  ‘interp’  function 
without  the  frequency  domain  interpolation.  The  image  in  Figure  45e  is  the  equivalent  of 
viewing  the  original  pixel  image  in  (a)  at  1/10  the  original  size. 
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Figure  45.  THz  image  of  a  12.5  mm  diameter  circle  etched  in  glass  fiber  showing  (a) 
the  original  1X1  mm  pixel  image,  (b)  image  after  MATLAB  finterp’  function,  (c) 
image  after  zero-padding  10  times  in  the  spatial  frequency  domain,  and  (d)  final 
image  after  applying  MATLAB  finterp’  function  to  part  (c).  (e)  Comparison  of  (c) 
with  l/10th  scale  pixel  image. 
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One  of  the  consequences  of  using  the  zero-padding  approach  was  that  certain 
features,  like  sharp  edges,  caused  a  periodic  modulation  in  the  image.  An  example  of  this 
effect  can  be  seen  in  Figure  46a.  The  “checker  board”  pattern  could  be  reduced  by  zero¬ 
padding  with  fewer  zeros  and  then  using  a  Wiener  filter  on  the  image  to  reduce  the 
modulation  as  is  shown  in  Figure  46b. 


0  2  16  8  10  Peak  Pulse  0  2  4  0  S  10  Peak  Pulse 


ititti  Position  (psec)  mm  Position  (psec) 

Figure  46.  Image  comparison  showing  the  (a)  the  zero-padding  technique  (10  times) 
and  the  (b)  zero-padding  technique  (4  times)  with  the  addition  of  a  Wiener  filter. 

4.4.3  2-D  Imaging  Results 

The  first  sample  was  prepared  to  calibrate  the  system  based  on  material  thickness. 
A  THz  image  was  taken  of  a  1  X  1  cm  section  of  overlapping  milled  out  areas  with 
different  thicknesses.  The  image  is  shown  in  Figure  47  where  (a)  represents  the  peak 
pulse  amplitude  and  (b)  shows  the  pulse  position  at  each  pixel  of  the  sample  shown  in  (c) 
[48].  The  decrease  in  amplitude  of  the  pulse  at  the  edges  is  a  result  of  THz  scattering 
[48],  The  periodic  modulation  in  Figure  47b  is  a  result  of  the  interpolation  technique  to 
smooth  out  the  pixels  and  does  not  match  the  pattern  of  the  glass  weave  within  the 
composite. 
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Similar  images  can  be  shown  using  the  area  under  the  curve  from  the  amplitude 
spectrum  in  the  frequency  domain.  The  image  in  Figure  48  shows  pixels  that  represent 
the  sum  of  the  amplitudes  over  the  entire  frequency  spectrum  that  could  penetrate  the 
glass  fiber  composite.  It  closely  resembles  the  peak  pulse  amplitude  image  shown  in 
Figure  47a.  The  differing  heights  of  the  edges  of  the  milled  out  areas  causes  a  frequency 


Figure  47.  THz  TDS  transmission  images  showing  a  section  of  the  glass  fiber  that  had 
been  milled  to  two  different  thicknesses  using  peak  pulse  amplitude  (a)  and  peak 
pulse  position  (b)  techniques.  A  picture  of  the  scan  area  on  the  calibration  sample  (c). 

dependence  in  the  amplitude  intensity.  This  can  be  seen  in  Figure  48  (b)  and  (c).  The 
upper  left  quadrant  is  about  0.92  mm  lower  than  the  right  half  while  the  lower  left 
quadrant  is  about  0.69  mm  lower  than  the  right  half.  The  ridge  between  the  two  left 
quadrants  has  a  height  difference  of  0.23  mm.  The  THz  spot  size  is  wavelength 
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dependent  and  was  spread  across  the  milled  ridges  differently  depending  on  the  position 
of  the  pixel  during  this  raster  scan.  The  lower  amplitude  areas  represent  the  positioning 
of  the  THz  spot  so  that  a  majority  of  the  radiation  at  that  frequency  was  scattered. 


Figure  48.  THz  images  formed  by  summing  the  area  under  the  curve  of  the 
amplitude  spectrum  within  a  given  frequency  range  for  each  pixel:  (a)  0.1  -  1.3 
THz,  (b)  0.5  -  0.7  THz,  (c)  0.7  -  0.9  THz. 

4.4.3. 1  Imaging  of  Burn  Damage 


The  next  task  was  to  image  the  bum  samples.  The  results  of  the  bum  dot  images 
formed  by  the  peak  amplitude  of  the  THz  pulses  are  shown  in  Figure  49  with  the  visibly 
burned  areas  within  the  circles.  Figure  49a  shows  the  sample  that  was  burned  at  440°C 
for  4  minutes  over  ~  2  X  1 .5  cm  area.  It  has  a  visibly  noticeable  bubble  or  blister  on  its 
exterior  which  roughly  corresponds  with  the  blue  area  within  the  black  oval. 
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Figure  49.  THz  TDS  images  and  photos  for  three  burn  areas  on  glass  fiber  samples: 
(a),  (d)  440°C  for  4  minutes;  (b),  (e)  430°C  for  6  minutes;  and  (c),  (f)  425°C  for  20 
minutes. 
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Comparisons  of  the  THz  time  domain  signal  of  pixels  from  the  bum  area  were  made  with 
those  from  outside  the  bum  area.  A  comparison  is  shown  in  Figure  50  contrasting  a  THz 
signal  from  an  undamaged  portion  of  the  sample  to  one  from  the  center  of  the  large  bum 
area  (440°C  for  4  minutes).  There  appeared  to  be  time  domain  reflections  in  several  of 
the  pixels  within  the  bum  area,  showing  evidence  of  air  gaps,  but  without  a  consistent 
pattern  between  pixels.  These  reflections  caused  a  reduction  in  the  peak  pulse  amplitude 
based  on  the  Fabry-Perot  etalon  effect  from  the  air  -  composite  interfaces,  which  reduces 
the  initial  peak  pulse,  and  because  of  scattering  from  the  small  air  gaps  and 
discontinuities.  In  the  frequency  domain,  there  were  no  consistent  spectral  changes  to  the 
THz  signals  through  the  bum  areas. 


0  2  4  6  S  10  Peak  Pulse 

mni  Amplitude 


Figure  50.  THz  TDS  (a)  pulse  image  and  (b)  amplitude  spectra  for  undamaged  glass 
fiber  sample  and  for  an  area  with  burn  damage  (440°C  for  4  minutes),  (c)  THz  TDS 
amplitude  image  of  high  resistivity  silicon. 
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The  other  two  bum  areas  are  shown  in  Figure  49b  and  49c,  neither  of  which 
showed  much  visual  evidence  of  blistering.  The  dark  blue  area  within  the  circle  in  Figure 
49b  is  roughly  equivalent  to  the  position  of  the  residue  or  coating  loss.  Since  the 
absorption  spectrum  remained  higher  than  the  unbumed  sample,  even  when  adjusted  to 
the  same  optical  path  length,  it  is  more  likely  that  the  dark  blue  area  is  a  form  of  residue. 
The  red  circular  dots  in  Figure  49c  and  less  noticeable  orange  dots  in  Figure  49b 
correspond  to  the  white  marker  dots  made  on  the  samples  to  show  the  extent  of  the  bum 
area.  The  THz  image  of  the  damage  area  in  Figure  49c  was  inconclusive  in  showing 
evidence  of  burning. 

During  many  of  the  transmission  scans,  the  laser  output  power  varied,  and  this 
variation  was  not  stabilized  until  after  the  laser  optics  were  cleaned.  To  demonstrate  this 
variation,  a  scan  was  made  of  the  thin  wafer  of  high  resistivity  silicon.  This  scan  can  be 
seen  in  Figure  50c  showing  the  approximate  color  variation  for  the  pulse  amplitudes  that 
were  demonstrated  in  49c.  The  mean  amplitude  was  measured  at  1.306  mV  while  the 
standard  deviation  was  0.017  mV  (1.33%)  and  the  maximum  variation  was  ±4%  from  the 
mean.  Over  the  same  scan  period,  the  laser  power  varied  by  3%,  roughly  matching  the 
pulse  amplitude  pattern  shown  in  the  image.  If  one  considers  that  in  Figure  49c  that  the 
amplitude  varied  by  more  than  ±7%  across  the  image  (not  including  the  marker  dots)  then 
one  could  say  that  the  signal-to-noise  of  the  transmission  setup  contributed  to  the  image 
variations,  but  was  not  the  most  significant  factor. 
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4.4.3.2  Imaging  of  Bending  Damage 


An  attempt  was  made  to  show  damage  caused  by  mechanical  fatigue  as  a  result  of 
6240  bending  cycles.  Visually,  one  could  observe  a  thin  area  of  discoloration  on  one  side 
of  the  glass  fiber  strip  and  one  could  see  a  small  amount  of  cracking  and  buckling  on  the 
back  side.  In  Figure  51,  the  image  of  the  sample  shows  an  area  of  lesser  amplitude 
corresponding  roughly  to  where  the  axis  of  bending  occurred.  Electrical  tape,  in  the 
shape  of  an  ‘X’,  was  attached  at  the  top  of  the  image  area  for  reference  as  shown  in  the 
first  photo  in  5  lb.  The  second  photo  in  5  lc  shows  the  location  of  cracking  on  buckling 
on  the  back  side  of  the  sample  that  corresponds  to  the  pattern  in  the  THz  image.  The 
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Figure  51.  THz  TDS  image  showing  bend  damage  across  the  central  bend  axis  (a). 
Photographs  of  the  (b)  front  side  of  the  composite  strip  showing  the  scan  area  with  the 
location  of  the  taped  6X’  and  the  (c)  back  side  showing  the  location  of  the  cracking  and 
buckling  and  the  hexagonal  structure,  (d)  Ultrasound  image  of  the  bend  damage. 
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cracking  and  buckling  follows  the  hexagonal  pattern  faintly  visible  in  the  picture. 
Finally,  the  image  in  5  Id  shows  an  ultrasound  image  of  the  same  damage,  showing  the 
same  damage  area  and  the  electrical  tape. 


4.4.3.3  Imaging  of  Hidden  Defects 


Voids  were  also  investigated,  which  simulated  either  manufacturing  defects  or 
damage  caused  by  stress  over  time.  THz  images  are  shown  in  Figure  52a  and  52b  for  a 


circular  void  (3X3  mm)  and  a  crack  or  slit  void  (6  mm  length)  in  52c.  For  the  circular 
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Figure  52.  THz  TDS  images  showing  3  mm  diameter  milled  area  hidden  between 
two  glass  fiber  strips  using  (a)  peak  pulse  amplitude  and  (b)  peak  pulse  position. 
Linear  slit  void  (6  mm  length)  (c)  also  hidden  between  two  strips  of  glass  fiber.  The 
final  picture  is  an  x-ray  computed  tomography  scan  of  the  laminated  sample 
showing  the  hidden  circular  and  slit  voids  in  the  sample. 
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void,  a  simple  time  domain  amplitude  or  phase  technique  was  sufficient  to  detect  the 
void.  For  the  slit  void,  a  specific  frequency  range  was  required  to  isolate  the  approximate 
position,  and  was  more  difficult  to  isolate  on  multiple  attempts  [48].  The  voids  show  an 
area  of  reduced  amplitude,  most  likely  due  to  the  multiple  reflections  from  the 
air/composite  interface.  A  larger  number  of  samples  per  area  could  be  useful  in  isolating 
this  slit  void  in  space  by  providing  more  complete  coverage  at  the  higher  frequency 
components.  Another  candidate  would  be  to  use  a  rotate  the  sample  at  45°  to  the  THz 
polarization  and  then  collect  pulses  with  the  detector  at  horizontal  polarization  and  then 
vertical  polarization. 

4.4.4  Depth  of  Delamination  Analysis 

The  previous  experiments  demonstrated  the  imaging  of  voids  in  a  two- 
dimensional  plane.  An  attempt  was  made  to  use  THz  TDS  to  isolate  a  void  in  the  third 
dimension  of  depth  using  the  time  domain.  The  goal  was  to  use  THz  TDS  to  find  Fabry- 
Perot  reflections  in  a  delaminated  area  of  the  laminated  sample  (#5  in  Figure  39)  where 
the  epoxy  had  not  adhered  the  two  glass  fiber  strips.  The  THz  scans  are  shown  in  Figure 
53  a,  showing  the  THz  pulse  after  it  had  transmitted  through  the  adhered  area  and  the 
delaminated  area.  The  autocorrelation  was  taken  for  each  of  the  two  pulses 
independently  of  each  other  in  the  attempt  to  show  the  Fabry-Perot  reflections  in  Figure 
53b.  Then  the  sample  was  flipped  over  and  the  process  was  repeated  in  the  opposite 
direction  Figure  53c  and  53d.  The  timing  of  the  Fabry-Perot  reflections  was  calculated 
with  TFP_re/  =  InJ j c  where  is  the  index  of  the  glass  fiber  and  /  is  the  thickness  of  a 
delaminated  piece.  Given  that  the  index  of  refraction  is  2.05  and  that  the  thickness  of  the 
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two  layers  is  approximately  0.23  mm  and  0.92  mm,  the  predicted  times  would  be  3.12 
psec  and  12.5  psec.  The  first  Fabry-Perot  reflection  occurred  at  approximately  3.3  psec 
after  the  main  pulse  in  the  autocorrelations.  The  approximate  location  of  the  second 
Fabry-Perot  reflection  is  shown;  however,  the  autocorrelation  is  difficult  to  identify  on 
the  graph.  The  Fabry-Perot  reflections  show  up  similarly  in  both  directions,  indicating 
the  ambiguity  in  determination  of  depth  using  a  transmissive  setup. 


Figure  53.  THz  pulses  after  propagating  through  (a)  laminated  and  delaminated 
portions  of  a  glass  fiber  strip,  (b)  Autocorrelation  of  each  of  the  two  pulses  showing 
the  approximate  location  of  Fabry-Perot  reflections.  Thz  pulse  propagation  in  the 
opposite  direction  showing  pulses  and  their  autocorrelations  (c),  (d). 

The  attenuation  of  the  THz  signal  can  be  estimated  based  on  the  measured 
absorption  coefficient  of  the  glass  fiber,  TFP  (&>)  =  R(co)exp(-a(oj)l)  ?  where  R(co)  is  the 
reflection  of  the  signal  off  of  the  two  air-sample  interfaces.  The  amount  of  terahertz 
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radiation,  relative  to  the  initial  pulse,  that  is  transmitted  through  the  glass  fiber  composite 
after  undergoing  a  single  Fabry-Perot  reflection  is  shown  in  Figure  54,  with  a  diagram  of 


the  two  pulses  in  54a,  calculations  for  TFP(co)  in  54b,  and  a  picture  of  the  0.23  mm 
delamination  in  54c.  One  can  observe  in  the  figure  that  unless  a  delamination  is  very  thin 
(~  0.23  mm)  the  reflection  is  unlikely  to  be  recognizable  above  the  system  response. 


Figure  54.  (a)  Diagram  showing  the  main  pulse  and  the  first  Fabry-Perot  reflection 
in  transmission  configuration,  (b)  A  chart  showing  the  relative  amplitude  of  the 
first  Fabry-Perot  reflection  to  the  main  pulse  after  traveling  through  various 
thicknesses  of  composite,  (c)  A  photograph  showing  the  delamination. 
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4.5  Reflection  Mode  Imaging  of  Composites 


This  section  shows  the  THz  imaging  results  collected  using  the  reflection  setup 
shown  in  Figure  29.  The  reflection  configuration  is  a  more  likely  technique  for 
maintenance  than  transmission.  The  samples  were  the  same  as  those  used  for  the 
transmission  setup,  providing  a  direct  comparison  of  the  two  techniques.  This  section 
also  looks  at  the  same  types  of  damage  in  two-dimensional  imaging  as  well  as  a  depth  of 
delamination  analysis.  Finally,  results  are  shown  for  reflection  scans  on  entire  exterior 
panels.  All  of  the  reflection  mode  images  are  taken  with  1X1  mm  step  size  in  the  raster 
scan. 

4.5.1  Thickness  Calibration 

A  comparison  of  the  thickness  calibration  sample  imaged  with  a  reflection  setup 
was  compared  to  the  transmission  setup.  The  results  of  the  peak  pulse  amplitude  and 
peak  pulse  position  images  are  shown  in  Figure  55a  and  55b.  The  resolution  of  the 
images  is  not  as  good  as  the  transmission  images,  but  they  still  clearly  show  the  three 
distinct  areas.  The  scattering  of  the  THz  signal  is  also  shown  in  dark  blue  in  Figure  55a 
along  the  sharp  edge  between  the  different  milled  areas.  The  difference  in  amplitude  in 
each  area  is  based  on  the  geometry  of  the  THz  optical  system  and  not  based  on  any 
characteristic  of  the  glass  fiber.  The  phase  image  in  Figure  55b  shows  the  opposite 
pattern  of  the  pulse  arrival  from  the  images  for  transmission  configuration.  Thicker 
sections  have  a  shorter  arrival  time  in  reflection,  but  a  longer  arrival  time  in  transmission. 
An  additional  scan  was  made  by  lengthening  the  pixel  scan  time  from  25  psec  to  40  psec 
to  ensure  that  the  pulse  reflected  off  the  back  side  of  the  glass  fiber  was  always  collected. 
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From  this  data,  a  thickness  image  Figure  55c  could  be  assembled  using  the  optical  path 
length  of  the  THz  in  glass  fiber. 
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Figure  55.  THz  TDS  reflection  images  showing  a  section  of  the  glass  fiber  that  had 
been  milled  to  two  different  thicknesses  using  peak  pulse  amplitude  (a)  and  peak 
pulse  position  (b)  techniques.  The  final  image  (c)  shows  a  thickness  estimation  for  a 
separate  scan. 

4.5.2  2-D  Imaging  Results 

The  following  section  shows  the  image  results  for  the  various  types  of  damage. 
One  of  the  issues  that  makes  reflection  mode  imaging  different  from  transmission  setup  is 
that  the  tip  and  tilt  of  the  sample  can  vary  as  the  sample  is  moved  through  the  THz  beam. 
This  is  caused  by  a  slight  angle  misalignment  (±2°)  of  the  sample  to  the  raster  scanning 
mechanism  and/or  a  warping  or  curvature  of  the  sample.  This  can  create  an  artificial 
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delay  of  the  pulse  from  top-to-bottom  or  from  left-to-right.  It  can  have  the  same  effect  on 
the  amplitude  due  to  a  changing  of  the  geometric  alignment  of  the  reflected  beam  through 
the  collection  optics.  Linear  correction  techniques  were  implemented  in  post  processing 
to  attempt  to  correct  these  deficiencies.  The  bandwidth  of  the  signal  reflected  from  the 
surface  was  much  greater  than  the  bandwidth  of  any  signal  that  penetrated  the  glass  fiber 
at  any  depth. 

4.5.2.1  Imaging  of  Burn  Damage 

The  bum  damage  was  imaged  using  the  reflection  configuration  on  the  three  areas 
analyzed  earlier.  The  first  bum  area  (440°C  for  4  minutes),  which  had  the  most 
noticeable  blister,  was  imaged  first  and  the  amplitude  and  phase  images  are  shown  in 
Figure  56.  The  rough  area  of  the  blister  is  encompassed  by  the  oval.  The  amplitude 
image  shows  the  sensitivity  of  the  geometry  of  the  sample  within  the  TDS  setup.  One 
part  of  the  bubble  actually  reflected  more  THz  than  the  undamaged  portion  due  to  the 
initial  alignment  with  the  upper  right  comer  of  the  image.  The  angular  dependence  of  the 
reflection  to  the  sample  geometry  causes  a  variation  in  the  reflected  amplitude  due  to  a 
changing  alignment  of  the  THz  beam  within  the  THz  optics.  The  phase  plot  shows  a 
more  accurate  representation  of  the  bubble  where  the  true  shape  of  the  bubble  is  more 
apparent. 

One  can  also  observe  the  pulse  shape  characteristics  in  each  of  the  pixels.  For 
example,  a  pixel  was  chosen  on  the  right  side  of  the  bubble  and  is  shown  in  Figure  56c. 
One  can  see  the  double  reflection  off  the  bubble  surface  and  off  a  hidden  delamination 
below.  The  time  of  flight  difference  of  6.6  psec  can  be  used  to  calculate  the  distance 
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between  the  pulses.  Assuming  an  index  of  refraction  of  2.05,  the  distance  between  the 
reflections  was  483  pm.  The  index  is  a  rough  estimate,  because  the  bubble  could  contain 
some  combination  of  air,  coating,  and  glass  fiber.  However,  if  this  estimate  is  close,  then 
it  would  make  sense  that  the  delamination  is  close  to  where  the  original  surface  was 
located.  An  x-ray  image  of  the  sample,  showing  the  bubble  from  the  side  profile,  can  be 
seen  in  Figure  56d. 
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Figure  56.  THz  TDS  reflection  images  for  burn  damage  from  heating  at  440°C  for  4 
minutes  for  (a)  peak  pulse  amplitude  and  (b)  peak  pulse  position.  Pulse  return  from 
the  right  side  of  the  bubble  showing  the  double  pulse  return  (c).  X-ray  computed 
tomography  image  showing  a  side  profile  of  the  large  burn  blister  (d). 

Additionally,  one  can  determine  the  approximate  height  of  the  bubble  by 
measuring  the  difference  in  the  pulse  phase  between  the  center  of  the  bum  area  and  the 
unbumed  area.  From  the  diagram,  one  could  choose  a  pulse  position  of  13  psec  at  the 
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center  of  the  bubble  and  another  of  approximately  16  psec  from  outside  the  bubble  for  a 
difference  of  3  psec.  This  equates  to  a  height  difference  of  450  pm  which  is  the 
approximate  thickness  difference  between  the  bubble  (1.80  mm)  and  the  undamaged 
glass  fiber  (1.35  mm). 

The  sample  with  the  other  two  bum  spots  was  imaged  next.  First,  the  area  that 
had  been  heated  at  430°C  for  6  minutes  was  imaged  and  is  shown  in  Figure  57a.  Instead 
of  using  the  peak  pulse  amplitude,  the  minimum  pulse  amplitude  was  used  since  it 
showed  the  residue  coating  better.  This  might  have  been  due  to  the  THz  pulse  reflecting 
off  both  the  exterior  paint  and  the  surface  of  the  glass  fiber,  which  could  have  flipped  the 
polarity  of  the  pulse.  The  dot  in  the  center  of  the  circled  bum  spot  corresponds  with  the 
residue.  Next,  the  area  that  was  heated  at  425°C  for  20  minutes  was  imaged  with  the 
reflection  setup  and  is  shown  in  Figure  57b.  Like  transmission  mode,  there  was  no  hint 
of  damage  to  the  glass  fiber,  however,  the  white  marker  dots  showed  up  again  quite  well. 
After  adjusting  the  contrast,  the  apparent  resolution  of  the  dots  is  almost  as  precise  as  the 
transmission  setup  images. 
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Figure  57.  THz  TDS  reflection  images  showing  the  (a)  minimum  peak  amplitude  for 
burning  at  430°C  for  6  minutes  and  the  (b)  peak  pulse  amplitude  for  burning  at 
425°C  for  20  minutes. 
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4.5.1.2  Imaging  of  Bending  Damage 


Next,  the  sample  having  the  mechanical  damage  caused  by  bending  the  glass  fiber 
6240  cycles  was  imaged  in  reflection  setup.  The  image  is  shown  in  Figure  58.  There  is 
perhaps  a  slightly  greater  density  of  lower  amplitude  (1.2  -  1.4  THz)  radiation  around  the 
bend  axis,  but  it  is  difficult  to  categorize  as  conclusive.  Further  attempts  were  made  to 
measure  the  thickness  using  a  longer  40  psec  scan  time,  but  no  additional  detail  was 
obtained. 
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Figure  58.  THz  TDS  reflection  image  using  the  amplitude  of  frequencies  (1.2  -  1.4 
THz)  showing  bend  damage  across  the  central  bend  axis. 

4.5.2.3  Imaging  of  Hidden  Defects 

The  laminated  glass  fiber  coupon  with  hidden  defects  was  imaged  next.  The  3 
mm  diameter  circular  void  was  scanned  first  and  shown  in  Figure  59.  Since  the  void  was 
hidden  below  the  surface,  a  technique  for  imaging  using  the  second  THz  pulse  after  the 
main  pulse  was  used  instead.  This  was  in  an  effort  to  measure  the  first  Fabry-Perot 
reflection  located  at  the  position  of  the  adhesive  and/or  the  air  void.  In  this  case  the  air 
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void  had  a  larger  index  of  refraction  difference  with  the  air  than  the  glue,  allowing  it  to 
appear  in  the  image.  Attempts  were  made  to  show  an  image  of  the  linear  slit  void,  but 
none  of  the  techniques  were  successful  in  isolating  its  location. 
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Figure  59.  THz  TDS  reflection  image  showing  3  mm  diameter  milled  area  hidden 
between  two  glass  fiber  strips  using  the  amplitude  of  the  Fabry-Perot  reflection 
from  the  void. 

4.5.2.4  Imaging  and  Analysis  of  Entire  Panels 

Up  to  this  point,  all  of  the  THz  imaging  has  been  on  thin  coupons  of  glass  fiber 
with  various  forms  of  damage.  On  an  aircraft,  the  glass  fiber  strips  are  attached  to  both 
sides  of  a  much  thicker  honeycomb  structure  and  attached  to  the  external  part  of  the 
aircraft.  A  THz  TDS  transmission  scan  was  attempted  on  an  entire  panel  (~3  cm  thick), 
however,  very  little  of  the  energy  penetrated  through.  Obviously,  TDS  reflection  mode 
would  be  required  to  detect  damage  in  a  panel.  A  couple  of  panels  were  prepared  with 
bum  damage,  holes,  and  paint  remover  on  their  surface  for  imaging  and  analysis  with 
THz  radiation. 

The  first  panel  that  was  scanned  had  holes  (4.5  mm  diameter)  punched  through 
the  glass  fiber  representing  puncture  damage  and  it  also  had  bum  damage  between  the 
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holes.  The  scanned  area  of  the  panel  included  one  hole  and  two  bum  areas.  A  very 
noticeable  burn  area  was  on  the  right  half.  The  other  bum  area  was  much  less  visibly 
noticeable;  in  fact,  I  didn’t  visibly  notice  it  until  it  showed  up  on  the  THz  TDS  scan. 
Images  of  the  panel  can  be  seen  in  Figure  60,  where  the  (a)  minimum  peak  amplitude,  (b) 
maximum  peak  position,  and  (c)  amplitude  spectrum  (1.4  -  1.6  THz)  are  used  to 
construct  the  images.  Finally,  a  plot  of  the  time  domain  pulse  reflected  from  an 
undamaged,  coated  surface  is  compared  to  a  pulse  reflected  from  a  partially  uncoated, 
blistered  portion  of  the  panel  in  60d. 


Figure  60.  THz  TDS  reflection  setup  scan  on  entire  panel  with  puncture  hole  (4.5 
mm  diameter)  and  two  burn  blisters.  Images  were  constructed  using  the  (a) 
minimum  peak  pulse  amplitude,  (b)  maximum  peak  pulse  position,  and  (c)  the 
frequency  spectrum  amplitude  added  under  the  curve  (1.4  -  1.6  THz).  THz  pulses 
reflected  from  coated  surface  and  partially  uncoated  burn  blister  (d). 
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The  next  panel  that  was  scanned  had  a  solvent,  called  butanone,  applied  to  the 
exterior  surface  of  the  panel  to  dissolve  away  the  paint.  The  area  that  was  dissolved  was 
about  7  mm  wide  and  can  be  seen  on  either  of  the  images  in  Figure  61.  One  could  also 
estimate  the  thickness  of  the  paint  that  was  dissolved  away  if  the  index  of  refraction  for 
the  paint  coating  was  known. 
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Figure  61.  THz  images  constructed  using  the  (a)  peak  pulse  amplitude  and  (b)  peak 
pulse  position  of  a  panel  with  its  coating  dissolved  away  by  butanone. 

The  two  past  panels  show  that  THz  TDS  setup  in  reflection  mode  is  capable  of 
imaging  changes  to  the  exterior  of  glass  fiber  on  a  panel.  For  this  research  there  were  no 
samples  available  with  known,  hidden  defects,  but  one  could  assume  that  hidden  defects 
would  reflect  THz  energy  in  the  same  way  that  the  thin  glass  fiber  coupons  would. 
However,  another  important  question  was  how  the  honeycomb  structure  under  the  glass 
fiber  reflects  the  THz  signal.  Long,  time  duration  scans  were  taken  of  an  entire  panel  to 
see  whether  there  were  any  reflections  within  the  honeycomb.  A  typical  plot  of  the  first 
half  of  a  long  scan  is  shown  in  Figure  62a  where  the  time  domain  reflections  off  of  a 
coated  and  uncoated  glass  fiber  surface  are  shown.  Apart  from  the  ZnTe  etalon  effect, 
there  is  nothing  remarkable  about  the  pulse  characteristics  even  when  it  was  extended  to 
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200  psec.  Apparently,  the  honeycomb  either  traps  the  energy  or  it  gets  attenuated 
sufficiently  so  that  it  is  not  detected.  There  is  also  no  back  reflection  off  of  the  back  side 
of  the  glass  fiber,  indicating  a  smooth  transition  from  the  glass  fiber  to  the  honeycomb. 

A  graph  of  the  corresponding  frequency  spectra  are  shown  in  Figure  62b. 


Figure  62.  THz  TDS  (a)  time  domain  plots  and  (b)  frequency  spectra  for  a  coated 
and  uncoated  surface  of  an  aircraft  panel. 

4.5.2.5  Paint  Thickness  Approximations 

An  observation  that  was  noticed  was  that  the  coating  on  the  glass  fiber  caused  a 
double  reflection  of  the  THz  pulse  on  the  front  surface.  A  time  domain  plot  of  this  can  be 
seen  in  Figure  63a  where  there  is  a  double  peak  in  the  positive  polarity  of  the  pulse.  It  is 
also  apparent  that  the  back  reflection  of  the  coated  glass  fiber  sample  is  non-existent, 
making  any  attempt  to  measure  the  material  properties  with  the  back  pulse  to  be 
impossible.  A  plot  of  the  THz  frequency  spectrum  is  shown  in  63b,  where  the  spectrum 
shows  a  fast  roll-off  after  0.4  THz  and  1.6  THz  due  to  the  Fabry-Perot  effect  with  the 
paint.  One  could  use  the  time  domain  plot  to  determine  the  thickness  of  the  coating  given 
the  index  of  refraction  of  the  coating  or  vice  versa.  Since  the  index  of  refraction  for  the 
coating  was  estimated  to  be  around  2.05,  the  thickness  of  the  coating  on  the  glass  fiber 
samples  was  55  pm. 
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Figure  63.  THz  TDS  reflection  configuration  plots  showing  the  (a)  time  domain 
pulse  with  a  double  peak  and  the  (b)  THz  amplitude  spectrum. 

4.5.3  Depth  of  Delaminations  Analysis 


Similar  to  the  transmission  configuration,  an  analysis  was  performed  on  data 
collected  on  the  depth  of  delaminations  for  reflection  configuration.  During  the  data 
collections  it  was  quickly  apparent  that  the  adhesive  -  glass  fiber  interface  provided  a 
large  enough  index  discontinuity  to  create  a  back  reflection.  The  air  -  glass  fiber 
interface  was  even  better  for  reflections  because  the  greater  An  resulted  in  a  higher 
reflection  coefficient.  A  diagram  of  the  reflection  configuration  is  shown  again  in  Figure 
64a.  To  show  the  relative  strength  of  the  first  Fabry-Perot  etalon  reflection  to  the  surface 
reflection,  we  used  the  following  equation: 
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where  the  Fresnel  coefficients  were  shown  previously  in  Equations  83  and  84.  The  graph 
in  Figure  64b  shows  the  relative  strength  of  the  first  Fabry-Perot  reflection  at  various 
depths  off  of  an  air  -  glass  fiber  discontinuity  to  the  strength  of  the  front  pulse.  A  similar 
chart  is  shown  for  a  glue  -  glass  fiber  interface  in  Figure  64c,  including  an  air  gap 
delamination  at  0.93  mm  and  the  back  side  interface  at  1 .25  mm.  These  relative 
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FP  Reflection/  Sample  Pulse 


amplitudes  of  the  Fabry-Perot  reflections  compared  favorably  with  the  amplitudes  of  the 
measured  data  from  the  laminate  sample. 


(a) 


(b)  (c) 

Figure  64.  (a)  Diagram  of  a  delaminated  sample  showing  the  first  Fabry-Perot 
reflection  in  reflection  configuration,  (b)  A  chart  showing  the  relative  strength  of 
the  first  Fabry-Perot  reflection  after  traveling  through  various  thicknesses  of  glass 
fiber  composite  before  the  delamination,  (c)  A  chart  showing  the  relative  strength 
of  the  first  Fabry-Perot  reflection  from  discontinuities  including  both  glue-glass 
fiber  (G-GF)  and  air-glass  fiber  (A-GF). 

The  next  step  was  to  analyze  data  that  was  collected  at  the  various  depths. 
Examples  of  TDS  time  domain  scans  in  reflection  mode  are  shown  in  Figure  65.  The 
scan  in  65a  shows  the  very  distinct  reflections  that  occur  when  the  THz  pulse  reflects 
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from  the  delaminated  air-glass  fiber  interface  and  also  the  back  air  -  glass  fiber  interface. 
The  scan  in  65b  shows  an  example  of  a  Fabry-Perot  reflection  from  a  glue  -  glass  fiber 
interface.  Both  discontinuities  were  located  at  0.23  mm  below  the  sample  surface  and  the 
back  side  reflection  can  also  be  seen  in  both. 


Tune  (psec) 


Figure  65.  THz  TDS  time  domain  plots  showing  reflections  from  discontinuities: 
(a)  air  and  (b)  adhesive. 

The  technique  that  was  used  for  depth  determination  was  to  use  a  Fourier 
deconvolution  D(t)  in  the  time  domain  instead  of  the  autocorrelation  used  for 
transmission  mode.  This  is  demonstrated  in  the  following  equation: 


D(t)= 


(87) 


in  which  fft  is  the  fast  Fourier  transform  and  ifft  is  the  inverse  fast  Fourier  transform.  A 
comparison  of  the  returns  from  discontinuities  at  various  depths  is  shown  in  Figure  66  a- 
f.  The  dashed  lines  were  added  to  show  the  similarity  of  the  backside  reflections  in  each 
scan  along  with  the  pulse  reflection  from  the  ZnTe  etalon  effect.  A  separate  calculation 
was  made  to  remove  the  TDS  system  effects  from  the  deconvolution.  This  was 
accomplished  by  making  a  separate  measurement  of  the  THz  pulse  reflection  from 
aluminum  and  then  subtracting  out  the  aluminum  from  the  sample.  This  removed  most 
of  the  main  pulse  and  the  ZnTe  reflection.  A  series  of  depth  measurements  can  be  seen  in 
Figure  66  g-1  for  the  various  thicknesses.  Each  series  of  plots  was  measured  from  a 
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Backside  ZnTe 


Backside 


Figure  66.  Fourier  deconvolution  of  the  discontinuities  present  in  a  glass  fiber 
sample  laminated  at  various  thicknesses:  (a)  0.23  mm,  (b)  0.46  mm,  (c)  0.69  mm,  (d) 
0.92  mm,  (e)  0.92  mm  delaminated,  and  (f)  no  lamination.  Fourier  deconvolution 
with  the  TDS  system  characteristics  subtracted  are  shown  next.  The  various 
thicknesses  were  measured  in  the  opposite  direction:  (g)  0.23  mm,  (h)  0.23  mm 
delaminated,  (i)  0.46  mm,  (j)  0.69  mm,  (k)  0.92  mm,  and  (1)  no  lamination. 
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different  side  of  the  sample,  showing  minimal  difference  from  each  side.  A  plot  of  the 
calculated  times  and  the  autocorrelation  times  for  the  depth  of  delaminations  measured 
with  transmission  and  reflection  configuration  is  shown  in  Figure  67.  The  measured 
autocorrelation  time  for  transmission  configuration  is  questionable  beyond  a 
depth  of  0.23  mm.  Both  of  the  techniques  are  capable  of  imaging  defects  on  the  order  of 
magnitude  of  a  THz  wavelength.  A  comparison  of  THz  images  between  transmission 
and  reflection  configuration  is  shown  in  Figure  67. 
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Figure  67.  Comparison  of  the  calculated  time  for  F-P  reflections  through  glass  fiber 
to  the  measured  autocorrelation  values  measured  in  both  transmission  and 
reflection  configuration. 
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5.  Conclusions  and  Recommendations 


This  chapter  provides  a  summary  of  the  research  that  was  presented  in  the 
previous  chapters  followed  by  directions  for  future  research.  It  explains  how  THz  TDS 
has  the  potential  to  be  a  superior  NDE  technique  for  glass  fiber  composites  over  other 
NDE  techniques.  It  then  compares  and  contrasts  the  advantages  and  disadvantages  of  the 
TDS  transmission  or  reflection  setup.  Finally,  it  provides  recommendations  for  follow- 
on  research  projects. 

5.1  Major  Research  Contributions 

This  research  was  novel  in  a  variety  of  different  aspects.  First,  it  is  one  of  the  first 
research  efforts  that  uses  THz  TDS  as  an  NDE  technique  to  inspect  glass  fiber  aircraft 
composite  samples  with  representative  forms  of  damage.  Second,  it  demonstrates  how  to 
use  THz  TDS  to  measure  the  material  properties  of  a  bum  damage  area  on  a  composite  in 
an  effort  to  determine  if  the  material  has  been  chemically  altered.  Third,  it  shows  how 
various  types  of  surface  and  subsurface  defects  can  be  imaged  and  detected  by  either  a 
transmissive  or  reflective  TDS  configuration.  Fourth,  it  also  provides  an  analysis  on  how 
effective  each  configuration  would  be  at  isolating  the  depth  of  discontinuities  within  the 
composite.  Finally,  it  also  presents  a  technique  for  measuring  material  parameters  in 
reflection  mode. 

5.2  Conclusions 

THz  TDS  has  several  potential  advantages  over  other  nondestructive  evaluation 
methods  for  inspection  of  aircraft  glass  fiber  composites.  Both  ultrasound  and  eddy 
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current  techniques  require  that  the  source  and  detector  remain  in  contact  with  the  aircraft. 
In  addition,  eddy  current  testing  requires  that  the  material  is  conductive,  which  is  not 
applicable  to  glass  fiber.  X-ray  techniques  use  ionizing  radiation  to  penetrate  the  sample 
which  can  cause  safety  issues  for  operators.  X-rays  can  also  require  the  use  of  penetrates 
to  help  detect  delaminations.  Video  inspection  can  be  time  consuming  and  prone  to 
human  error.  Microwaves  have  the  ability  to  penetrate  composites,  but  their  resolution  is 
more  limited  due  to  their  longer  wavelengths.  Infrared  techniques  can  have  better 
resolution,  but  less  penetration  depth.  THz  can  penetrate  glass  fiber  without  contacting  it, 
with  submillimeter  transverse  resolution,  and  can  detect  surface  defects,  hidden  voids, 
delaminations,  and  bending  damage  in  composites.  The  depth  resolution  for  defects  is 
potentially  much  better  than  the  transverse  resolution,  and  thin  composite  samples  can  be 
analyzed  with  THz  TDS.  Additionally,  it  can  also  be  used  to  evaluate  whether  the 
aircraft  composite  has  been  chemically  altered  from  engine  bum  damage  by  measurement 
of  its  index  of  refraction  and  absorption  coefficient  spectrum. 

Material  parameter  measurements  were  made  with  both  transmission  setup  and 
reflection  setup.  The  dielectric  properties  of  low  resistivity  silicon  wafers  were  modeled 
with  the  Drude  model  and  showed  a  close  comparison  with  measured  results.  Fabry- 
Perot  etalon  removal  techniques  were  moderately  effective  in  removing  the  oscillations 
from  the  index  of  refraction  and  absorption  spectra.  It  was  more  difficult  to  measure 
material  parameters  with  reflection  configuration  because  the  technique  of  comparing  the 
sample  pulse  to  a  reference  pulse  did  not  provide  accurate  results.  For  thicker  samples,  a 
time  domain  separation  of  the  front  and  back  reflections  provides  a  useful  technique,  if 
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one  assumes  a  flat  refractive  index  response.  It  becomes  more  difficult  when  the  samples 
are  thinner  and  their  absorption  coefficients  are  higher. 

Effective  medium  approximations  were  used  to  estimate  the  refractive  index  and 
absorption  coefficient  spectrum  of  the  glass  weave  after  measuring  the  polyimide  and 
glass  fiber  composite  separately.  These  indices  of  refraction  were  then  used  to  calculate 
the  scattering  properties  of  the  composite  and  estimates  of  the  scattering  coefficient  were 
low  for  the  THz  spectrum  due  to  the  similarity  of  the  indices  of  refraction  for  the  two 
composite  materials  and  the  density  of  inclusions. 

An  aircraft  glass  fiber  composite  with  various  forms  of  damage  was  examined 
using  both  a  transmissive  and  a  reflective  THz  TDS  system  for  NDE.  Indices  of 
refraction  and  absorption  coefficients  in  the  terahertz  frequency  range  were  measured 
using  THz  TDS  in  transmission  configuration  for  comparison  of  damaged  and 
undamaged  material  states.  Results  showed  that  localized  heat  damage  did  not  change 
the  material  properties  of  the  composite  sample  noticeably;  however,  changes  in  the 
terahertz  signal  were  observed  due  to  material  blistering  and/or  residue  deposited.  The 
approximate  depth  of  a  delaminations  could  be  determined  in  the  time  domain  by 
measuring  the  timing  of  a  Fabry-Perot  reflection  through  a  thin  slice  of  the  composite. 
Reflection  configuration  was  much  more  effective  in  showing  discontinuities,  even 
between  an  adhesive  layer  and  the  glass  fiber.  A  hidden  circular  void  was  imaged  and 
there  was  also  evidence  that  areas  of  damage  from  mechanical  bending  stress  and 
simulated  hidden  cracks  could  be  detected  with  terahertz  TDS  imaging  in  both 
configurations. 
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THz  TDS  in  transmission  setup  was  effective  in  locating  voids  at  any  depth  or 
thickness,  but  had  difficulty  in  finding  damage  that  was  smaller  than  THz  wavelengths. 
Bum  spots  could  be  detected,  but  the  dielectric  properties  of  polyimide  appeared  to 
remain  unaltered  at  bum  temperatures  below  480°C.  Additional  bum  tests  at 
temperatures  approaching  480°C  would  have  to  be  attempted  to  see  if  and  when  the 
material  properties  of  the  composite  would  change.  A  THz  TDS  system  in  reflection 
mode  would  be  a  more  likely  candidate  for  inspection  of  glass  fiber  composites  on 
aircraft  and  would  likely  be  more  effective  in  determining  the  depth  of  damage  areas. 

The  glass  fiber  composite  absorbs  a  significant  amount  of  the  THz  radiation,  attenuating 
a  transmission  signal,  and  limiting  the  use  of  Fabry-Perot  reflections  for  depth 
measurements.  In  reflection  mode,  the  magnitude  of  the  first  surface  reflection  and 
subsequent  reflections  were  more  similar,  which  enabled  a  greater  ability  to  locate  the 
depth  of  a  void. 

There  were  advantages  and  disadvantages  to  using  either  the  TDS  transmission  or 
reflection  configuration.  It  was  easier  to  get  a  smaller  THz  spot  size  onto  the  sample 
using  a  transmission  configuration.  This  was  primarily  due  to  the  large  ratio  of  focal 
length  to  diameter  in  the  parabolic  mirrors  and  the  cumbersome  stages  used  to  steer  them. 
Custom  optical  parts  could  alleviate  this  problem  for  reflection  configurations. 

Reflection  mode  is  certainly  a  more  likely  candidate  for  a  maintenance  technique,  as  was 
mentioned  earlier,  but  transmission  setup  could  aid  in  the  design  and  quality  control  of 
glass  fiber  composites.  Reflection  mode  would  also  be  a  more  effective  technique  in 
determining  the  depth  of  delaminations.  This  is  due  to  reflection  mode  signals  traveling  a 
shorter  distance  through  the  attenuating  glass  fiber  and  due  to  the  relative  similarity  in  the 
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amplitude  of  Fresnel  reflections  between  a  void  and  the  initial  pulse.  This  can  be  easily 
seen  when  one  compares  the  results  in  Figure  64  to  Figure  54,  which  is  shown  again  as  a 
comparison  in  Figure  68.  The  relative  amplitude  of  the  first  Fabry-Perot  reflection  is 
calculated  to  be  an  order  of  magnitude  larger  in  reflection  configuration  than  in 
transmission. 


(a)  (b) 

Figure  68.  Summary  charts  showing  the  relative  amplitude  of  the  first  Fabry-Perot 
reflection  to  the  main  pulse  after  traveling  through  various  thicknesses  of  glass  fiber 
composite  in  (a)  transmission  configuration  and  (b)  reflection  configuration. 

Finally,  a  comparison  of  the  transmission  and  reflection  configuration  images  are 

shown  next  to  a  photograph  of  the  damage  sample  in  Figure  69.  THz  images  are  shown 

for  the  hidden  circular  void,  the  calibration  sample  using  both  phase  and  amplitude 

techniques,  and  the  three  bum  damage  areas. 


5.3  Future  Work 

One  of  the  areas  of  research  that  is  currently  on-going  is  to  compare  the  results  of 
THz  TDS  on  a  set  of  glass  fiber  composites  samples  to  other  NDE  techniques,  such  as 
ultrasound,  x-ray  tomography,  and  thermography.  This  comparison  could  be  very 
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Transmission 


Reflection 


Figure  69.  Comparison  of  various  types  of  glass  fiber  damage/alterations  between 
transmission  and  reflection  configuration. 
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valuable  in  defining  the  potential  strengths  and  weaknesses  of  using  THz  TDS  for  NDE 
of  glass  fiber  composites.  THz  TDS  could  also  be  investigated  as  an  NDE  technique  for 
two  other  types  of  composites:  Kevlar  and  carbon  fiber.  The  ability  of  THz  to  penetrate 
carbon  fiber  would  probably  be  limited  due  to  its  higher  conductivity,  but  it  could  be 
useful  to  determine  the  penetration  depth  and  surface  defects.  Kevlar  and  glass  fiber 
samples,  with  various  forms  of  damage,  are  also  available  from  AFRL/RXLP. 

Another  task  that  could  be  performed  in  future  research  would  be  to  calibrate  the 
absorption  coefficient  measurements  with  the  THz  TDS  transmission  setup.  The 
measurements  on  the  absorption  coefficient  in  this  dissertation  appeared  to  be 
accurate,  but  did  not  explore  coefficients  much  larger  than  0  cm'1.  The  silicon  wafer  or 
other  semiconductor  could  be  measured  with  THz  TDS  and  compared  to  the  Drude 
Model  prediction.  For  this  dissertation,  the  only  types  of  wafers  that  were  obtained  in 
single  quantities  were  those  with  a  high  resistivity.  An  attempt  was  made  to  find  a 
semiconductor  with  a  higher  doping  level  that  could  be  purchased  in  small  quantities 
economically,  but  none  could  be  located. 

An  area  that  could  be  explored  in  more  detail  would  be  the  modeling  of  the  THz 
pulse  shape.  In  the  laboratory,  the  pulse  shape  characteristics  could  be  changed 
drastically  by  altering  the  position  of  any  of  the  optics,  including  the  parabolic  mirrors 
and  hemispherical  lenses,  and  the  alignment  of  the  pump  or  probe  beam  on  either  the 
emitter  or  detector.  One  could  characterize  the  effects  that  each  of  the  parts  of  the  THz 
quasi-optical  system  has  on  the  final  pulse  shape  and  frequency  spectrum.  One  would 
have  to  be  able  to  isolate  each  component,  vary  the  alignment,  and  measure  the  output. 
These  measurements  could  then  be  compared  to  a  model  of  the  THz  Gaussian  pulse 
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shape  in  space  and  time.  Perhaps  the  most  difficult  items  to  characterize  would  be  the 
emitter  and  detector,  and  a  special  arrangement  with  their  manufacturer  may  be  required. 

The  analysis  of  the  effective  medium  approximations  could  be  expanded.  A  solid 
piece  of  the  glass  used  in  the  composite  would  have  to  be  obtained  so  that  it  could  be 
accurately  measured  with  the  THz  TDS  transmission  configuration.  Then  various 
volume  concentrations  of  the  glass  and  polyimide  could  be  mixed  to  form  different 
composite  mixtures.  The  material  properties  of  each  composite  mixture  could  then  be 
measured  and  compared  directly  to  the  effective  medium  approximations. 

One  of  the  quantitative  NDE  analyses  used  to  assess  the  value  of  the  technique  is 
to  measure  its  resolution  in  detecting  and  identifying  both  surface  and  hidden  damage. 
Another  area  of  future  research  would  be  to  measure  damage  of  various  sizes,  shapes,  and 
depths  to  define  the  limits  of  either  a  THz  TDS  transmission  or  reflection  setup.  The 
THz  spot  size  would  have  to  be  minimized  in  either  case,  using  short  focal  length  optics 
and  small  translation  stages.  Samples  would  have  to  be  prepared  with  defects  that  have 
been  carefully  manufactured  and  measured,  and  with  those  defects  varying  in  size  from 
the  largest  THz  wavelengths  down  to  a  fraction  of  the  smallest  THz  wavelengths. 

Finally,  the  last  area  of  suggested  future  research  would  be  to  design  a  setup  in 
reflection  mode  which  could  accurately  measure  the  index  of  refraction  and  absorption 
coefficient.  The  reflection  technique  that  was  used  in  this  research  was  not  robust  or 
precise  enough  to  consistently  measure  material  parameters  accurately.  One  would  have 
to  use  an  alignment  system  to  ensure  that  the  sample  and  reference  were  within  1  pm 
across  their  entire  reflected  surface,  if  the  setup  in  Figure  29  was  used.  An  alternate 
technique  would  be  to  use  ellipsometry  to  measure  the  material  properties  in  which  no 
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reference  measurement  is  needed.  One  would  place  a  wire  grid  polarizer  after  the  emitter 
to  split  the  THz  into  s-  and  p-  polarizations  followed  by  a  wire  grid  analyzer  before  the 
detector  to  detect  the  two  polarizations  with  equal  sensitivity.  The  sample  would  then  be 
placed  at  an  angle  near  the  Brewster  angle  for  a  good  signal-to-noise  ratio,  the 
ellipsometric  angles  measured  and  fitted,  and  the  complex  refractive  index  determined. 
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